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The  primary  objective  of  the  Test  Bed  program  was  to  bridge  the  gap 
between  the  fundamental  studjes  and  the  application  efforts  relating  to 
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physical  Test  Bed,  as  well  as  th£  fabrication  of  the  ultrasonic  phased 
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quantitative  NPE  techniques.  The  .techniques  that  Mere  Implemented 
consisted  of  a  ;$ubset  of  the  deft^fpclasslfl cation  and  Imaging  tech- 
nlques  that  were  developed  In  the  Interdisciplinary  Program  for  Quan¬ 
titative  Flaw  Definition.  The  third  phase  consisted  of  a  detailed 
system  evaluation  on  real  and  simulated  problems.^ 

The  work  with  the  quantitative  Inversion  techniques  demonstrated 
the  feasibility  of  combining  data  obtained  from  lowy  medium  and  high 
frequency  portions  of  the  flaw  spectrin  to  obtalnan  accurate  estimate 
of  the  size,  shape  orientation  and  composition  or  the  flaw.  The  Bom 
Inversion  technique  for  classification  of  defeats  was  studied  exten¬ 
sively.  The  effects  on  the  accuracy  of  theBorn  Inversion  due  to 
signal -to-nolse  ratio,  bandwidth,  andgebmetry  were  all  Investigated. 
This  work  culminated  In  a  successful  demonstration  with  a  naturally 
occurring  flaw  In  a  nickel  based  alloy  supplied  by  Rolls-Royce  Ltd. 

The  phased  array  system  utilizes  two,  32  element  transducer  arrays 
operating  at  2.25  MHz.  It  features  parallel  processing  of  the  16  re¬ 
ceiving  channels,  qn  array  processor  for  “beam- forming",  and  mini-com¬ 
puter  control  of  the  system.  A  custom  hybrid  pul ser/ receiver  was  de¬ 
signed  and  fabricated  specifically  for  this  system.  The  system  has  been 
used  for  Imaging  known  flaws  In  metal  parts  and  for  obtaining  scatter¬ 
ing  data  from  these  flaws. 

^The  Test  Bed  program  has  addressed  a  wide  range  of  topics  that  are 
Important  to  a  quantitative  NDE  program.  It  has  successfully  demon¬ 
strated  that  quantitative  NDE  techniques  can  be  applied  to  naturally 
occurring  flaws.  It  provides  a  sound  technological  base  for  the  exten¬ 
sion  of  these  techniques  to  practical  NDE  applications.^ 
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1.0  INTRODUCTION 


1.1  Objective 

The  primary  objective  of  this  program  was  to  bridge  the  gap  between 
the  fundamental  studies  and  the  application  efforts  relating  to  ultrasonic 
NDE.  To  accomplish  this  objective,  the  program  created  a  test  facility 
containing  state-of-the-art  ultrasonics  equipmet  that  is  sufficiently 
flexible  to  allow  the  Implementation  and  evaluation  of  new  techniques  and 
equipment.  This  evaluation  has  been  done  using  both  simulated  and  real 
defects  which  occur  in  laboratory  and  real  aircraft  samples.  The  facility  is 
accessible  to  researchers  not  associated  with  the  contractor  to  facilitate  the 
development  of  new  NDE  methods  and  techniques. 

Interested  parties  may  use  the  facility  at  a  nominal  cost.  For 
information  regarding  procedures  for  accessing  the  facility  and  business 
arrangements,  contact: 

P.H.  Milham 

Director  of  Contracts  and  Pricing 
Rockwell  International  Science  Center 
1049  Camlno  Dos  Rios 
Thousand  Oaks,  CA  91360 
PH  (805)  498-4545 

Inquiries  concerning  technical  capabilities  should  be  directed  to  R.C. 

Addison,  Jr.,  at  the  same  address. 

The  program  was  composed  of  three  phases.  The  first  phase  consisted 
of  the  system  design,  assembly  and  shakedown.  This  phase  encompassed  the 
construction  of  the  self  contained  unit  comprising  the  physical  Test  Bed,  as 
well  as  the  fabrication  of  An  ultrasonic  phased  array  system.  The  second 
phase  Involved  the  reduction  to  practice  of  quantitative  NDE  techniques.  The 
techniques  that  were  Implemented  consisted  of  a  subset  of  the  long  wavelength. 


defect  classification  and  Imaging  techniques  that  were  developed  In  the 
Interdisciplinary  Program  for  Quantitative  Flaw  Definition. *  The  third  phase 
consisted  of  a  detailed  system  evaluation  on  real  and  simulated  problems. 

1.2  Conceptual  Design 

The  Ultrasonic  Test  Bed  program  was  Initiated  to  complement  the 
AF/DARPA  Interdisciplinary  Program  for  Quantitative  Flaw  Definition.1  This 
Interdisciplinary  program  has  fostered  the  development  of  a  variety  of  new 
techniques  for  the  quantitative  characterization  of  flaws.  The  program  was 
and  Is  currently  concerned  with  characterization  methods  encompassing  both 
ultrasonic  and  electromagnetic  techniques.  A  major  thrust  of  the  program  has 
been  directed  towards  the  development  of  ultrasonic  inversion  techniques. 

These  techniques  permit  the  characterization  of  a  flaw  based  on  the 
measurement  of  the  ultrasonic  fields  Incident  on  and  scattered  by  the  flaw. 

The  utility  of  these  techniques  has  been  categorized  according  to  the  ratio  of 
a  representative  flaw  dimension  to  the  wavelength  of  the  ultrasound.  For 
convenience  this  ratio  Is  referred  to  In  terms  of  the  dimensionless  quantity 
ka,  where  k  *  and  a  Is  a  representative  radius  of  the  flaw.  Thus  the 

Inversion  technique  commonly  referred  to  as  Imaging  Is  useful  at  large  ka 
values  In  the  regime  where  the  flaw  size  is  larger  than  the  wavelength.  In 
the  medium  ka  range  where  0.5  <  ka  <  6,  the  flaw  size  Is  less  than  or 
comparable  to  the  wavelength.  Techniques  applicable  In  this  ka  range  are 
herein  referred  to  as  model-based  reconstruction  techniques.  The  low  ka  range 
refers  to  values  of  ka  less  than  0.5.  The  inversion  methods  used  In  this 
region  are  simply  called  long  wavelength  techniques. 

At  the  Inception  of  the  Ultrasonic  Test  Bed  program  the  variety  of 
techniques  that  had  been  developed  to  operate  In  each  of  these  three  regimes 
had  been  shown  to  work  theoretically  and  In  some  cases  experiments  with 
Idealized  geometries  had  been  conducted.  What  was  needed  was  some  means  to 
test  these  Inversion  algorithms  under  conditions  that  were  representative  of 
those  found  In  real  parts  and  to  Integrate  those  algorithms  that  were  suitable 
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into  a  coherent  system  that  could  collect  data  over  the  entire  range  of  ka 
values. 

The  Ultrasonic  Test  Bed  program  was  designed  to  carry  out  these 
tasks.  The  Test  Bed  program  was  specifically  concerned  with  the  ultrasonic 
techniques  and  of  those  available  It  has  concentrated  on  the  subsets  that  are 
useful  in  charcterlzlng  internal  flaws  such  as  voids  and  Inclusions.  It  is 
Important  to  note  that  the  procedure  for  testing  these  techniques  proved  to  be 
a  lengthier  task  than  had  been  initially  estimated.  This  resulted  in  a  selec¬ 
tion  of  the  available  techniques  to  produce  a  flaw  analysis  capability  that 
extended  over  all  wavelength  ranges  and  yet  permitted  the  program  to  be  com¬ 
pleted  within  the  available  time  frame.  As  an  example,  a  long  wavelength 
technique  for  accurately  locating  the  center  of  a  flaw  was  Implemented  because 
it  was  useful  In  improving  the  accuracy  of  the  model -based  reconstruction 
technique  known  as  the  inverse  Born  approximation.  However  the  long  wave¬ 
length  technique  for  measuring  the  flaw  related  parameter  known  as  A2  was  not 
Implemented. 

The  work  statement  specifically  mentioned  that  the  defect  classifica¬ 
tion  technique  that  utilized  nonlinear  networks  was  to  be  implemented.  This 
technique  was  not  Included  among  those  that  were  tested.  A  preliminary  inves¬ 
tigation  of  this  technique  revealed  that  the  applicable  ka  range  for  these 
networks  is  1  <  ka  <  3.  The  range  of  applicability  of  the  inverse  Born 
approximation  completely  overlaps  this  range.  Further,  the  nonlinear  networks 
require  measurements  of  the  flaw  to  be  made  from  a  minimum  of  12  different 
angles  with  19  Independent  measurements  being  preferred.  This  is  an  Imprac- 
tlcally  large  number  of  measurements  to  characterize  a  single  flaw.  In  addi¬ 
tion  to  the  time  required  to  acquire  the  data,  the  geometry  of  most  real  parts 
would  restrict  the  spatial  window  to  a  relatively  few  measurement  directions. 
Despite  these  drawbacks  there  was  still  Interest  In  testing  the  utility  of  the 
networks.  The  software  for  using  the  non-linear  networks  did  not  arrive  until 
the  program  for  analyzing  flaws  as  well  advanced  and  much  of  the  available 
time  had  already  been  allocated.  When  the  software  was  received  it  was  not 
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compatible  with  the  computer  used  for  the  Test  Bed.  A  significant  effort 
would  have  been  required  to  correct  the  software  so  that  It  would  run  on  the 
Test  Bed  computer.  Given  the  other  drawbacks  of  this  particular  technique,  a 
decision  was  made  to  omit  this  technique  from  consideration. 

A  useful  function  of  the  Test  Bed  program  has  been  to  communicate 
Information  back  to  the  Interdisciplinary  program  concerning  areas  where 
practical  constraints  degraded  the  results  and  further  fundamental  work  was 
required. 

To  achieve  these  goals,  the  Test  Bed  was  built  as  a  self-contained 
unit,  complete  with  a  large  ultrasonic  water  tank,  transducer  array,  mechan¬ 
ical  positioning  devices,  analog  and  digital  signal  generation  and  processing 
circuitry,  dedicated  mnlcomputer  for  control  and  data  Interpretation,  and 
display  units.  Fundamental  to  the  system  Is  the  capability  to  scan  tran¬ 
sducers  over  parts  and  acquire  data.  Important  differences  between  the  Test 
Bed  and  many  of  the  conventional  systems  are  that  the  multi  axis  scanning 
system  Is  controlled  by  a  microprocessor  system  and  Is  capable  of  contour 
following  modes  of  scanning.  The  data  acquisition  techniques  differ  In  that 
entire  waveforms  are  digitized  and  stored  on  the  disk  memory  of  the  mini¬ 
computer.  The  chief  distinction  of  the  Test  Bed  system,  however.  Is  the 
presence  of  algorithms  within  the  minicomputer  for  carrying  out  the  inversion 
techniques  outlined  above.  The  capabilities  of  the  Test  Bed  for  Imaging 
studies  are  enhanced  by  an  ultrasonic  phased  array  system  that  was  constructed 
at  the  Science  Center  as  part  of  the  program. 

1.3  Structure  of  Report 

This  report  covering  the  activities  and  accomplishments  of  the  Test 
Bed  program  Is  divided  Into  three  major  sections.  The  first  of  these,  the 
Test  Bed  System,  provides  a  description  of  the  various  hardware  subsystems 
that  are  Interconnected  to  form  the  physical  system.  The  major  hardware  Items 
covered  are  the  scanning  system,  the  data  acquisition  device,  the  display 
processor,  and  the  minicomputer  system.  The  software  for  controlling  the 
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operation  of  this  system  Is  also  detailed.  This  software  Includes  the  main 
multitasking  program,  as  well  as  the  specific  programs  for  directing  the 
multi  axis  controller  and  for  displaying  data  via  the  display  processor. 

The  next  section  describes  the  inspection  procedures.  A  thorough 
explanation  Is  given  of  the  protocol  used  in  examining  a  part.  The  operation 
of  the  one-dimensional  Born  inversion  algorithm  is  explained  and  the 
procedures  used  in  testing  the  algorithm  are  presented  in  conjunction  with  the 
results  of  these  tests.  Finally  experimental  results  with  known  flaws  and 
unknown  flaws  In  various  specimens  are  provided. 

The  final  major  section  details  the  design,  fabrication  and  testing 
of  the  ultrasonic  phased  array  system.  All  of  the  major  subsystems  are 
described  from  the  transducers  to  the  array  processor.  The  general  aspects  of 
phased  array  operation  and  the  fundamental  advantages  and  limitations  are 
discussed.  Results  obtained  with  the  system  for  flaws  within  metal  parts  are 
presented,  for  both  Inversion  methods  Involving  imaging  and  Born  inversion 
techniques. 

The  report  ends  with  a  section  providing  conclusions  and  recom¬ 
mendations  for  future  work  In  the  area  encompassed  by  the  Test  Bed  program. 
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2.0  TEST  BED  SYSTEM 


2.1  Overview 

The  Test  Bed  Is  a  multifaceted  system  with  many  of  Its  unique 
capabilities  residing  In  especially  written  algorithms  and  signal  processing 
techniques.  There  Is,  however,  a  significant  amount  of  hardware  that  com¬ 
prises  the  physical  portion  of  the  system.  The  photograph  of  the  Test  Bed 
laboratory  In  Fig.  2.1  shows  most  of  this  hardware  and  the  block  diagram  In 
Fig.  2.2  shows  how  It  Is  interconnected.  The  major  component  of  the  system  is 
the  six-axis  scanning  system  and  the  large  water  tank  that  dominate  the 
foreground.  Closely  associated  with  this  Is  the  microprocessor  based  multi¬ 
axis  controller  shown  to  the  right  of  the  water  tank.  In  the  background,  the 
Data  General  ECLIPSE  S/200  minicomputer  can  be  seen.  The  color  raster-scan 
image  display  unit  (TV  monitor)  associated  with  the  Genlsco  display  processor 
Is  shown  to  the  left  of  the  S/200.  Adjacent  to  It  Is  one  of  the  Tektronix 
4006  terminals  associated  with  the  system.  Behind  the  terminal  can  be  seen 
the  Data  General  disk  memory  unit.  Not  visible  In  the  photograph  are  the 
Blomatlon  8100  A/D  converter  used  for  data  acquisition,  the  pulser/  receiver, 
and  the  Tektronix  453A  oscilloscope  used  with  the  system.  There  are  also  two 
hard  copy  units  available,  a  Versatec  unit  for  the  computer  terminals  and  an 
Image  Resources  Videoprint  500  for  the  Image  display  unit. 

This  hardware  provides  a  capability  of  scanning  rotational ly  sym¬ 
metric  parts  using  a  contour  scanning  technique  as  well  as  flat  parts  using  a 
raster  scanning  technique.  The  pul ser/recel ver  Is  a  unit  made  by  Panametrlcs 
that  satisfies  the  needs  of  this  system.  The  received  waveforms  are  digitized 
with  a  Blomatlon  data  acquisition  unit  and  stored  on  the  disk  memory  of  the 
Oata  General  minicomputer.  With  the  possible  exception  of  the  waveform 
digitization,  these  tasks  are  accomplished  using  procedures  that  are 
widespread  and  to  a  large  extent  can  be  considered  as  conventional. 

The  system  has  the  capability  of  displaying  the  conventional  C-scan 
and  B-scan  Images  using  a  variety  of  amplitude  coding  schemes  Involving  the 
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Fig.  2.2  Block  diagram  of  Test  Bed  system. 


gray  scale  or  pseudo-color  rendition  of  the  data.  These  techniques  are  not 
generally  available  in  conventional  Inspection  systems  but  significantly 
enhance  the  information  content  conveyed  by  the  Images. 


The  unique  capabilities  of  the  system  reside  in  its  ability  to 
manipulate  the  data  with  the  minicomputer.  This  permits  an  array  of  signal  - 
processing  techniques  Including  time  gating,  frequency  filtering,  and  trans¬ 
ducer  deconvolution.  These  techniques  are  frequently  used  to  prepare  signals 
for  processing  by  algorithms  designed  to  characterize  detected  flaws  In  terms 
of  their  size,  shape,  orientation,  and  composition. 

2.2  Scanning  Sy:tem 

2.2.1  Hardware 

The  scanning  system  for  the  Test  Bed  provides  five  degrees  of  freedom 
for  the  transducer  plus  a  turntable  that  can  be  used  to  rotate  circularly 
symmetric  parts.  The  X  and  Y  motions  are  based  on  an  assembly  that  existed  at 
the  Science  Center  prior  to  the  Test  Bed  program.  This  assembly  was  upgraded 
so  that  the  bridge  could  support  a  manipulator  arm  and  the  frame  work  would 
fit  around  the  water  tank.  Although  no  special  efforts  were  made  to  assure 
precise  mechanical  positioning  accuracies,  the  rigidity  Is  probably  adequate 
to  Insure  an  accuracy  of  about  0.040  Inches  over  the  entire  range  of  motion. 
The  X  and  Y  axes  are  driven  by  stepping  motors  having  200  steps/revolution  via 
ball  screws  that  provide  a  linear  precision  of  0.001  Inches/step.  The 
relative  accuracy  over  short  distances  of  less  than  10  Inches  Is  about  0.010 
Inches.  The  maximum  speed  of  the  X  or  Y  axis  Is  2  Inches  per  second  If  the 
motor  Is  ramped  up  to  the  maximum  speed.  Without  ramping  the  maximum  speed  Is 
0.3  Inches  per  second. 

The  Z  axis  motion  and  the  two  glmbal  axes  are  provided  by  a  manipula¬ 
tor  that  Is  available  commercially  from  Automation  Industries.  The  Z  axis  has 
a  maximum  travel  distance  of  48  Inches.  It  Is  driven  by  a  stepping  motor  that 
provides  a  linear  precision  of  0.001  Inches/step.  The  accuracy  for  the  full 


travel  range  Is  specified  to  be  0.010  Inches.  The  angular  motions  are  pro¬ 
vided  by  a  primary  and  secondary  gimbal  that  are  driven  by  stepping  motors. 

The  primary  gimbal  rotates  about  an  axis  that  is  parallel  to  the  Y  axis  of  the 
system  and  Is  designated  the  B  axis.  Similarly  the  secondary  gimbal  rotates 
about  an  axis  that  Is  orthogonal  to  the  B  axis  and  is  designated  the  A  axis. 
The  primary  gimbal  can  travel  over  a  total  angular  range  of  300°  with  a  pre¬ 
cision  of  0.1°  per  step.  The  secondary  gimbal  can  travel  over  a  total  angular 
range  of  ±  40°  with  the  same  precision.  The  accuracy  of  these  axes  has  not 
been  determined  but  appears  to  be  significantly  worse  than  the  precision  and 
Is  affected  by  an  Interaction  between  the  axes.  The  inaccuracies  in  these 
axes  are  the  chief  sources  of  error  in  the  overall  pointing  accuracy  of  the 
system. 

The  sixth  degree  of  freedom  is  provided  by  a  turntable  that  is  avail¬ 
able  commercially  from  Automation  Industries.  The  turntable  has  a  three-jaw 
chuck  capable  of  holding  parts  with  diameters  up  to  30  inches.  The  unit  Is 
driven  by  a  stepping  motor  with  a  precision  of  0.1°  per  step  and  a  maximum 
speed  of  15  rpm. 

The  scanning  assembly  Is  built  around  a  water  tank  with  overall 
dimensions  of  5  feet  long  x  3.25  feet  wide  *  3.5  feet  deep.  The  tank  Is  made 
of  mild  steel  with  a  PVC  liner.  There  Is  a  recirculating  pump  and  filters  for 
keeping  the  water  clean. 

2.2.2  Multi axis  Controller 

The  multiaxis  controller  Is  a  microprocessor  based  unit  that  controls 
the  movements  of  each  of  the  six  axes.  A  block  diagram  of  the  unit  Is  shown 
In  Pig.  2.3.  The  system  Is  capable  of  scanning  a  transducer  over  a  two- 
dimensional  profile  while  maintaining  an  angle  of  Incidence  perpendicular  to 
the  surface  and  maintaining  a  constant  water  path  distance.  The  system  can 
also  be  used  to  scan  a  tansducer  In  a  raster  pattern  for  generating  C-scan 


UUViflM  OiftV)  IMS  MX 


2.J  Block  diagram  of  multlaxls  controller  with  sample  profile. 


The  software  for  the  microprocessor-based  manipulator  control  sub¬ 
system  was  created  by  that  unit's  manufacturer.  It  Is  resident  In  an  EPROM, 
and  Is  therefore  relatively  easily  changed.  The  function  of  this  software  Is 
to  receive  command  characters  and  data  strings  from  the  minicomputer,  and 
translate  these  Into  pulses  to  the  stepping  motors  and/or  responses  to  the 
minicomputer.  There  are  13  command  characters  recognized  by  the  Z80  software, 
and  13  response  characters  which  It  generates,  some  of  which  are  followed  by 
signed  Integers.  Table  2.1  shows  a  list  of  these  26  characters  with  a  des¬ 
cription  of  their  meanings.  Data  Is  passed  to  the  Z80  in  a  sequence  of  ASCII 
characters  called  a  data  block.  An  example  of  this  Is  MN5F10X10Y1009B91$II< 
which  Is  interpreted  in  Fig.  2.4  and  can  be  read  "on  the  fifth  move  go  from 
the  present  position  to  X  *  0.010  Inches,  Y  *  1.000  Inches,  with  a  B  angle  of 
91°  at  a  speed  of  10%  of  the  maximum  (although  the  actual  maximum  speed  of  the 
axis  Is  limited  by  the  motor  to  2  Inches/second),  the  maximum  repetition  rate 
of  the  pulser  driving  the  motors  Is  3  KHz  corresponding  to  a  speed  of  3 
Inches/second.  It  Is  this  maximum  of  3  Inches/second  that  is  referred  to  as 
100%  of  maximum  speed.)  The  data  block  Is  decoded  and  stored  in  a  RAM  memory 
until  It  Is  time  for  execution  of  the  move.  Up  to  256  separate  data  blocks 
may  be  stored.  If  necessary,  before  movement  Is  Initiated  with  an  "I" 
character.  Execution  consists  of  calculating  the  distance  required  In  each 
coordinate,  converting  It  Into  stepping  motor  pulses  for  the  axes  Involved, 
and  sending  out  the  pulses  to  each  axis,  appropriately  Interleaved  to  provide 
maximum  accuracy  In  following  the  contour  of  the  move. 

The  Z80  also  may  receive  Input  from  the  "joystick,"  the  handheld 
manual  control  unit.  If  no  data  block  move  is  in  progress.  In  this  mode.  It 
Is  possible  to  reposition  any  of  the  6  axes  by  closing  a  toggle  switch  on  this 
unit  for  a  period  of  time.  Then  a  *G"  command  from  the  minicomputer  to  the 
190  will  return  the  correct  new  position  of  the  6  axis  system.  However,  It  Is 
necessary  after  one  or  more  such  manual  moves,  to  send  to  the  Z80  a  special 
"N0"  data  block  In  order  to  update  the  memory  of  the  Z80. 
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Fig.  2.4  Interpretation  of  data  block  command 


Table  2.1 

Command  Characters  for  Multi  axis  Controller 


Name 

To 

uP 

From 

u? 

Description 

HOLD 

D 

— 

Halt  exislting  processes 

CONTINUE 

E 

— 

Continue  halted  process 

GET  X 

X 

— 

Request  X  position 

GET  Y 

Y 

— 

Request  Y  position 

GET  Z 

Z 

— 

Request  Z  position 

GET  A 

A 

— 

Request  A  position 

GET  B 

B 

— 

Request  B  position 

GET  TT 

T 

— 

Request  TT  position 

X  RESPONSE 

— 

X 

Response  header  (sign  and  5  digits) 

Y  RESPONSE 

— 

Y 

Response  header  (sign  and  5  digits) 

Z  RESPONSE 

— 

Z 

Response  header  (sign  and  5  digits) 

A  RESPONSE 

— 

A 

Response  header  (sign  and  5  digits) 

B  RESPONSE 

— 

B 

Response  header  (sign  and  5  digits) 

TT  RESPONSE 

— 

T 

Response  header  (sign  and  5  digits) 

GET  ALL 

G 

— 

Get  all  axis  positions 

DATA 

N 

— 

Data  block  to  follow 

STATUS 

S 

— 

Request  status  of  microprocessor 

STATUS  RESPONSE 

— 

M 

Returned  status 

FAULT  RESPONSE 

— 

V 

Non-sol icited  hardware  error 

BLOCK  DATA 

REQUEST  RESPONSE 

___ 

W 

Non-sol i cited  request  for  data 

START 

I 

— 

Commence  movement  at  block  0000 

DATA  VALID 

RESPONSE 

Q 

Non-sol icited,  good  data  field 

DATA  NON- VALID 
RESPONSE 

— 

R 

Non-sol icited,  ignore  data  until  data 

valid  command  received 

ERROR  RESPONSE 

U 

U 

Unrecognized  message,  arguments 
may  or  may  not  exist 

INITIALIZE 

0 

— 

Initialize  to  reference  switches 

ACKNOWLEDGE 

RESPONSE 

— 

K 

Acknowledge  receipt  of  command 
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A  number  of  problems  have  been  encountered  with  the  operation  of  this 
system  which  require  more  operator  patience  and  skill  than  is  desired.  Many 
of  these  items  have  reputedly  been  corrected  on  newer  versions  of  the 
controller.  In  view  of  the  growing  acceptance  of  such  systems  by  the  NDE 
community,  it  seems  appropriate  to  list  some  of  the  items  that  are  sources  of 
error  or  frustration  so  that  others  who  need  to  select  or  specify  similar 
systems  can  avoid  the  same  pitfalls. 

A.  Problems  relating  to  limit  switches: 

1.  When  the  system  activates  a  limit  switch  under  manual 
control  it  does  not  stop. 

2.  When  the  system  activates  a  limit  switch  under  computer 
control,  it  stops  but  the  system  has  to  be  reset  before  it 
can  be  driven  out  of  the  limit  under  manual  control. 

3.  When  the  system  activates  a  limit  switch  it  dee?  sot  kee?  *• 
count  of  how  far  it  moved  beyond  the  limit. 

4.  The  limit  switches  on  the  manipulator'  do  not  work  well  and 
inhibit  the  operation  of  the  manipulator  for  some 
applications.  They  have  been  removed  from  our  system. 

B.  Problems  relating  to  the  reset  switch: 

1.  The  reset  switch  has  to  be  activated  twice  to  ensure  that 
the  system  is  completely  Initialized. 

2.  When  the  reset  switch  Is  activated,  the  axes  tend  to  move  a 
distance  corresponding  to  about  one  stepping  motor  step. 
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C.  Problems  related  to  the  manual  control: 


1.  If  an  N0  block  has  been  executed  and  an  initiate  (I) 
command  has  been  given,  the  system  will  keep  track  of  its 
position  accurately  as  long  as  it  is  driven  under  computer 
control.  However,  if  the  manual  control  is  used  to  move 
the  system,  the  Z-80  is  not  updated  at  the  end  of  the  move 
and  further  moves  under  computer  control  do  not  account  for 
the  manual  moves  that  have  occurred. 

2.  The  speed  of  an  axis  move  with  the  "fast"  manual  control  is 
too  slow  for  moves  greater  than  about  1  in.  It  would  be 
preferable  to  be  able  to  obtain  maximum  speeds  while  under 
manual  control. 

3.  There  are  n*/  position  encoders  to  monitor  the  motion  of  the 
axis  making  it  necessary  to  query  the  computer  to  determine 
position  when  moving  the  transducer  with  manual  control. 

4.  There  is  no  on-line  display  of  position. 

D.  Problems  related  to  turntable: 

1.  The  turntable  cannot  be  moved  to  a  specific  position  under 
computer  control . 

E.  Problems  related  to  passing  blocks  from  minicomputer  to  micro: 

1.  Carriage  returns  during  data  block  sending  are  not  always 
s  recognized.  Systematic  Interruptions  occur  approximately 

every  8  or  16  blocks,  although  the  events  seem  to  be 
aperiodic. 

* 
e 
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F.  Problems  related  to  contour  following  mode  of  operation: 

1.  Curvilinear  moves  are  Inaccurate.  The  stepping  motors  for 
the  Involved  axes  do  not  always  Interleave  the  pulses 
properly.  For  Instance,  the  X  and  B  axes  will  reach  their 
terminal  positions  while  Z  Is  continuing  to  move  to  reach 
Its  final  position.  This  means  the  final  position  is 
correct,  but  the  path  followed  to  reach  It  Is  not  correct. 

2.  Curvilinear  moves  cannot  cross  90°  quadrant  boundaries. 
Furthermore,  for  moves  that  have  to  cross  such  boundaries. 
It  Is  necessary  to  Insert  a  block  that  produces  a  linear 
move  between  the  curvilinear  moves  on  each  side  of  the 
boundary. 

3.  It  Is  necessary  for  the  coordinates  of  a  curvilinear  move 
to  be  provided  to  the  precision  of  the  machine  otherwise 
confusion  results  and  the  machine  moves  totally  incor¬ 
rectly.  The  machine  has  no  way  of  checking  data  and 
Indicating  that  a  problem  exists.  Thus  the  system  Is  not 
robust  to  operator  error. 

G.  Problems  related  to  error  messages: 

1.  There  are  a  variety  of  situations  In  which  the  system  falls 
to  function  or  will  function  incorrectly  but  there  is  no 
error  message  that  would  tell  the  operator  that  he  has  made 
an  error. 

H.  Problems  related  to  Initialization  of  system: 
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1.  Home  values  are  stored  In  an  EPROM  and  the  EPROM  has  to  be 
reprogrammed  to  change  them. 

I.  Problems  related  to  C-scans  and  "data  dropout": 

1.  There  are  no  trigger  pulses  available  that  can  be  used  to 
gauge  the  distance  that  the  transducer  has  moved  and  thus 
synchronize  a  pulser. 

The  personnel  who  have  used  the  system  believe  that  the  most  pressing 
need  Is  robustness  In  the  face  of  operator  error  and  better  diagnostics  when 
errors  are  encountered. 

2.3  Display  Processor 

The  Genlsco  display  processor  Is  specifically  Intended  to  accept 
command  and  data  strings  sent  from  the  minicomputer  and  to  translate  these 
Into  arrangements  of  pixels  in  a  480  *  512  RAM  memory.  The  8  bit  word 
associated  with  each  pixel  can  be  coded  as  one  of  16  levels  of  gray  or  as  one 
of  256  different  colors  and  displayed  on  the  Conrac  color  raster-scan  Image 
display  unit  (hereafter  referred  to  as  the  TV  monitor).  The  TV  monitor 
displays  a  525  line  raster  on  a  cathode  ray  tube  (CRT)  measuring  11.4  inches 
high  x  15.5  inches  wide. 

The  software  In  the  display  processor  was  created  by  Genisco  and  is 
for  all  practical  purposes  hardwired  and  unchangeable.  Twenty  basic  commands 
are  defined.  A  complete  description  of  them  Is  beyond  the  scope  of  this  text, 
but  they  Include  these  functions: 

1.  Display  of  alphanumeric  characters  with  any  combination  of  four 
different  orientations  and  four  different  character  sizes; 

2.  Replacement  of  the  bit  pattern  describing  any  pixel; 
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3.  Display  of  a  general  rectangle; 

4.  Display  of  a  vector  drawn  between  any  two  points; 

5.  Cursor  display  and  control; 

6.  Return  to  minicomputer  of  data,  such  as  the  position  of  the 
cursor; 

7.  Hardware  zoom  and  scroll; 

8.  Replacement  of  the  video  look-up  table  which  provides  the 
correspondence  between  bit  pattern  and  displayed  color. 

Communication  between  the  Genlsco  unit  and  the  Data  General 
minicomputer  takes  place  over  the  I/O  bus  of  the  minicomputer  at  Instantaneous 
data  rates  In  excess  of  4  megabits  per  second. 

2.4  Data  Acquisition  System 

The  data  acquisition  system  Is  composed  of  a  number  of  Interconnected 
units  shown  In  Fig.  2.5.  The  transducers  used  with  the  system  span  the  fre¬ 
quency  range  of  approximately  1  to  25  Wz.  Typically  these  are  wideband  units 
that  provide  short  transient  responses.  Both  focused  and  unfocused  units  are 
used.  The  pul ser/ receiver  Is  a  commercially  available  unit  from  Panametrlcs. 
None  of  Its  functions  can  be  computer  controlled  but  It  produces  a  short  pulse 
and  the  receiver  can  provide  up  to  40  dB  of  gain.  The  receiver  recovers  from 
saturation  quickly  and  an  Internal  attenuator  can  provide  up  to  68  dB  of 
attenuation.  It  has  functioned  well  for  the  applications  encountered  with  the 
Test  Bed.  The  output  of  the  receiver  goes  to  a  Tektronix  453A  oscilloscope 
used  to  observe  the  received  waveforms.  The  output  also  goes  to  a  waveform 
acquisition  unit  made  by  Blomatlon  that  digitizes  selected  portions  of  the 
received  signals.  It  contains  a  nominal  8  bit  analog  to  digital  converter 
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Fig.  2.5  Block  diagram  of  data  acquisition  system 


TO  S/200 
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with  a  maximum  sampling  rate  of  100  Wz.  Tests  of  this  unit  show  that  Its 
accuracy  decreases  with  Increasing  frequency  at  frequencies  above  500  kHz.  In 
order  to  obtain  sufficient  accuracy  In  the  sampled  data,  a  signal  averaging 
technique^  was  used.  This  technique  Involved  averaging  multiple  acquisitions 
of  the  waveform  while  varying  the  dc  offset  of  the  instrument  between 
acquisitions.  A  typical  procedure  which  provided  more  than  adequate  accuracy 
consisted  of  averaging  250  waveforms  In  groups  of  10  at  each  of  25  different 
dc  offsets. 

Some  special  software  has  been  designed  for  Interfacing  the  waveform 
acquisition  unit  to  the  minicomputer.  This  code  which  Is  partially  In 
assembly  language,  arms  the  unit,  waits  for  a  trigger,  collects  the  waveform 
and  then  processes  the  waveform  In  a  manner  which  is  specified  by  the 
operator.  Although  the  Inherent  speed  of  the  Biomatlon  unit  will  permit  500 
waveforms  of  1024  samples  each  to  be  acquired  each  second,  this  rate  has  never 
been  achieved  In  practice.  A  more  typical  waveform  acquisition  rate  Is  about 
100  per  second. 

The  ultrasonic  phased  array  system  Is  also  used  for  data 
acquisition.  This  system  will  be  fully  described  In  Section  4.0. 

2.5  Minicomputer  System 

The  minicomputer  for  the  Test  Bed  is  a  Data  General  Eclipse  S/200. 

The  unit  contains  a  162  K  byte  memory  and  can  be  partitioned  to  allow 
operators  at  two  separate  terminals  to  access  It.  It  has  been  equipped  with  a 
50  M  byte  disk  memory  and  9  track  magnetic  tape  unit. 

The  software  for  the  minicomputer  has  been  written  entirely  at  the 
Science  Center.  It  Is  largely  In  Fortran,  although  certain  portions  which 
access  the  Biomatlon  A/D  converter  are  In  assembly  language.  It  Includes  four 
different  categories  of  programs,  all  of  which  operate  Interactively  with  the 
user.  The  categories  are  listed  below: 
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1.  programs  for  generating  and  checking  disk  files  of  data  blocks 
which  may  be  interpreted  by  the  Z80; 

2.  multitasking  control  programs  for  managing  the  manipulator,  data 
acquisition,  and  storage  on  disk; 

3.  programs  for  displaying  the  data  In  color  or  black  and  white  via 
the  Genlsco  display  processor; 

4.  ISP,  a  signal  processing  and  display  program  of  enormous 
flexibility. 

All  but  ISP  have  been  created  for  the  Test  Bed  project.  A  brief  description 
of  each  of  the  categories  of  minicomputer  programs  with  the  exception  of  ISP 
follows. 

In  category  (1)  there  are  two  kinds  of  programs  for  generating  disk 
files  of  data  blocks.  One  kind  which  generates  a  file  corresponding  to  a 
raster  pattern  In  X  S  Y  contains  two  slightly  different  programs.  In  each 
program  the  pattern  of  moves  Illustrated  In  Fig.  2.6  Is  generated.  The 
pattern  begins  with  the  move  (AB,  BC)  and  then  repeats  the  pair  (CD,  DE)  until 
the  area  within  the  dashed  lines  Is  traversed.  One  program,  RASTER,  Is 
completely  flexible  and  allows  any  number  of  CO  moves,  and  any  length  of  DE 
and  CD.  Figure  2.7a  shows  the  user  terminal  screen  display  which  occurs 
during  operation  of  RASTER  with  both  prompts  from  the  minicomputer  and  user 
Inputs.  User  Inputs  are  outlined  with  boxes.  The  other  program,  RASTMAG,  Is 
limited  to  creating  one  of  seven  raster  patterns.  Each  of  these  patterns  has 
125  (CD,  DE)  moves.  The  pattern  is  exactly  square  and  mqy  have  one  of  seven 
different  edge  lengths  0.25  Inches,  0.50  Inches,  1.0  Inches,  2.5  Inches,  5 
Inches,  10  Inches,  and  25  Inches.  This  represents  a  variation  of  10*  In  area 
In  order  to  service  parts  or  samples  ranging  In  size  from  a  transistor 
assembly  to  a  turbine  disk.  Figure  2.7b  shows  the  screen  display  for 
operation  of  RASTMAG. 
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SCHEMATIC  CROSS  SECTION 
OF  ULTRASONIC  BEAM  ^ 


Fig.  2.6  A  schematic  of  the  pattern  of  moves  followed  In  an  X-Y  raster  scan. 
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I  RASTER 


Y: 

f5~~5 

ENTER  FINAL  X.  Y 

NUMBER  OF  RASTERS  (I  E.  #  OF  X-SCANS) 
ENTER  NAME  OF  NEW  DISK  FILE:  (7  CHAR  LOB 
1  EDGE  010  DB 
R 


(a) 


1 R  ASTMAG  I 

FOR  SQUARE  OF  25.00 
MAGNIFICATION  *  00.40 


10.00  5.00  2.50  1.00 
1.00  2.00  4.00  10.00 


00.50  00.25  INCHES 

20.00  40.00 


PLEASE  ENTER  MAGNIFICATION  DESIRED:  fTo] 
PLEASE  ENTER  RASTER  CENTER  IN  X,  Y:  |  0,0 1 
ENTER  NAME  OF  NEW  DISK  FILE:  (7  CHARLDB 

Imagnfio.dbI 

R 


(b) 


Fig.  2.7  User  terminal  display  after  operation  of  RASTER  and  RAST  MAG.  User 
Inputs  are  delineated  by  boxes. 
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The  second  kind  of  data  block  generation  program  Is  represented  by 
PROFILE,  which  Is  primarily  for  creation  of  data  blocks  corresponding  to  the 
profile  of  a  complex  part  shape  with  rotational  symmetry  such  as  a  turbine 
disk.  Figure  2.8  shows  the  screen  display  for  operation  of  PROFILE  to  create 
a  disk  file  corresponding  to  the  contour  shown  In  Fig.  2.9a.  The  resulting 
disk  file  Is  shown  In  Fig.  2.9b.  At  a  few  points  In  the  contour.  It  Is 
necessary  to  change  the  entry  angle  of  the  acoustic  beam  without  changing  the 
entry  point  on  the  part  surface.  During  this  change,  it  may  be  Important  not 
to  record  data.  Ordinarily,  as  execution  of  each  data  block  is  begun,  an  "R" 
Is  sent  from  the  micro,  rapidly  followed  by  a  “Q”  after  the  final  speed  has 
been  reached.  However,  If  “M601'  has  been  specified  in  the  data  block,  the  "Q" 
Is  never  sent  and  data  Is  not  taken  during  that  move.  Notice  that  PROFILE  has 
automatically  Inserted  non-data-taking  moves  at  the  appropriate  places.  These 
can  be  easily  modified  to  omit  the  "M60"  If  data-taking  is  desired  during 
those  moves.  The  arc  centers  MI"  and  "K"  are  specified  If  a  circular  contour 
must  be  followed  during  the  move. 

In  order  to  provide  a  check  on  the  data  blocks  generated  with 
PROFILE,  a  program  named  PROPLOT  has  been  coded  that,  for  any  given  contour 
scan,  makes  a  plot  of  the  part  profile  and  two  paths  associated  with  the 
transducer  manipulator.  It  needs  as  Input  from  the  operator  only  the  name  of 
the  data  file  containing  the  data  blocks  and  the  length  of  the  transducer. 

One  path  is  the  path  the  center  of  the  transducer  face  makes  during  the  scan; 
the  other  Is  the  path  made  by  the  manipulator  pivot  point.  Figure  2.10  shows 
one  example  of  the  user-computer  Interaction  with  user  Inputs  shown  In 
boxes.  Data  non-valld  moves  are  flagged  and  the  values  of  two  Internal 
variables  are  displayed  for  diagnostic  purposes. 

The  lower  half  of  Figure  2.11  shows  the  result  of  executing  PROPLOT; 
the  part  profile  appears  with  the  two  paths  of  Interest.  Use  of  color  allows 
an  operator  to  make  an  easy  check  for  Interference:  the  part  profile  Is 
white,  the  pivot  point  path  Is  red,  and  the  transducer  path  Is  green.  For 
purposes  of  Illustration,  the  data  blocks  used  to  generate  the  graphic  display 
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Fig.  2.8  User  terminal  display  after  operating  PROFILE.  User  Inputs  are 
delineated  by  boxes. 
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Fig.  2.10  User  terminal  display  during  operation  of  PROPLOT.  User  inputs  are 
delineated  by  boxes. 
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.11  Output  display  of  PROPLOT  showing  part  profile,  transducer  path 
and  pivot  point  path. 
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have  been  added  to  the  upper  half  of  this  monitor  photograph;  this  information 
does  not  normally  appear.  The  labels  under  the  drawing,  however,  are  standard 
for  this  display. 

In  category  2)  we  have  the  multitasking  control  program,  named  TBCON , 
which  manages  the  manipulator,  data  acquisition,  and  storage  on  disk.  It 
maintains  three  tasks.  Task  1  communicates  with  the  operator  via  the  user 
terminal;  its  standard  prompt  is  "WAITING  FOR  CONTROL  WORD."  It  allows  one  to 
send  characters  directly  to  the  Z80  to  construct  and  send  data  blocks,  to 
start  and  stop  data-taking,  or  to  perform  other  useful  tasks  such  as 
calculating  the  refraction  of  the  acoustic  beam  through  a  flat  part  surface. 
Task  2,  which  has  priority  over  all  the  other  tasks,  reads  and  Interprets 
response  characters  from  the  Z80.  Task  3  operates  the  Biomatlon  A/D  converter 
as  described  in  Section  2.4,  and  stores,  at  intervals,  pertinent  data  on  disk 
for  later  access  by  display  programs.  Communication  between  tasks  is  via  an 
area  of  memory  common  to  all  three.  Table  2.2  shows  a  list  of  the  commands 
recognized  by  Task  1  along  with  their  functions.  The  commands  "RI"  and  "RF" 
are  explained  later.  The  command  "TW"  calls  a  subroutine  which  varies  the 
angles  A  and  B  separately  over  a  small  range  and  picks  the  values  which  yield 
the  maximum  signal.  That  is,  it  "tweaks  up"  A  and  B.  TBCON  provides  the  user 
with  three  modes  of  acquisition  which  are  detailed  separately  below,  and 
correspond  to  C-scans,  B-scans,  and  A-scans.  B-  and  C-scans  are  performed  in 
a  raster  fashion,  either  rectangular  or  circular,  i.e.,  in  increments  of  X  and 
Y  of  a  cartesian  coordinate  system,  or  increments  of  radius  and  angle  of  a 
polar  coordinate  system.  If  the  user  specifies  a  circular  scan,  additional 
position  Information  is  stored  on  disk  In  Task  3. 

MODE  1;  C-scan.  During  a  scan,  only  the  minimum  and  maximum  peak 
amplitudes  of  each  waveform  are  actually  stored  on  disk;  this  provides  a  rapid 
scan  of  a  horizontal  plane  of  an  entire  sample.  Figure  2.12  shows  a  typical 
C-scan  start-up  for  TBCON,  with  user  inputs  outlined  by  boxes.  In  this  case, 
the  user  chose  to  take  data  with  one  of  the  seven  pre-determl ned  rectangular 
raster  patterns.  After  sending  over  the  data  blocks  contained  in  MGNF010.DB, 
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2.12  User  terminal  display  during  operation  of  TBC0N.  User  Inputs  are 
delineated  by  boxes. 


32 


Table  2.2 

Commands  Recognized  by  Task  1  of  TBCON 


User 

Task  1 

Command 

Response 

Any  single  letter:  A,  8,  etc. 

Sends  that  character  directly  to  the  Z80 

08 

Allows  entry  of  a  data  block  for  transmission 
to  the  Z80. 

OF 

Reads  a  file  of  data  blocks  from  disk  and 
sends  them  to  the  Z80. 

8E 

Starts  Task  3,  plus  sends  "I"  to  Z80. 

HA 

Suspends  Task  3,  sends  "D"  to  Z80. 

CO 

Restarts  Task  3,  plus  sends  "E"  to  Z80. 

FI 

Suspends  Task  3,  sends  “D,"  closes  the 
output  file. 

CH 

Allows  a  change  In  threshold  while  data 
taking  and  movement  are  momentarily  halted 

RI,  RF 

See  text 

TW 

See  text 

N0 

Gets  current  position  of  axes;  returns  a  data 
block  beginning  with  NjJ  with  current  position 
specified.  This  updates  the  Z80's  memory 
after  a  joystick  move. 

X+l  m. 

This  type  of  command  constructs  a  data  block 

or  Y-1232, 

consisting  of  N1X1(J|J0  or  N1  Y-1232  or 

or  Z-5jJ,  etc. 

NIZ-50  and  sends  it  to  the  Z80. 

XX 

Closes  output  file  and  suspends  all  tasks, 
effectively  ending  the  program. 
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the  standard  prompt  recurs,  and  the  user  types  In  HBE.H  From  there  on,  the 
scan  Is  automatic  and  the  presence  of  an  operator  Is  not  necessary  until  the 
data  Is  ready  for  display.  For  this  program,  start-up  takes  less  than  five 
minutes  and  a  scan  of  a  5  Inch  by  5  Inch  area  with  125  raster  lines  takes 
about  20  minutes. 

While  sending  the  data  blocks  to  the  Z80,  Task  1  Is  also  storing  them 
In  the  output  file  (on  disk).  When  movement  begins.  Task  2  goes  Into  a  con¬ 
tinuous  loop  which  Is  broken  only  by  occasional  disk  storage  activities.  The 
rate  of  this  loop  Is  set  by  the  Panametrlcs  Pulser/Recelver  which  provides  the 
trigger  for  the  digitizer  unit.  All  the  data  pairs  for  a  given  move  (See 
Fig.  2.6)  are  stored  In  a  sequence  on  the  disk.  In  this  manner,  the  resulting 
data  output  file  has  as  many  data  pair  sequences  as  It  has  J3)  data  blocks. 

This  rectangular  scan  Is  the  fastest  possible  way  of  obtaining  a  full 
Image  or  C-scan  because  It  does  not  query  the  Z80  after  movement  Is  Initiated. 
The  current  data-taklng  rate  In  this  mode  can  be  as  high  as  150  Hz,  which 
corresponds  to  a  part  sampling  rate  of  417  points  per  inch  at  the  normal  scan 
speed  of  0.36  Inches  per  second.  A  circular  scan  queries  the  Z80  and  stores 
position  Information  for  each  waveform.  Although  its  maximum  data-taklng  rate 
Is  only  80  Hz,  this  mode  lends  Itself  to  more  complicated  scans  with  other 
than  pure  raster  moves.  For  example,  this  Is  used  during  scans  Involving 
rotation  of  the  turntable. 

In  C-scan  operation,  a  further  option  allows  the  user  to  specify  an 
amplitude  threshold  whose  value  can  vary  as  a  function  of  depth  below  the 
surface  of  the  part.  The  amplitude  of  a  waveform  is  only  stored  If  It  exceeds 
the  threshold  value.  A  varying  threshold  permits  adjustment  of  the  detection 
sensitivity  to  compensate  for  material  attenuation,  focusing  effects,  front 
surface  "ring-down",  etc. 

MODE  2:  B-scan.  During  a  scan.  Information  about  the  entire  waveform 
Is  kept,  as  opposed  to  only  the  maximum  and  minimum  values  that  are  kept  In  a 
C-scan.  Scans  may  be  performed  In  either  a  rectangular  or  circular  fashion. 
Data  storage  occurs  more  frequently  than  during  the  C-scan  mode,  and  Includes 
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position  and  amplitude  Information,  Identified  by  X,  Y,  and  Z  coordinates 
where  Z  represents  the  position  In  time  of  a  point  within  the  waveform. 


The  speed  of  a  B-scan  decreases  proportional  to  the  amount  of  data 
recorded.  As  In  the  previous  mode,  the  user  may  specify  a  time-varying 
threshold  which  must  be  exceeded  before  the  data  will  be  stored.  The  thres¬ 
hold  may  be  described  as  constant  or  exponentially  decaying.  This  option 
provides  the  capability  of  storing  only  portions  of  waveforms  and  thus 
Increases  the  data  acquisition  speed. 

MODE  3:  A-scan.  The  third  mode  of  operation  Involves  using  TBCON  as 
a  sub-task  of  ISP.  In  this  mode,  TBCON  provides  position  control  and  calcu¬ 
lates  the  Fresnel  reflection  loss  at  the  water/metal  Interface.  ISP  performs 
data  acquisition  and  analysis  of  the  data.  This  mode  Is  employed  when  It  is 
necessary  to  record  entire  waveforms  for  the  detailed  examination  of  a  flaw. 
In  the  course  of  such  an  examination.  It  Is  desirable  to  acquire  data  from 
arbitrary  angles  relative  to  the  surface  of  the  part.  Figure  2.13  shows  the 
refraction  of  the  acoustic  beam  at  the  surface  of  the  sample.  The  angle  of 
Incidence  Is  designated  0e  and  the  angle  of  refraction  Is  designated  e^.  The 
TBCON  commands  “RIM  and  “RFM  are  used  to  call  a  subroutine  which: 

1.  calculates  the  position  of  the  pivot  point  in  X,  Y,  and  Z, 

2.  calculates  the  glmbal  angles  A  and  B, 

3.  sends  the  appropriate  data  block  to  the  Z80, 

4.  calculates  the  Fresnel  reflection  loss  at  the  water/metal 
Interface. 

Then  data  Is  acquired  at  that  angle  with  ISP  and  can  be  stored  on  a  disk  file 
for  further  analysis. 
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The  third  category  of  programs  comprise  the  routines  that  have  been 
written  so  the  S/200  can  control  the  Genlsco  display  processor.  Twenty  of 
them  are  basic  to  the  display  functions  and  are  gathered  into  a  library, 

PGP. LB,  which  Is  listed  In  Table  2.3.  Two  programs  have  been  written  to 
display  the  results  of  the  first  two  data-taklng  modes  of  TBCON.  C-scan  data 
obtained  In  MODE  1  Is  displayed  by  TBOISC,  while  MODE  2  B-scan  data  is  dis¬ 
played  by  TBDISZ.  The  two  programs  are  very  similar  and  an  explanation  of  the 
operation  of  TBDISC  should  suffice  for  both  programs. 

The  interaction  between  the  user  and  computer  to  obtain  the  display 
of  a  C-scan  image  Is  shown  In  Fig.  2.14;  user  Inputs  are  outlined  with 
boxes.  The  program  has  a  number  of  options  Including: 

1.  Clearing  or  not  clearing  the  monitor; 

2.  Rectangular  raster  or  spiral  format  data; 

3.  Choice  of  seven  magnifications; 

4.  Arbitrary  reference  point  to  establish  an  absolute  coordinate 
system; 

5.  Use  of  most  positive  or  most  negative  signal  acquired; 

6.  Relative  sizing  options  such  that  more  than  one  Image  may  be 
displayed  at  once  (up  to  nine); 

7.  Alteration  of  the  contrast  and  brightness  from  the  calculated 
values; 

8.  Offset  of  alternate  rasters  In  the  X  direction; 
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Table  2.3 


Library  of  Display  Processor  Routines 


PGP. LB 

Function  Name 

Description  of  Function 

VLTGS 

Generates  a  gray  scale  and  updates  the  VLT 
with  It. 

VLTNS 

Generates  a  blue,  green,  red,  white  scale  and 
updates  the  VLT  with  It. 

COLOR 

Allows  selection  of  color  to  apply  to  other 
functions,  with  optional  pre-empt. 

WRSTR 

Sends  an  alphanumeric  string  to  the  display 
processor  for  display  on  the  monitor 

XYZOO 

Zooms  display  In  a  region  centered  at  an 
arbitrary  X,  Y,  to  one  of  four  factors:  1,  2, 

4,  or  8. 

CUREN 

Enables  a  cursor  of  any  color  combined  with 
four  different  shapes. 

WRZRA 

Writes  a  raster  of  data. 

WRVEC 

Writes  a  vector  given  beginning  and  end  points 

CURD  I 

Disables  a  previously  enabled  cursor. 

CURTR 

Allows  user  to  track  with  cursor  to  a  desired 
position  and  returns  X,  Y  coordinates  to 
minicomputer. 

SETPL 

Turn  on  selected  memory  planes. 

SCROL 

Scrolls  display  In  X  and  Y. 

WRREC 

Writes  a  cartesian  rectangle  at  an  arbitrary 

X,  Y  with  optional  color  fill. 

SETUP 

Clears  all  parameters  to  their  default 
values. 

WRVLT 

Writes  a  user  defined  video  lookup  table  to 
the  display  processor. 

LENGT 

Function  returning  the  length  of  an  alpha 
numeric  array. 

ZOOM 

Zooms  top  left  area  of  screen  to  desired 
factor. 

CLRPL 

Clears  selected  memory  planes. 

ARGCN 

Counts  the  number  of  arguments  In  a 
subroutine  call. 

PGPDR 

Interface  between  minicomputer  CPU  and 
display  processor. 
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Fig.  2.14  User  terminal  display  after  operation  of  T8DISC. 
delineated  by  boxes. 


User  inputs  are 
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9.  Alteration  of  VLT; 


10.  Use  of  cursor  to  zoom  any  portion  and  measure  any  point  in 
absolute  coordinates. 

Two  examples  of  Images  resulting  from  the  display  of  MODE  1  data  are  shown  In 
Fig.  2.15.  It  displays  the  amplitude  of  reflected  sound  from  two  different 
voids  In  titanium,  translated  into  a  gray  scale  video  display.  The  transducer 
was  a  focused  beam  piezoelectric  unit  operating  at  5  MHz.  The  image  on  the 
left  Is  from  a  0.50  inch  by  0.50  inch  region  centered  on  a  spherical  void  of 
800  micron  diameter;  was  0°.  The  Image  on  the  right  Is  from  a  0.25  Inch  by 
0.25  Inch  region  centered  on  an  oblate  spheroid,  250  microns  by  800  microns 
(see  Fig.  2.13).  The  Internal  angle,  0^,  was  60°.  The  file  name  and  the  key 
to  the  color  vs  Intensity  code  or  video  look  up  table  (VLT)  are  displayed 
automatically  at  the  bottom  of  the  Image. 

In  addition,  a  large  variety  of  display  software  has  been  created  for 
testing  various  subroutines,  for  displaying  data  in  3D  perspective,  for  chang¬ 
ing  the  VLT,  for  using  the  cursor,  for  labeling,  and  other  auxiliary  func¬ 
tions. 

The  software  that  Is  used  for  processing  the  signals  acquired  during 
the  Investigation  of  a  particular  flaw  Is  known  as  the  interpretive  j>1gnal 
Processing  (ISP)  program.  This  software  was  created  at  the  Science  Center 
separately  from  the  Test  Bed  program  and  Is  unique  in  the  field  of  NDE.  ISP 
Is  a  large  multitasking  program  that  runs  on  the  S/200  under  the  command  of  an 
external  terminal.  It  can  be  used  to  string  together  a  number  of  signal 
processing  operations  In  a  group  that  Is  designated  a  "procedure".  There  are 
over  40  distinct  operations  that  can  be  used  In  a  "procedure".  Those 
operations  pertinent  to  the  Test  Bed  Include: 

1.  Waveform  acquisition  (in  conjunction  with  an  A/D  converter). 
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2.  Waveform  transformation  Including:  addition,  subtraction, 
Fourier  transforming,  time  shifting,  multiplication,  division, 
etc. 

3.  Feature  extraction  including:  maxima,  minima  centroids,  energy, 
statistical  data,  etc. 

4.  Display  of  data  including:  automatic  scaling,  linear  plots,  log 
plots,  etc. 

5.  Position  control  by  calling  TBCON  as  a  sub-task. 

The  "procedures"  employing  these  operations  can  read  single  or  multiple  record 
files  from  the  disk  memory  of  the  S/200  and  modify  them  with  sequential  pro¬ 
cessing  operations.  There  can  also  be  internal  looping  within  a  "procedure". 
The  algorithms  used  for  the  inversion  technique  to  be  described  in  Section  3.0 
are  designed  using  the  ISP  procedures  to  take  advantage  of  the  flexibility  and 
ease  of  use.  Once  the  utility  of  an  algorithm  has  been  established,  the 
required  element  of  ISP  can  be  abstracted  and  Incorporated  into  a  stand-alone 
software  routine. 

In  the  case  of  the  Born  Inversion  algorithm,  a  combination  of  three 
ISP  procedures  entailing  over  100  separate  operations  was  constructed  and 
employed  for  Inversion  measurements  under  this  contract.  Afterwards,  In  a 
separately  funded  effort, ^  the  required  ISP  functions  were  extracted  and 
Incorporated  into  a  new  program  called  BIP  which  Is  an  Interactive  stand-alone 
program. 
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3.0  INSPECTION  PROCEDURES 


3.1  System  Protocol 

The  protocol  for  the  system  is  designed  to  permit  the  acquisition  of 
data  from  a  part  in  a  logical  fashion  that  will  permit  the  outputs  of  the 
available  analysis  techniques  to  be  combined  to  produce  the  best  estimate  of 
the  flaw  characteristics.  In  other  words  the  goal  is  to  maximize  the 
synergism  that  can  be  obtained  from  the  available  analysis  techniques. 


There  is  a  hierarchy  associated  with  these  analysis  techniques.  The 
fundamental  inversion  techniques  each  use  a  specific  portion  of  the  flaw's 
frequency  spectrum  to  extract  estimates  of  various  flaw  characteristics  such 
as  size,  shape,  orientation  and  composition.  These  characteristics  are  then 
combined  using  a  universal  Inversion  algorithm.  This  algorithm  also  uses  £ 
priori  Information  such  as  the  probabilities  of  having  particular  flaw  types, 
noise  levels,  etc.  The  output  Is  a  probabilistic  estimate  of  the  flaw 
characteristics  that  is  based  on  all  the  data  available  from  that  flaw.  This 
estimate  Is  then  used  as  an  Input  to  a  fracture  mechanics  based  accept/reject 
criteria  to  decide  whether  the  part  is  acceptable. 


Although  the  general  structure  for  the  protocol  Is  generic  In  nature 
and  can  be  applied  to  many  different  systems  and  flaw  types,  the  Test  Bed  has 
been  specifically  applied  to  the  ultrasonic  examination  of  bulk  flaws  In 
materials.  Consequently  the  protocol  has  been  adapted  to  these  specific 
needs.  The  technique  consists  of  a  search  mode  for  the  detection  of  tentative 
flaw  Indications  followed  by  a  detailed  examination  of  these  tentative  flaws 
to  extract  as  much  quantitative  Information  about  them  as  possible.  A  diagram 
outlining  the  protocol  Is  shown  in  Fig.  3.1. 


The  description  of  the  protocol  Is  logically  coherent  only  when  taken 
as  a  whole.  Although  as  many  features  of  the  protocol  have  been  implemented 
as  time  r  /mltted,  all  features  have  not  been  Implemented.  The  descriptions 
contained  In  the  following  sections  point  out  the  current  state  of 
Implementation  of  a  particular  technique  or  procedures. 


Fig.  3.1  Inspection  protocol. 


3.1.1  Search  Mode 

The  search  mode  is  implemented  by  scanning  a  transducer  over  the  part 
in  some  regular  pattern.  For  a  flat  part  this  can  be  a  back  and  forth  raster 
pattern.  In  the  case  of  a  shaped  disk,  the  transducer  follows  the  contour  of 
the  disk  as  it  moves  radially.  The  disk  is  mounted  on  the  turntable  and 
rotates  beneath  it  to  produce  a  concentric  circle  pattern.  The  transducer 
used  for  this  inspection  is  either  unfocused  or  weakly  focused  to  provide 
complete  coverage  of  a  part  with  a  coarse  scanning  pattern.  A  tentative  flaw 
indication  is  identified  as  any  signal  that  exceeds  a  preset  threshold.  When 
a  flaw  indication  is  detected  its  location  is  stored  in  memory  so  that  a 
detailed  inspection  of  each  tentative  flaw  can  be  made  at  the  conclusion  of 
the  search  mode  scan. 

The  transducer  selected  for  the  search  mode  depends  on  the  nature  of 
the  part  that  is  to  be  Inspected,  the  types  of  flaws  that  are  anticicpated, 
and  the  required  speed  of  the  Inspection.  The  transducer  frequency  depends  on 
the  thickness  of  the  part  and  the  attenuation  of  the  material.  The  focusing 
properties  of  the  transducer  depend  on  the  depth  range  over  which  flaw  indi¬ 
cations  must  be  detected  as  well  as  the  expected  size  of  the  flaws  and  the 
time  allotted  for  the  scan. 

3.1.2  Detailed  Examination  Mode 

At  the  conclusion  of  the  search  mode,  each  site  of  a  tentative  flaw 
Is  returned  to  In  turn  and  a  detailed  examination  of  that  flaw  signal  is  per¬ 
formed.  The  detailed  examination  utilizes  several  techniques  that  have 
validity  In  different  spectral  regions  or  for  different  ranges  of  ka  where  a 
is  a  representative  dimension  of  the  flaw  and  k  Is  the  wavenumber  (2it/x)  in 
the  material.  The  short  wavelength  region  or  high  ka  range  Is  covered  by 
imaging  techniques  which  produce  B-  and  C-scan  displays.  The  medium 
wavelengths  are  covered  by  model  based  reconstruction  techniques.  The  Test 
Bed  has  Implemented  the  particular  technique  known  as  the  one  dimensional 
Inverse  Born  algorithm  for  obtaining  flaw  Information  In  this  domain.  The 
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long  wavelength  technique  provides  information  about  the  location  of  the 
center  of  the  flaw  as  well  as  the  coefficient  of  the  term  in  the  expansion 
of  the  scattering  spectrum  of  the  flaw  which  is  commonly  designated  Ap.  The 
location  of  the  flaw  center  is  essential  as  an  input  for  the  inverse  Born 
technique  and  the  coefficient  A2  can  provide  information  about  the  composition 
of  the  flaw.  The  determination  of  A2  has  not  been  implemented. 

3. 1.2.1  Imaging 

The  detailed  inspection  begins  with  the  imaging  of  a  localized  region 
surrounding  the  flaw  site.  Either  a  single  element  focused  transducer  or  a 
multielement  phased  array  transducer  can  be  used.  If  a  single  element  focused 
transducer  selected,  it  should  operate  at  the  highest  frequency  that  is  con¬ 
sistent  with  the  required  penetration  depth  and  the  transducer  availability. 

In  practice  the  focal  length  of  the  transducer  is  limited  by  the  finite  number 
of  available  transducers.  However  given  the  availability  of  suitable  trans¬ 
ducers  the  focal  length  should  be  determined  by  the  expected  size  and  depth  of 
the  flaws  as  well  as  by  the  depth  of  focus  required  for  the  Inspection.  At 
the  present  time,  no  provision  Is  made  for  Imaging  through  curved  surfaces. 
Cylindrlcally  curved  and  spherically  curved  surfaces  could  be  compensated  for 
by  fitting  the  transducer  with  auxiliary  lenses  but  no  efforts  have  been  made 
to  do  this.  This  localized  scan  can  provide  either  a  C-scan  with  an  operator 
selectable  depth  slice  displayed  on  the  monitor  or  a  B-scan  corresponding  to  a 
single  scan  line  passing  through  the  flaw  Indication.  If  a  multielement 
phased  array  transducer  Is  selected.  It  will  be  able  to  rapidly  scan  an  ultra¬ 
sonic  beam  to  acquire  B-scan  data.  It  will  be  able  to  change  the  focus  of  the 
beam  In  one  dimension.  This  will  permit  the  accommodation  of  various  curva¬ 
tures  In  a  cylindrical  surface. 

The  results  of  the  Imaging  Inspection  can  be  classified  according  to 
the  completeness  with  which  the  flaw  Is  resolved.  If  the  flaw  Is  sufficiently 
large  that  It  Is  completely  resolved  by  the  focal  spot  of  the  scanning 
transducer  then  the  Image  can  provide  a  complete  description  of  the  flaw's 
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size,  shape,  and  orientation.  The  protocol  would  then  move  directly  to  the 
accept/reject  criteria  to  determine  the  suitability  of  the  part  for  additional 
usage.  It  may  only  be  possible  to  resolve  selected  dimensions  of  the  flaw.  In 
which  case  the  rough  shape  of  the  flaw  and  its  orientation  can  still  be 
determined  from  the  image  and  used  in  the  next  step  of  the  protocol.  If  none 
of  the  flaw  dimensions  are  resolvable,  then  the  image  can  still  provide  a 
mapping  function  to  determine  how  many  individual  flaws  are  represented  by  the 
one  flaw  Indication.  Since  the  Inversion  techniques  Involving  medium  and  long 
wavelengths  only  work  with  Isolated  flaws,  it  is  essential  to  have  this 
Information  before  applying  them. 

3. 1.2. 2  Selection  of  Observation  Directions  and  Estimation  of  Flaw  Size 

To  obtain  an  Image  It  Is  necessary  to  acquire  many  waveforms  but 
these  are  acquired  with  very  regular  spacing  of  the  transducer  positions  and 
the  data  Is  obtained  quickly.  Some  of  the  model  based  reconstruction 
techniques  also  require  a  large  number  of  waveforms.  However  these  waveforms 
usually  have  to  sample  a  sphere  or  hemisphere  surrounding  the  flaw  in  a  very 
specific  way  to  insure  a  good  Inversion.  When  the  spatial  window  available  to 
Inspect  the  flaw  Is  restricted  by  the  geometry  of  the  part,  these  techniques 
Incur  large  errors  or  fall  to  work  altogether.  In  contrast,  the  Born 
Inversion  technique  requires  only  a  single  waveform  to  obtain  a  single  size 
estimate  corresponding  to  the  projection  of  the  flaw  along  the  propagation 
direction.  Although  a  reasonably  detailed  image  of  a  flaw  can  be  obtained  by 
using  the  Born  inversion  technique  with  waveforms  taken  at  many  different 
angles,  this  Is  unlikely  to  be  done  in  a  practical  situation.  The  size  of  the 
flaw  along  two  or  three  critical  directions  is  probably  all  that  will  be 
required.  These  directions  are  selected  on  the  basis  of  the  C-  and  B-scan 
Images  as  well  as  the  available  spatial  windows  allowing  access  to  the  flaw. 

Once  the  observation  directions  are  selected,  waveforms  are  acquired 
from  the  flaw  with  transducers  whose  bandwldths  encompass  both  the  long  and 
medium  wavelength  regime.  These  waveforms  are  then  processed  to  obtain  the 
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flaw  size  for  each  observation  direction.  The  exact  procedure  for  doing  this 
currently  requires  some  iteration  of  the  measurements  and  the  judgment  of  an 
operator.  The  prospects  for  making  this  procedure  completely  automatic  are 
favorable  but  will  require  additional  development. 

3. 1.2. 3  Combination  of  All  Data 

After  these  flaw  parameters  are  determined,  all  of  the  available 
information  about  the  flaw  including  any  a  priori  information  is  combined 
using  an  algorithm  that  produces  an  estimate  of  the  flaw  size,  shape, 
orientation,  and  composition. 

The  Test  Bed  program  has  been  chiefly  concerned  with  the  imaging 
techniques  and  with  the  one-dimensional  Born  inversion  technique  for  esti¬ 
mating  flaw  sizes.  The  long  wavelength  portion  of  the  spectrum  provides 
information  about  the  location  of  the  center  of  the  flaw  and  the  long  wave¬ 
length  scattering  coefficient  A2.  The  information  about  the  center  of  the 
flaw  has  been  used  extensively  in  conjunction  with  the  Born  inversion 
algorithm  but  little  work  has  been  done  to  measure  A2.  The  algorithms  for 
combining  the  data  to  obtain  a  probabll istlc  answer  for  the  flaw  size,  shape 
and  composition  has  been  incorporated  into  the  Test  Bed  software  resident  In 
the  mini -computer.  Because  complete  Input  data  was  not  available  for  any  of 
the  flaws  that  were  examined  l.e.,  A2,  signal  to  noise  data  for  the  signals, 
and  probabilistic  data  concerning  different  flaw  compositions,  the  algorithm 
was  not  utilized.  An  Informal  approach  has  been  used  to  combine  data  from  the 
Imaging  technique,  the  Born  inversion  and  data  obtained  about  the  Impedance  of 
the  flaw  relative  to  the  host  material.  From  this  size,  shape  orientation  and 
some  Information  about  the  flaw's  composition  was  deduced.  Additional  Infor¬ 
mation  on  this  experiment  Is  presented  In  3. 2. 4. 4. 


3.2  Quantitative  Flaw  Sizing  Techniques 

The  Inversion  technique  required  for  the  estimation  of  flaw  size 
encompasses  the  medium  wavelength  regime  as  defined  above.  The  flaw  sizing 


technique  Is  the  one-dimensional  Born  Inversion  algorithm  which  requires  an 
Input  from  the  long  wavelength  measurement  regarding  the  location  of  the 
center  of  the  flaw.  This  technique  Is  strictly  valid  only  for  ellipsoidal, 
weakly  scattering  flaws;  however,  In  practice  It  appears  to  work  reasonably 
well  for  nonelllpsoldal  strongly  scattering  flaws  that  have  a  convex  shape. 
This  particular  algorithm  has  been  found  to  be  robust  In  the  presence  of  low 
signal  to  noise  ratio  and  limited  transducer  bandwidth.  Before  using  the 
algorithm  It  Is  necessary  to  ascertain  if  the  bandwidth  of  the  transducer 
encompasses  the  appropriate  spectral  region  of  the  flaw  being  examined.  This 
Is  a  transducer  selection  problem  that  Is  solved  by  examining  the 
experimentally  measured  spectrum  of  the  flaw.  The  spectrum  Is  obtained  by 
sequentially  acquiring  signals  from  the  flaw  with  transducers  having  suc¬ 
cessively  lower  center  frequc  ides.  The  transducer  characteristics  are 
deconvolved  from  the  measured  flaw  signals  and  the  magnitude  spectrum  is 
examined  for  peaks.  Below  some  frequency  there  will  be  no  additional  peaks  in 
the  spectrum  but  only  an  f*  dependence  of  the  spectral  magnitude.  The 
transducer  whose  bandwidth  encompasses  the  lowest  peak  of  the  magnitude 
spectrum  Is  the  appropriate  one  to  use  for  acquiring  Born  inversion  data. 


3.2.1  Signal  Preprocessing  Techniques 

There  Is  some  pre-processing  required  for  the  signals  acquired  from 
the  flaw  and  this  Is  diagramed  In  the  flow  chart  shown  In  Fig.  3.2.  The  main 
objective  of  this  pre-processing  Is  to  Isolate  the  flaw  signal,  and  force  the 
signal  spectrum  to  zero  in  those  regions  where  the  transducer  has  insufficient 
energy.  To  accomplish  this  a  reference  signal  Is  acquired  from  the  backface 
of  the  sample  or  from  a  flat  surface  on  a  piece  of  similar  material.  This 
signal  is  taken  to  be  the  Impulse  response  of  the  transducer  system  In  the 
absence  of  a  flaw.  The  portion  of  the  acquired  waveform  containing  the  flaw 
signal  Is  separated  by  multiplying  the  waveform  by  a  rectangular  window  whose 
length  and  position  are  selected  by  the  operator.  The  resulting  waveform  1$ 
then  shifted  In  time  so  that  t  ■  0  occurs  at  the  centroid  of  the  pulse.  This 
eliminates  the  excess  phase  that  Is  associate'*  with  a  pulse  that  Is  located 


Preprocessing  of  Acquired  Flaw  Signal 


Fig.  3.2  Preprocessing  of  acquired  flaw  signal. 


SO 


far  from  t  =  0.  This  waveform  is  then  Fourier  transformed  to  obtain  the 
frequency  spectrum  of  the  measured  flaw  signal.  The  reference  (transducer) 
spectrum  signal  is  obtained  in  a  similar  manner. 

The  measured  flaw  spectrum  is  equal  to  the  product  of  the  flaw  spec¬ 
trum  and  the  transducer  spectrum,  plus  various  sources  of  noise.  The  flaw 
spectrum  is  obtained  by  dividing  the  measured  spectrum  by  the  transducer  spec¬ 
trum.  The  division  is  a  noise  vulnerable  operation.  The  energy  in  the  trans¬ 
ducer  spectrum  is  only  above  the  noise  level  over  a  certain  bandwidth.  Out¬ 
side  of  this  bandwidth  a  straight  forward  division  would  produce  a  result  that 
was  dependent  on  the  noise  and  had  undefined  portions  because  of  division  by 
zero.  To  avoid  problems  caused  by  the  noise,  a  desensitized  division  tech¬ 
nique  is  used.  This  is  essentially  the  same  as  the  Wiener  filtering  technique 
described  by  Murakami  et  al.4  This  technique  forces  the  quotient  of  the  two 
spectra  to  zero  when  the  amplitude  of  the  transducer  spectra  becomes  compar¬ 
able  to  that  of  the  noise  associated  with  the  transducer  spectra.  The  Wiener 
filter  has  a  bandpass  amplitude  that  is  near  one  over  most  of  the  bandwidth  of 
the  transducer.  At  the  band  edges  it  decreases  smoothly  to  zero.  It  is 
essentially  zero  at  a  bandwidth  corresponding  to  the  -20  dB  bandwidth  of  the 
transducer  for  the  transducer  spectra  used  in  the  Test  Bed  work. 

The  result  of  this  division  is  an  estimate  of  the  complex  spectral 
back  scattering  amplitude  of  the  flaw  alone.  Frequently  there  are  small 
artifact  signals  produced  outside  the  -20  dB  bandwidth  of  the  transducer  and 
these  are  eliminated  by  replacing  those  portions  of  the  waveform's  spectral 
magnitude  that  lie  below  the  transducer's  bandwidth  with  a  parabolic  curve 
that  passes  through  zero  at  zero  frequency  and  coincides  with  the  flaw 
spectral  magnitude  at  approximately  the  low  frequency  -20  dB  level.  Those 
artifacts  occurring  above  the  bandwidth  of  the  transducer  are  suppressed  by 
multiplying  the  spectral  magnitude  by  a  cosine-squared  function  (Hann  window) 
that  Is  centered  at  zero  frequency  and  goes  to  zero  at  a  frequency  that  Is 
operator  selected  but  usually  corresponds  to  a  frequency  that  Is  20-30*  higher 
than  the  high  frequency  -20  dB  level  of  the  transducer  bandwidth  function. 
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The  exact  value  Is  not  critical  since  the  Born  Inversion  algorithm  tends  to 
suppress  high  frequency  data. 

The  complex  frequency  spectrum  that  results  from  these  calculations 
is  essentially  ready  to  be  Input  to  the  Born  Inversion  algorithm  once  It  has 
been  time  shifted  so  that  the  centroid  of  the  flaw  Is  at  t  =  0.  This  requires 
shifting  the  signal  by  an  amount  t  equal  to  the  difference  between  the 
centroid  of  the  pulse  and  the  centroid  of  the  flaw.  This  quantity  can  best  be 
determined  by  examining  the  long  wavelength  data  obtained  from  the  flaw.  The 
value  of  r  can  then  be  found  from  the  slope  of  the  phase  In  the  long 
wavelength  region.  This  slope  is  found  approximately  by  using  the  a  priori 
Information  that  the  phase  Intercept  at  zero  frequency  must  be  an  Integral 
multiple  of  180  degrees.  The  slope  Is  found  by  fitting  a  straight  line  to  the 
long  wavelength  portion  of  the  phase  curve  of  the  deconvolved  frequency 
spectrum  and  having  the  zero  frequency  intercept  be  an  integral  multiple  of 
180  degrees.  Since  the  slope  of  the  phase  curve  can  be  directly  related  to  a 
time  shift  In  the  time  domain,  the  time  shift  required  to  reduce  the  slope  to 
zero  can  be  calculated.  Generally  this  Is  done  Iteratively  to  produce  the 
best  estimate  of  the  required  time  shift. 

In  some  cases,  the  transducer  used  to  acquire  the  data  may  not  have 
sufficient  bandwidth  to  cover  the  long  wavelength  region.  In  these  cases  a 
lower  frequency  transducer  may  be  used  to  acquire  the  required  long  wavelength 
data.  An  empirical  technique  has  also  been  found  which  works  with  good 
accuracy  when  the  bandwidth  available  is  suitably  broad  and  Is  well  centered 
on  the  first  peak  of  the  magnitude  spectrum.  This  technique  which  depends  on 
certain  features  of  the  characteristic  function  produced  by  the  Born  inversion 
Is  an  Iterative  technique  that  Is  more  easily  explained  after  the  Born 
Inversion  algorithm  has  been  described  and  will  therefore  be  postponed  until 
this  algorithm  Is  presented. 
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3.2.2  Born  Inversion 

The  particular  Born  inversion  technique  that  is  referred  to  in  this 
report  is  the  one-dimensional  inverse  Born  approximation  algorithm  that  was 
developed  by  Rose5  as  part  of  the  Interdisciplinary  Program  for  Quantitative 
NDE.  It  has  been  theoretically  investigated  by  a  number  of  people6  and  its 
operation  has  been  examined  in  both  the  frequency  domain7  and  the  time 
domain.®  The  operation  is  most  easily  explained  in  the  time  domain  and  this 
will  be  done  here.  The  particular  algorithm  that  has  been  used  experimentally 
performs  the  same  processing  on  the  data,  but  for  reasons  of  computational 
convenience  performs  them  in  the  frequency  domain. 

The  theory  of  the  Born  inversion  is  strictly  valid  only  for  weakly 
scattering  or  singly  scattering  objects.  In  practice  it  works  for  many 
strongly  scattering  objects  such  as  voids  and  inclusions  in  metals  and 
inclusions  in  ceramics.  The  operation  Is  illustrated  In  Fig.  3.3.  Consider 
an  ultrasonic  wave  that  is  propagating  toward  a  weakly  scattering  ellipsoidal 
object  as  shown  at  the  top  of  Fig.  3.3.  The  band-limited  impulse  response  of 
this  object  can  be  shown  to  be  a  function  like  that  in  the  center  of  Fig.  3.3 
with  a  small  constant  positive  amplitude  over  the  region  corresponding  to  the 
Interior  of  the  object  and  two  negative  going  peaks  of  short  duration  at  each 
edge  of  the  object.  The  one-dimensional  Inverse  Born  algorithm  generates  a 
characteristic  function  by  Integrating  this  Impulse  response  in  a  particular 
way.  The  Integral  extends  over  a  window  that  Is  symmetrically  located  about 
the  center  of  the  object.  This  center  point  is  taken  as  the  t  ■  0  point.  As 
the  Integration  window  expands.  It  produces  a  function  of  time  (or  distance 
since  the  two  are  related  through  the  velocity  of  sound  in  the  medium).  This 
characteristic  function  Is  shown  schematically  In  the  lower  portion  of  Fig. 
3.3.  The  function  has  zero  slope  at  the  origin  and  then  decreases  to  zero 
amplitude  smoothly  as  the  window  expands.  The  size  of  the  object  along  the 
propagation  direction  of  the  Incoming  wave  is  given  by  the  distance  from  the 
center  of  the  scatterer  to  the  edge  of  the  characteristic  function. 
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dimensional  Bom  inversion 


The  algorithm  used  in  the  Test  Bed  operates  in  the  frequency  domain 
but  performs  all  of  the  equivalent  functions  described  for  the  time  domain 
version.  The  expression  for  the  characteristic  function  is  given  by: 


y(t)  =  y  J  ~  A(<d)  sin  wt  doi 


(3.1) 


where  A(w)  is  the  conditioned  input  spectrum  obtained  from  the  flaw.  The  key 
steps  in  the  algorithm  are  listed  below: 

1.  Divide  spectrum  by  frequency:  A’ (id)  =  A-^-. 

2.  Obtain  inverse  transform  of  A'  (to)  designated  y'(t). 

09 

y1  (t)  =  /  A'  ((D)  sin  u>(t)du> 

.00 

3.  Rescale  the  resulting  time  axis  by  (acoustic  veloc1ty)/2  to 
convert  the  time  units  to  distance  units  measured  from  the 
center  of  the  flaw.  This  radius  coordinate  is  designated  r. 

4.  Divide  the  resulting  function  by  r.  A  linear  extrapolation  of 
the  two  data  points  closest  to  r  *  0  is  used  to  obtain  the 
undefined  point  at  r  ■  0.  The  resulting  function  is  the 
characteristic  function. 

5.  The  radius  estimate  Is  obtained  by  calculating  the  area  under 
the  characteristic  function  and  dividing  It  by  Its  peak 
height.  This  has  been  found  to  work  better  than  a  variety  of 
other  radius  estimators  which  were  investigated. 
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6.  In  the  curves  that  are  plotted,  the  characteristic  function  is 
normalized  through  division  by  the  value  of  the  peak  amplitude 
point.  Further,  since  the  curve  is  symmetric  about  r  =  0,  only 
one  half  is  plotted. 

3.2.3  Testing  of  the  Algorithm 

Much  of  our  time  has  been  occupied  with  the  testing  and  refining  of  the 
Born  inversion  algorithm.  First  the  effect  of  several  parameters  on  the 
accuracy  of  the  algorithm  was  Investigated.  The  parameters  that  have  been 
identified  as  being  of  interest  Include: 

1.  the  signal  to  noise  ratio  of  the  measured  flaw  signal, 

2.  the  ratio  of  the  mid-band  wavelength  of  the  transducer  to  the 

radius  of  the  flaw, 

3.  the  bandwidth  of  the  transducer, 

4.  the  curvature  of  the  wavefronts  Incident  on  the  flaw, 

5.  the  geometric  shape  of  the  flaw. 

To  minimize  the  effects  of  unknown  or  unwanted  paramters,  the  theoretical 
scattering  spectrum  from  a  flaw  was  used  as  an  input  signal  and  the  effects  of 
the  desired  parametric  changes  on  the  predicted  flaw  radius  were  noted. 

To  test  the  sensitivity  of  the  Born  inversion  to  the  presence  of  noise 
In  the  data,  simulated  experimental  waveforms  were  synthesized.  These  waveforms 
consisted  of  the  calculated  scattering  signal  from  a  spherical  void  with  a  noise 
signal  added.  The  noise  signal  was  designed  to  simulate  the  scattering  of 
ultrasound  by  grains  In  the  host  material.  This  type  of  noise  has  a  power 
spectrum  proportional  to  f^  (f  =  frequency)  and  the  noise  waveform  spectrum  Is 
given  by: 

Ng(0  -  N(f)  ,  (3.2) 
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were  fQ  is  a  reference  frequency  and  N(f)  is  a  complex  valued,  white,  Gaussian 
noise  spectrum  consisting  of  samples  of  a  Gaussian  random  variable,  v(f),  for 
which 


E[v(f)  v(f’)*]  =  o*6 ff( 


(3.3) 


Where  E  represents  the  operation  of  averaging  over  the  members  of  the  ensemble, 

2 

v(f)  is  one  member  of  the  noise  ensemble,  a ^  is  the  variance  of  v(f)  and  is 
the  Kronecker  delta.  The  parameter  f0  is  chosen  arbitarily  to  be  the  maximum 
frequency  of  the  sampled  waveform.  The  mean  energy  in  Ng  is  therefore 


'o 

Af  l 
f-f. 


(3.4) 


The  energy  in  the  simulated  signal  is 


U  =  Af  l  |A(f) 
f“-fo 


(3.5) 


where  A(f)  is  the  calculated  spectrum  of  the  signal.  The  signal-to-noise  ratio 
of  the  simulated  experimental  waveforms  is  defined  to  be 


(c/m\2  =  energy  in  flaw  signal 

mean  energy  in  noise  signal 


(3.6) 


(S/N;  *  ^ 


(3.7) 


where  the  bandwidth  extends  from  0  to  f0. 
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The  simulated  waveforms  are  calculated  In  the  frequency  domain,  then 
multiplied  by  the  measured  spectrum  of  an  ultrasonic  transducer  and  finally 
put  into  the  time  domain.  These  operations  can  be  written  in  the  form 

ax(t)  iX(f)  •  [A(f )  +  Ng(f)]}  3.8 

where  ax(t)  is  the  simulated  flaw  waveform,  A(f)  is  the  calculated  spectrum  of 
the  signal,  Ng(f)  is  the  spectrum  of  the  grain  scattering  noise,  X(f)  is  the 
transducer  spectrum  and.^"*  refers  to  a  Fourier  transform  from  frequency  to 
time  domains.  The  Born  inversion  algorithm  is  then  applied  to  the  waveforms 
in  the  same  way  that  it  is  applied  to  experimental  signals. 

Figure  3.4  shows  simulated  and  experimental  flaw  waveforms.  The 
upper  curve  Is  an  experimental  waveform  recorded  for  an  800  pm  diameter 
spherical  void.  The  middle  and  lower  curves  are  simulated  flaw  waveforms. 

The  middle  curve  contains  no  noise  and  the  lower  curve  contains  enough  noise 
to  produce  a  signal  to  noise  ratio  of  0  dB  as  defined  in  Eq.  3.6.  The 
similarity  between  the  simulated  and  the  measured  waveforms  demonstrates  the 
accuracy  of  the  model  used  to  create  the  simulated  waveforms. 

The  transducer  X(f)  used  in  Eq.  3.8  was  assumed  to  be  a  cosine- 
squared  function  that  was  centered  at  5  MHz  and  had  zero  amplitude  at  below 
f  *  0  and  above  10  Miz.  The  calculated  spectrum  was  always  that  of  a 
spherical  void.  Three  different  sizes  of  spheres  were  used.  The  different 
sized  spheres  have  different  portions  of  their  spectrum  within  the  0-10  MHz 
bandwidth  of  the  transducer;  thus  providing  a  test  of  the  sensitivity  of  the 
Born  Inversion  algorithm  to  the  combined  effects  of  noise  and  limited 
bandwidth.  The  radius  estimates  obtained  for  ensembles  of  noisy  signals  as  a 
function  of  signal -to-nolse  ratio  are  shown  In  Figs.  3.5,  3.6  and  3.7.  Both 
the  average  radius  and  the  95%  confidence  limits  are  shown  for  spheres  having 
diameters  of  400  urn,  800  urn  and  1200  urn.  Each  radius  represents  the  average 
of  100  signals  each  composed  of  the  calculated  scattering  amplitude  and  an 


58 


1  I  I  ILL 

l  1  1-1- 

aanindwv 


3anindiAiv 


aanmdwv 


-20  -10  0  10  20  30  40 


S/N  (dB) 

Fig.  3.5  Radius  estimates  for  ensembles  of  noisy  flaw  waveforms  vs  signal  - 
to-nolse  ratio  for  a  200  radius  spherical  void. 
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Fig.  3.6  Radius  estimates  for  ensembles  of  noisy  flaw  waveforms  vs  signal 
to-nolse  ratio  for  a  400  ym  radius  spherical  void. 
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Fig.  3.7  Radius  estimates  for  ensembles  of  noisy  flaw  waveforms  vs  signal  - 
to- noise  ratio  for  a  600  urn  radius  spherical  void. 
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independent  noise  waveform  from  the  ensemble.  Each  of  the  noise  waveforms  has 
the  same  mean  energy.  In  Fig.  3.5,  3.6  and  3.7,  the  dashed  curve  shows  the 
mean  radius  estimate  vs  signal-to-noise  ratio.  The  solid  curves  show  the  95% 
confidence  levels  for  the  ensemble.  As  the  signal-to-noise  ratio  decreases, 
the  uncertainty  of  the  estimates  increases  and  the  mean  of  the  estimates 
eventually  becomes  inaccurate  too.  Note  that  the  95%  confidence  level  is 
within  ±10%  of  the  correct  answer  for  signal-to-noise  ratios  as  low  as  0  dB. 
Examples  of  signals  in  the  presence  of  this  level  of  grain  scattering  noise 
are  shown  in  Fig.  3.4.  The  systematic  decrease  In  the  mean  value  is 
conjectured  to  be  related  to  the  method  of  estimating  the  radius  from  the 
characteristic  function  described  above.  The  area  under  the  characteristic 
function  is  probably  not  significantly  affected  by  the  noise,  but  it  is  likely 
that  the  peak  height  is  increased  by  the  noise  and  causes  the  radius  estimate 
to  be  decreased.  The  values  of  the  mean  estimates  for  high  signal-to-noise 
ratios  are  close  to  the  correct  answer  for  the  800  urn  and  1200  urn  diameter 
spheres  where  the  bandwidth  of  the  transducer  is  encompassing  the  appropriate 
portion  of  the  sphere  scattering  spectrum.  For  the  400  urn  diameter  sphere, 
the  mean  radius  estimate  is  becoming  too  large  and  reflects  the  lack  of 
adequate  high  frequency  data.  The  reason  for  this  will  be  explained  below. 

The  effect  of  transducer  center  frequency  and  bandwidth  on  the  radius 
predictions  provided  by  the  Born  Inversion  can  be  significant.  A  suitable 
technique  for  selecting  a  transducer  was  previously  outlined  In  the  protocol 
section.  To  understand  how  the  bandwidth  affects  the  results,  it  is  helpful 
to  consider  the  magnitude  of  the  scattering  spectrum  from  an  ellipsoidal 
flaw.  The  spectrum  for  a  spherical  void  is  shown  In  Fig.  3.8  where  the 
amplitude  Is  plotted  vs  ka.  The  details  of  the  spectrum  will  change  as  the 
material  properties  of  the  host  material  and  the  flaw  are  varied;  however, 
certain  major  features  will  remain  essentially  unchanged.  In  the  low  ka 
region  the  amplitude  of  the  scattering  always  increases  in  proportion  to  the 
square  of  the  frequency.  The  first  peak  generally  occurs  at  approximately 
ka  *  1;  although  for  certain  Inclusion  types.  It  can  get  as  high  as  ka  =  2. 
Finally  this  first  peak  Is  always  followed  by  many  other  peaks  that  are  not, 
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RELATIVE  AMPLITUDE 


In  general,  equally  spaced  nor  do  their  amplitudes  have  any  simple  relation¬ 
ship  to  one  another.  The  Born  inversion  technique  is  sensitive  to  the  portion 
of  the  spectrum  that  includes  the  first  peak.  It  is  therefore  essential  that 
the  bandwidth  of  the  transducer  used  to  obtain  data  for  the  Born  inversion 
encompass  this  first  peak.  This  requirement  can  become  severe  for  ellipsoidal 
flaws  with  a  factor  of  two  or  more  difference  in  the  lengths  of  the  major  and 
minor  axes.  It  Is  sometimes  necessary  to  resort  to  the  use  of  two  or  more 
transducers  in  these  cases. 

Figure  3.9  shows  radius  estimates  obtained  using  calculated  scatter¬ 
ing  data  for  a  spherical  void  when  the  bandwidth  is  limited.  The  dashed  curve 
shows  the  effects  of  raising  the  minimum  ka  value  present  in  the  data.  The 
result  Is  an  underestimate  of  greater  than  20%  when  kamin  >  0.5.  The  solid 
curve  shows  the  effect  of  lowering  the  maximum  ka  value  present  In  the  data. 
The  result  Is  an  overestimate  of  greater  than  20%  when  kamax  <  2.  From  this 
somewhat  simplified  analysis  one  can  conclude  that  the  bandwidth  required  to 
obtain  radius  estimates  accurate  to  20%  must  extend  from  ka  =  0.5  to  ka  =  2. 

If  insufficient  low  frequency  data  is  available,  then  an  underestimate  is 
obtained  and  if  Insufficient  high  frequency  data  is  available,  an  overestimate 
is  obtained. 

An  ultrasonic  transducer  limits  both  the  high  and  low  frequencies. 
Figure  3.10  illustrates  the  effect  of  having  both  limits  on  the  radius 
estimates.  The  curves  shown  in  Fig.  3.10  are  calibration  curves  obtained  for 
a  particular  type  of  bandpass  window  that  corresponds  to  the  windows  used  to 
generate  Fig.  3.9.  The  high  frequency  limit  was  provided  by  a  cosine  squared 
window  centered  at  f  *  0  and  decreasing  to  zero  at  fmax.  The  low  frequency 
window  was  a  rectangular  window  that  truncated  the  calculated  scattering  data 
at  fmin.  For  this  particular  case  the  results  show  that  as  the  center 
frequency  of  the  transducer  decreases  all  of  the  calibration  curves  coalesce 
Into  one  curve  that  overestimates  the  radius  of  the  flaw.  As  the  center 
frequency  of  the  transducer  Increases,  the  calibration  curve  tends  to  first 
run  parallel  to  a  line  of  constant  radius  and  then  decrease.  An  analysis  of 


65 


S/a 


Fig.  3.10  Normalized  calibration  curves  for  the  Born  Inversion  algorithm  with 
different  bandwldths  and  center  frequencies. 


i 


the  curves  shows  that  they  decrease  below  80%  of  the  true  radius  when  the  low 
frequency  end  of  the  transducer  bandwidth  as  defined  above  exceeds  a  frequency 
corresponding  to  ka  2  0.75.  The  curves  increase  above  120%  of  the  true  radius 
when  the  high  frequency  end  of  the  transducer  bandwidth  defined  above  Is  reduced 
below  a  frequency  corresponding  to  ka  -  1.8.  Thus  for  a  simulated  transducer  we 
obtain  requirements  for  kam^n  and  kamax  that  are  not  too  different  from  the 
crude  estimates  that  were  made  above.  One  should  use  caution  In  applying  these 
results  to  an  experimental  situation  since  the  particular  bandpass  character¬ 
istics  associated  with  the  transducer  and  the  desensitized  deconvolution 
function  must  be  known.  When  these  characteristics  are  known  a  calibration 
curve  can  be  calculated.  This  curve  can  be  used  to  correct  the  estimate 
obtained  from  the  Born  inversion  for  the  errors  caused  by  a  lack  of  bandwidth  or 
because  the  center  frequency  of  the  transducer  Is  too  high  or  too  low. 

In  situations  where  the  true  radius  of  a  scatterer  Is  unknown,  a  set  of 
normalized  curves  like  those  shown  In  Fig.  3.10  are  not  too  useful.  Note, 
however,  that  these  curves  can  be  replotted  In  the  form  shown  In  Fig.  3.11. 

In  this  case.  It  Is  only  necessary  to  know  the  bandwidth  and  center  frequency 
of  the  transducer.  The  entry  point  on  the  vertical  axis  Is  obtained  from  the 
radius  predicted  by  the  Born  Inversion.  The  actual  radius  can  be  found  by 
drawing  a  horizontal  line  from  the  entry  point  to  the  particular  transducer 
curve  that  is  relevant  to  the  experimental  situation,  and  then  drawing  a 
vertical  line  from  this  Intersection  to  the  horizontal  axis.  If  the  inter¬ 
section  with  the  calibration  curve  occurs  In  the  region  where  the  curves  have 
coalesced,  the  center  frequency  of  the  transducer  should  be  Increased. 

Although  It  Is  not  shown  In  Fig.  3.11,  the  curves  all  bend  over  and  drop  to 
low  values  for  sufficiently  high  transducer  center  frequencies.  This 
apparently  occurs  because  the  minimum  frequency  of  the  transducer  Is  no  longer 
encompassing  the  lowest  peak  In  the  magnitude  scattering  spectrum  of  the 
flaw.  The  Born  Inversion  algorithm  Is  beginning  to  Interpret  the  second  peak 
as  If  It  were  the  first  peak.  Presumably,  the  calibration  curves  will  be 
cyclic  as  the  center  frequency  of  the  transducer  Is  raised.  From  an 
experimental  viewpoint,  It  Is  this  aspect  of  the  Born  Inversion  algorithm's 
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operation  that  requires  a  preliminary  look  at  the  magnitude  spectrum  of  the 
flaw  to  insure  that  the  bandwidth  of  the  transducer  j_s_  encompassing  the  first 
peak.  If  this  is  not  true,  the  radius  estimates  provided  by  the  algorithm  are 
likely  to  be  in  error  by  a  factor  equal  to  the  ratio  between  the  frequencies 
of  the  second  and  first  peak  of  the  scattering  spectrum. 

From  this  investigation,  it  can  be  concluded  that  the  bandwidth 
requirements  of  the  Born  inversion  are  not  very  severe.  The  required  ratio  of 
the  maximum  to  minimum  frequencies  is  between  2.4:1  and  4:1  depending  on  the 
exact  nature  of  the  transducer  bandshape.  Since  many  commercially  available 
transducers  have  bandwidths  with  fmax/^n1n  a  can  be  used  ^or  ac(lu^rin9 

scattering  data.  A  single  transducer  can  be  used  for  obtaining  accurate 
radius  estimates  of  flaws  with  radii  variations  of  3.3:1  to  2:1  if  it  has  the 
appropriate  center  frequency.  As  an  example,  a  transducer  with  a  center 
frequency  of  5  MHz,  fmjn  ■  1.1  MHz  and  f,^  =8.9  MHz  would  be  capable  of 
measuring  flaw  radii  In  titanium  with  maximum  values  ranging  from  460  urn  to 
690  vm  and  minimum  values  ranging  from  205  um  to  230  urn.  The  exact  range 
would  of  course  depend  on  the  detailed  shape  of  the  transducer  bandpass  and 
the  slgnal-to-nolse  ratio. 

The  next  area  of  concern  about  the  performance  of  the  Born  Inversion 
technique  involves  the  distortion  of  the  Incident  wave  by  the  entry  surface  of 
the  part.  This  can  occur  when  flaws  are  located  below  curved  surfaces  or  when 
the  Incident  wave  must  propagate  obliquely  through  a  flaw  surface.  There  are 
two  potential  problems  Involved.  First,  the  theory  for  the  Born  Inversion 
assumes  that  a  plane  wave  Is  Incident  on  the  flaw.  It  is  not  known  how  the 
accuracy  Is  affected  If  this  wavefront  becomes  curved  or  otherwise  distorted. 
Second,  the  deconvolution  technique  uses  a  reflection  off  a  flat  surface  as  a 
reference.  In  the  first  case,  the  focusing  or  distortion  of  the  wave  at  the 
flaw  can  In  principle  be  removed  via  a  backward  wave  propagation  technique; 
although  In  practice,  this  Is  not  done.  In  the  second  case,  the  consequences 
of  the  deviation  of  the  wavefront  from  a  plane  wave  on  the  character  of  the 
deconvolved  waveform  are  not  known. 
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The  errors  caused  by  these  wavefront  distortions  can  be  calculated 
for  specific  geometries  and  experimental  situations,  but  this  was  beyond  the 
scope  of  this  program.  Furthermore,  this  calculation  did  not  seem  essential 
to  the  primary  goal  of  determining  how  the  inversion  technique  would  work  in 
practical  situations.  In  lieu  of  a  theoretical  calculation,  an  empirical 
approach  has  been  used  to  evaluate  errors  that  are  encountered  with  some 
common  geometries.  The  reference  waveform  is  typically  the  signal  obtained 
from  reflection  at  normal  incidence  of  a  wave  from  the  backface  of  the 
specimen  being  examined.  In  the  cases  that  have  been  tried,  the  backface  is 
about  10  to  20  wavelengths  further  along  the  path  of  propagation  than  the 
flaw.  In  a  few  cases,  a  backface  echo  from  a  separate  specimen  made  of  the 
same  material  has  been  used  as  a  reference.  The  signal  waveforms  back- 
scattered  from  the  flaw  have  been  obtained  using  both  waves  that  are  incident 
normally  on  the  front  surface  of  the  specimen  and  waves  that  are  incident  at 
various  angles,  with  the  maximum  angle  corresponding  to  an  internal  refraction 
angle  of  60°.  In  all  cases,  the  deconvolution  has  been  computed  using  the 
backface  echo  obtained  from  the  normally  incident  wave.  Both  the  backface  and 
the  flaw  were  in  the  far  field  of  the  transducer.  A  number  of  flaws  having 
known  sizes  were  examined  and  the  data  was  analysed  using  the  Born  inversion 
algorithm.  There  was  no  evidence  that  there  were  any  errors  being  generated 
via  the  deconvolution  that  could  be  attributed  to  wavefront  distortion.  Some 
backscattered  signals  have  also  been  acquired  through  a  concave  cylindrically 
curved  interface  with  a  radius  of  curvature  of  3.5  inches.  The  data  has  not 
been  exhaustively  tested  but  the  Initial  indications  are  that  an  underestimate 
of  between  10%  and  25%  is  obtained.  Although  these  errors  are  not  considered 
significant,  the  possibility  that  curved  wavefronts  could  cause  systematic 
errors  should  be  considered  when  using  the  Born  inversion.  In  these  tests, 
the  reference  waveform  was  obtained  from  a  flat  surface.  No  investigation  was 
made  of  the  effect  of  focused  beams  on  the  accuracy  of  the  Born  inversion. 
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Experimental  Results 


Experiments  involving  the  acquisition  of  data  from  flaws  and 
requiring  the  use  of  various  aspects  of  the  protocol  have  been  conducted 
concurrently  with  the  other  investigations.  Most  of  the  experiments  have  been 
conducted  using  known  flaw  types  that  have  been  machined  into  titanium  pieces 
which  were  subsequently  diffusion  bonded.  This  procedure^  for  fabricating 
samples  with  known  defects  was  developed  at  the  Science  Center  and  has  been 
extensively  used  to  produce  a  variety  of  flaw  types  including: 

1.  Ellipsoidal  voids 

2.  Spherical  inclusions 

3.  Simulated  fatigue  cracks. 

A  set  of  samples  were  specifically  prepared  for  this  program  that 
were  designed  to  simulate  turbine  bore  samples.  These  samples  are  described 
in  Appendix  A.  A  sample  was  also  obtained  from  Rolls-Royce  Ltd.  that 
contained  a  naturally  occurring  flaw.  This  sample  provided  an  example  of  the 
performance  of  the  Born  inversion  algorithm  when  the  characteristics  of  the 
flaw  were  unknown. 

3. 2. 4.1  Description  of  Experiments 

Much  of  the  data  was  acquired  from  spherical  flaws  located  in 
titanium  disks.  The  disks  are  4  in.  in  diameter  and  1  in.  thick.  They  are 
made  from  extruded  titanium  bar  stock.  The  apparent  grain  size  looking  along 
the  extrusion  direction  is  relatively  small  and  creates  a  low  amplitude 
ultrasonic  grain  scattering  noise.  In  the  transverse  direction,  the  grains 
are  considerably  larger  and  the  grain  scattering  noise  level  is  significant 
compared  to  the  signal  levels. 

Experimental ly  the  first  problem  is  the  detection  of  the  flaw. 
Although  there  are  a  variety  of  techniques  for  enhancing  the  reliability  of 
detection  through  detailed  analysis  of  a  t^tative  flaw  indication,  the  first 
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task  is  to  determine  that  there  is  something  present  in  the  waveform  to  be 
analyzed.  This  means  that  the  signal-to-noise  ratio  must  be  sufficiently  high 
that  the  amplitude  of  the  signal  is  above  the  noise  for  part  of  the  time  that 
it  is  present.  Therefore  the  initial  detection  is  accomplished  using  an 
amplitude  threshold  technique. 

It  is  desirable  to  obtain  as  high  a  signal  to  noise  ratio  as  pos¬ 
sible.  This  can  be  achieved  by  either  increasing  the  signal  level  or 
decreasing  the  noise  level.  The  major  noise  sources  are  either  thermal  noise 
generated  in  the  electronic  amplifier,  grain  scattering  noise  in  the  material 
being  investigated,  or  residual  transducer  "ring  down"  signals  from  the  front 
face  of  the  part.  The  thermal  noise  can  be  reduced  through  the  time  averaging 
of  multiple  signals.  The  grain  scattering  noise  can  in  principle  be  reduced 
through  the  spatial  averaging  of  multiple  signals,  but  this  has  not  been 
achieved  in  practice.  The  transducer  ’’ring-down"  noise  can  be  reduced  by 
subtracting  two  signals  acquired  from  different  spatial  regions  of  the  part 
but  having  identical  front  surface  echos.  For  the  Test  Bed  work,  the  time 
averaging  of  waveforms  and  subtraction  have  been  used  where  necessary  to 
reduce  thermal  and  "ring-down"  noise  below  the  level  of  the  grain  scattering 
noise.  Subtraction  has  been  able  to  suppress  front  surface  "ring-down"  noise 
by  25-30  dB.  This  makes  the  grain  scattering  noise  the  dominant  source  of 
noise  in  most  of  the  experimental  work.  The  strength  of  this  noise  source 
increases  as  the  square  of  the  frequency  and  this  must  be  considered  when 
selecting  a  transducer. 

The  other  approach  to  enhancing  the  signal-to-noise  ratio  is  to 
increase  the  signal  level.  If  grain  scattering  noise  is  dominant,  simply 
increasing  the  output  of  the  transducer  will  not  help,  since  the  amplitude  of 
the  grain  scattering  noise  rises  In  proportion  to  the  strength  of  the  incident 
signal.  The  signal  level  can  be  increased  by  focusing  the  transducer  provided 
the  size  of  the  defect  is  less  than  the  focal  spot  size  of  the  beam.  A 
decrease  of  the  spot  size  diameter  by  a  factor  of  two  increases  the  signal-to- 
noise  ratio  by  12  dB.  Since  the  minimum  spot  size  is  dependent  on  the 


wavelength,  higher  frequency  transducers  will  produce  smaller  spots  and 
correspondingly  higher  signals.  Note  that  this  frequency  dependence  and  that 
of  the  grain  scattering  tend  to  cancel.  In  some  cases  the  flaw  is  oriented  in 
an  unfavorable  way  and  the  backscattered  signal  amplitude  is  weak.  The 
situation  can  sometimes  be  improved  by  tilting  the  transducer  relative  to  the 
flaw  and  thereby  obtaining  a  larger  backscattered  signal. 

The  four  inch  diameter  titanium  disk  samples  allow  high  signal  to 
noise  ratio  signals  to  be  obtained.  Consequently,  the  samples  containing 
400  pm,  800  pm  and  1200  pm  diameter  spherical  voids  have  been  extensively  used 
to  investigate  experimentally  the  effects  of  transducer  bandwidth  on  the 
accuracy  of  the  Born  Inversion.  Although  the  flaw  location  in  these  samples 
is  known,  a  search  mode  with  a  C-scan  image  output  was  employed  using  the 
techniques  outlined  above  to  enhance  flaw  detectability.  An  example  of  one  of 
these  is  shown  in  Fig.  2.15. 

As  explained  in  Section  3.2.1,  the  signals  for  the  Born  inversion  are 
deconvolved  to  remove  the  transducer  dependent  components.  In  principle  the 
reference  signal  used  in  the  deconvolution  should  be  obtained  from  a  plane 
surface  that  is  the  same  distance  from  the  front  face  of  the  part  as  the 
flaw.  In  practice,  it  has  been  determined  that  a  reference  signal  from  the 
back  face  of  the  part  or  from  a  plane  surface  in  another  part  made  from  the 
same  material  works  satisfactori ly.  This  allows  a  single  reference  signal  to 
be  used  for  the  deconvolution  of  signals  acquired  at  different  angles  relative 
to  the  flaw.  There  is  some  question  about  the  problem  that  may  occur  for 
flaws  located  below  curved  surfaces;  i.e.,  is  it  necessary  to  obtain  a 
reference  signal  from  a  plane  below  the  same  curved  surface?  For  modest 
curvatures  the  errors  introduced  appear  to  be  small;  however,  an  extensive 
study  of  the  effect  was  not  conducted. 

When  the  Born  Inversion  algorithm  Is  used  to  obtain  the  size  of  a 
flaw,  the  location  of  the  center  of  the  flaw  has  to  be  precisely  determined. 
The  centroid  of  the  pulse  that  Is  backscattered  from  the  flaw  Is  used  as  an 
Initial  estimate  of  the  center  of  the  flaw  and  Is  sufficiently  accurate  to 
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obtain  a  deconvolved  waveform  during  the  pre-processing  of  the  signal  as 
indicated  in  Section  3.2.1.  Two  techniques  have  been  used  to  obtain  a  more 
precise  estimate  of  the  flaw  center.  The  first  technique  requires  that  the 
characteristic  function,  defined  in  Section  3.2.2,  have  a  flat  top  when  it 
intersects  the  axis  at  r  =  0.  This  is  an  empirical  technique  that  does  not 
have  a  strong  theoretical  basis.  It  does  provide  accurate  answers  in  cases 
where  the  bandwidth  of  the  transducer  is  centered  on  the  portion  of  the  flaw 
scattering  spectrum  containing  the  first  peak.  There  are  two  steps  involved 
in  obtaining  the  time  shift.  First  the  data  is  analysed  using  the  Born 
inversion  algorithm  with  a  range  of  time  shifts  both  positive  and  negative 
relative  to  the  centroid  of  the  pulse.  This  results  in  a  curve  of  the  type 
shown  in  Fig.  3.12.  The  fact  that  some  of  the  excursions  of  the  curve  are 
positive  and  some  are  negative  is  not  relevant.  The  sign  is  dependent  on  the 
pre-processing  of  the  signal  and  either  sign  is  acceptable.  The  important 
characteristic  of  this  curve  is  the  shape.  Note  that  in  the  region  from 
t  *  -0.05  to  «  0.15,  the  curve  has  a  sharp  minima,  a  broad  local  maxima,  and  a 
second  minima  that  is  less  sharp  than  the  first.  This  shape  forms  a  recogniz¬ 
able  signature  and  It  is  known  that  the  center  of  the  flaw  is  located  near  the 
second  minima  and  displaced  to  the  left  of  it.  The  exact  location  is  determined 
by  generating  a  sequence  of  characteristic  functions  for  time  shifts  ranging 
from  t  *  0  to  t  *  0.1  and  recording  which  curve  has  a  maximally  flat  top  at  the 
origin.  The  radius  estimate  associated  with  this  curve  is  taken  as  the  radius 
of  the  flaw.  The  curve  in  Fig.  3.12  was  generated  using  data  obtained  from  an 
800  pm  diameter  spherical  void  with  a  nominal  5  MHz  transducer.  Note  that  the 
radius  estimates  between  the  broad  local  maxima  and  the  second  minima  range  from 
•»  270  pm  to  -  390  pm.  Since  the  correct  value  is  400  pm,  the  lowest  estimate  is 
only  in  error  by  about  32%.  For  many  applications  this  may  be  acceptable,  in 
which  case  the  determination  of  the  maximally  flat  characteristic  function  is 
unnecessary.  Some  caution  should  be  used  with  this  latter  approach  since  the 
variation  of  the  maximum  error  as  the  flaw  spectrum  shifts  within  the  transducer 
bandwidth  has  not  been  Investigated.  Furthermore,  Inclusions  and  other  flaw 
types  may  produce  different  results. 

The  second  technique  for  determining  the  center  of  the  flaw  is 
rigorously  correct.  It  Involves  the  use  of  long  wavelength  data  to  determine 


FLAW  SIZE  (mm) 


TIME  SHIFT  (//sec) 


Flq.  3.12  Born  Inversion  estimates  of  flaw  radii  vs  time  shifts  of  the  flaw 
response  relative  to  Its  centroid. 
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the  phase  of  the  deconvolved  flaw  signal.  When  the  time  corresponding  to  the 
center  of  the  flaw  is  correctly  chosen,  the  phase  at  low  frequencies  (corres¬ 
ponding  to  ka  <  0.5)  is  constant  and  has  a  value  that  is  an  integral  multiple  of 
n.  This  method  was  described  in  Section  3.2.1 

3. 2. 4. 2  Spherical  Flaws 

Experimental  data  was  acquired  from  the  spherical  voids  on  many  occa- 
ions  and  using  several  different  transducers.  There  are  two  sets  of  data  that 
will  serve  to  summarize  the  results.  One  set  of  data  was  acquired  from  an  800 
urn  diameter  spherical  void  using  the  transducer  placement  pattern  shown  in  Fig. 
3.13b.  The  geometry  and  definition  of  the  angles  0^  and  <{>i  are  shown  in  Fig. 
2.13.  Note  that  the  angles  specified  in  Fig.  3.13b  are  internal  angles.  The 
signals  were  acquired  using  a  5  Wz  transducer  with  a  -20  dB  full  bandwidth  of 

5.50  MHz.  The  results  are  given  in  Table  3.1.  The  mean  diameter  is  390  un  and 

the  error  bars  to  the  95%  confidence  level  extend  ±  30  urn  about  this  mean. 

The  second  set  of  data  was  acquired  from  a  1200  urn  diameter  spherical 
void  using  the  transducer  placement  pattern  shown  in  Fig.  3.13a.  The  geometry 
and  definition  of  the  angles  are  the  same  as  for  the  800  pm  void.  In  this  case 
the  signals  were  acquired  using  a  2.25  MHz  transducer  with  a  -20  dB  full  band¬ 
width  of  2.65  MHz.  The  three  plots  shown  in  Fig.  3.14  summarize  the  results. 

The  Born  Inversion  assumes  that  the  object  measured  has  inversion  symmetry. 

This  Implies  that  the  diameter  Is  being  measured;  thus  each  measurement  produces 
two  points  on  the  polar  plots  in  Fig.  3.14.  To  show  a  full  cross  section,  the 
data  obtained  for  two  values  of  4>^  are  combined  in  each  of  the  diagrams.  Thus 
the  first  cross  section  includes  data  acquired  at  4>-|  =  0°  and  180°  for  0^  =  0, 
30°,  45°,  and  60°.  Similarly  the  data  at  4>^  =  120°  and  300°  is  combined  in  the 

second  cross  section  and  that  for  4^  =  60°  and  240°  in  the  third.  The  cross 

sections  shown  correspond  to  4>-|  =  0,  60°  and  120°  and  0^  =  0°,  ±30°,  ±45°,  and 
±60°.  The  mean  diameter  Is  1320  um  and  the  error  bars  to  the  95%  confidence 
level  extend  ±  210  um  about  this  mean.  The  mean  diameter  Is  approximately  10% 
high.  About  half  of  this  error  can  be  attributed  to  the  bandpass  of  the 
transducer  being  centered  somewhat  below  the  optimum  frequency,  while  the 
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urn  diameter  flaw. 


Table  3.1 

Results  of  Born  Inversion  Analysis  of  An 
800  um  Diameter  Spherical  Void 


e1 

(degrees) 

*1 

(degrees) 

Estimated 

Radius 

(um) 

0 

0 

390 

30 

30 

390 

30 

90 

390 

30 

150 

390 

30 

210 

390 

30 

270 

30 

330 

380 

52 

0 

420 

52 

60 

380 

52 

120 

390 

52 

180 

390 

52 

240 

440 

remaining  error  can  be  accounted  for  by  uncertainties  In  the  actual  size  and 
shape  of  the  defect. 

3. 2. 4. 3  Oblate  Spheroids 

A  test  of  the  Born  Inversion  algorithm  with  an  ellipsoidal  void  was 
more  difficult  because  of  the  wide  bandwidth  required  and  the  need  to  acquire 
data  with  two  transducers  having  different  center  frequencies.  Several 
attempts  were  made  to  acquire  data  from  ellipsoidal  voids  In  samples  with  flat 
entry  surfaces.  These  attempts  were  made  prior  to  the  determination  of  the 
bandwidth  requirements  of  the  Born  Inversion  algorithm  which  was  discussed  In 
Section  3.2.3.  Furthermore,  good  techniques  for  experimentally  determining 
the  location  of  the  center  of  the  flaw  had  not  been  established. 
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In  most  cases,  the  estimate  of  the  size  of  the  flaw  for  normal  inci¬ 
dence  of  the  ultrasonic  beam  was  accurate;  but  for  other  angles  the  results 
were  inaccurate.  This  inaccuracy  is  now  attributed  to  the  inadequate  band¬ 
width  in  some  cases  and  inadequate  slgnal-to-noise  ratio  in  others.  To 
demonstrate  that  the  Born  inversion  technique  could  be  used  with  ellipsoidal 
flaw  shapes,  a  somewhat  simpler  experimental  situation  was  selected.  A 
400  urn  x  800  pm  oblate  ellipsoid  was  available  in  a  "trailer  hitch"  sample  and 
the  contact  transducer  available  for  use  with  it  had  a  center  frequency  of  6.5 
MHz  and  a  -20  dB  full  bandwidth  of  13.5  MHz.  Thus  a  single  transducer  was 
adequate  for  measurements  at  all  angles.  A  set  of  measurements  was  made  for 
polar  angles  ranging  from  zero  degrees  to  90  degrees.  The  results  are 
displayed  in  Fig.  3.15. 

The  solid  curve  shows  the  outline  of  the  oblate  ellipsoid.  The 
dotted  curve  with  the  hourglass  shape  shows  the  results  that  would  be  produced 
by  a  true  Born  scatterer;  i.e.,  the  diameter  obtained  at  each  angle  would  be 
equal  to  the  distance  between  two  planes  that  were  perpendicular  to  the 
direction  of  propagation  and  tangent  to  the  near  and  far  surfaces  of  the 
ellipsoid.  The  complex  spectra  of  the  backscattered  signal  from  a  2:1  oblate 
ellipsoid  were  calculated  by  Opsal^  at  Lawrence  Livermore  Laboratories  and 
furnished  to  the  Science  Center.  This  theoretical  data  was  used  as  ?'»  input 
to  the  Born  Inversion  algorithm  and  the  results  are  given  by  the  open 
circles.  Note  that  the  resulting  curve  is  roughly  elliptical  In  shape  but  is 
slightly  larger  than  either  the  solid  or  dotted  curves.  The  cause  of  the 
discrepancy  has  not  been  identified. 

Finally,  the  experimental  results  are  given  by  the  open  squares  and 
open  triangles.  The  open  squares  represent  results  obtained  when  the  estimate 
of  the  location  of  the  center  of  the  ellipsoid  is  obtained  using  long 
wavelength  techniques  as  explained  In  Section  3.2.1.  The  open  triangles 
represent  results  obtained  when  the  location  of  the  center  of  the  ellipsoid  is 
obtained  using  the  empirical  technique  of  requiring  the  characteristic 
function  to  have  a  flat  top  at  the  origin  as  explained  above.  Note  that  these 
experimental  results  span  the  range  of  results  obtained  by  the  other 
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techniques.  The  overall  accuracy  Is  good  and  the  shape  is  definitely 
elliptical.  This  is  believed  to  be  the  first  time  that  extensive  measurements 
made  on  an  ellipsoid  of  revolution  and  used  with  an  inversion  technique  have 
resulted  in  accurate  size  estimates  of  the  ellipsoid. 

These  results  demonstrate  that  experimental  procedures  have  been 
devised  to  obtain  accurate  results  with  the  Born  inversion  algorithm.  No 
major  problems  are  being  encountered  with  the  deconvolution  procedures  or 
other  pre-processing  steps.  Further  the  results  demonstrate  that  the  Born 
inversion  can  be  used  to  obtain  the  shape  of  the  flaw.  In  the  cases  presented 
thus  far  the  shape  of  the  flaw  was  known  a  priori.  To  establish  the  validity 
of  the  technique,  it  is  important  to  obtain  the  characteristics  of  a  flaw 
whose  shape  is  unknown  before  the  experiment. 

3. 2. 4. 4  Rolls-Royce  Sample 

A  sample  made  of  nowder  metallurgy  nickel  based  superalloy  material 
was  obtained  from  Rolls-Royce  Ltd.  This  sample  is  shown  in  Fig.  3.16.  The 
sample  was  known  to  contain  a  flaw  but  none  of  the  characteristics  were 
known.  The  objective  was  to  apply  the  full  range  of  quantitative  flaw 
characterization  techniqes  to  this  sample  and  predict  the  characteristics  of 
the  flaw  as  completely  as  possible.  Upon  completion  of  this  task,  the  sample 
was  returned  to  Rolls-Royce,  Ltd.  and  mecallographically  sectioned  so  that  the 
true  characteristics  of  the  flaw  could  be  determined. 

The  protocol  described  in  Section  3.1  was  used  to  inspect  the 
sample.  Because  of  the  small  sample  size,  the  search  mode  was  combined  with 
the  imaging  that  is  required  as  the  first  step  of  the  detailed  examination. 

The  C-scan  and  B-scan  images  obtained  of  the  flaw  are  shown  in  Fig.  3.17. 

From  these  Images,  the  orientation  and  long  dimension  of  the  flaw  were 
determined.  The  Images  were  also  used  to  select  the  observation  directions 
for  the  backscattering  measurements  that  provide  input  data  for  the  Born 
Inversion  analysis.  This  data  was  analysed  and  the  short  transverse  dimen¬ 
sions  of  this  needle-like  flaw  were  estimated. 


Fig.  3.16  Rolls-Royce  sample  made  of  powder  metallurgy  nickel  based 
superalloy  material . 
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Fig.  3.17  (a)  C-scan  Image  of  flaw  In  Rolls-Royce  sample . 

(b)  B-scan  image  of  flaw  In  Rolls-Royce  sample. 
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Finally  the  phase  of  the  signals  returned  from  the  flaw  were  examined 
to  determine  that  the  acoustic  Impedance  of  the  flaw  was  greater  than  that  of 
the  surrounding  material.  These  measurements,  summarized  In  Table  3.2,  were 
returned  to  Rolls  Royce  along  with  the  sample. 


Table  3.2 

Results  of  Examination  of  Rolls-Royce  Sample 

X  Axis  *  0.172  Inch  (Imaging) 

Y  Axis  «  0.007  Inch  Born  Inversion 

Z  Axis  »  0.006  inch  Born  Inversion 

Angle  to  X  Axis  ■  4.3° 

Acoustic  Impedance  of  defect  >  metal 


The  shape  of  the  flaw  as  determined  by  metal lographlc  sectioning  Is 
shown  In  Fig.  3.18.  When  the  sectioning  was  being  done  the  presence  of  a 
prior  particle  boundary  region  surrounding  the  flaw  was  detected.  The  flaw 
was  a  calcium,  silicon,  aluminum,  titanium  rich  ceramic  material  and  the 
boundary  region  represents  a  ceramic/matrix  reaction  zone.  The  flaw 
dimensions  are  all  specified  In  Fig.  3.18.  Note  the  similarity  of  the  flaw 
shape  and  of  the  C-scan  Image  In  Fig.  3.17a  even  though  the  transverse 
dimension  Is  not  resolved  In  the  C-scan.  The  results  of  the  quantitative  NDE 
analysis  and  the  metal lographlc  sectioning  are  compared  In  Table  3.3.  The 
agreement  between  the  two  techniques  Is  excellent.  These  results  are  far  more 
precise  than  those  obtained  using  conventional  NDE  procedures. 

For  calibration  purposes  Rolls-Royce  measured  the  amptltude  of  the 
signal  returned  from  this  flaw  relative  to  the  amplitude  of  the  signal 
returned  from  a  0.050  Inch  diameter  flat  bottom  hole.  The  measurement  was 
repeated  using  three  different  transducers  and  two  different  pul ser/ receiver 
units.  For  each  transducer  and  pul ser/ receiver  combination,  the  measurement 
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Fig.  3.18  Results  of  metal lographic  sectioning  of  flaw. 


Table  3.3 

Comparison  of  Test  Bed  Results  and  Metallographlc 
Examination  Results 


Science  Center 
Test  Bed 

Rolls-Royce  Sectioning 

Flaw 

Flaw  Reaction  Zone 

X  Axis  0.172  Inch 

Y  Axis  0.007  inch 

Z  Axis  0.006  Inch 

Angle  4.3° 

0.161  0.178 

0.011  0.036 

0.002  0.009 

4°  -  6° 

Acoustic  Impedance  higher  than  the  surrounding  metal  Is 
consistent  with  a  ceramic  inclusion. 


was  made  from  each  of  two  opposite  farces  of  the  sample.  The  measured 
amplitudes  varied  by  5  dB.  Since  the  conversion  factor  between  flaw  size  and 
scattering  amplitude  Is  dependent  on  flaw  shape,  orientation  and  composition 
as  well  as  flaw  size,  one  cannot  precisely  state  the  flaw  size  variation  that 
would  correspond  to  a  5  dB  amplitude  variation.  If  the  shape,  orientation  and 
composition  remain  constant,  the  amplitude  will  be  proportional  to  flaw  size 
for  flaws  with  ka  greater  than  approximately  one.  Thus,  the  5  dB  variation 
would  correspond  to  nearly  a  factor  of  two  variation  In  flaw  size. 


4.0  PHASED  ARRAY  SYSTEM 


4.1  Concept  and  Purpose 

Even  though  the  use  of  phased  ultrasonic  array  systems  Is  widespread 
In  the  medical  community11"16  there  Is  currently  little  use  of  them  for  NDE 
applications.  There  Is  Interest  In  using  them  but  appropriate  applications 
must  be  found.  It  Is  also  Important  to  realize  that  the  requirements  placed 
on  an  array  for  medical  applications  are  quite  different  than  those  needed  for 
NDE  applications.16  The  Intent  of  the  design  of  this  system  Is  that  It  be 
useful  In  a  number  of  configurations  so  that  many  different  applications  can 
be  Investigated.  For  this  reason  the  system  Is  not  very  compact.  This  should 
not  be  viewed  as  a  shortcoming  since  the  essential  features  for  a  particular 
application  can  always  be  abstracted  to  produce  a  compact  system.  The  two 
Important  questions  are:  Is  the  system  reaching  Its  potential  as  a  data 
acquisition  device?  And,  how  does  Its  performance  compare  with  alternative 
techniques? 

The  usefulness  of  an  array  for  NDE  applications  stems  from  the 
capability  of  the  array  to  rapidly  scan  an  ultrasonic  beam  and  to  change  the 
focus  to  accommodate  the  curved  surfaces  of  different  parts.  This  latter 
characteristic  will  be  referred  to  as  beam  agility.  At  the  present  time,  the 
applications  of  arrays  In  NDE  are  limited.  These  limitations  are  due  to  the 
current  technology  relating  to  array  fabrication  and  to  physical  constraints 
Imposed  by  materials  that  are  to  be  Inspected.  A  thorough  discussion  of  this 
aspect  of  the  problem  has  been  given  In  the  paper  by  Addison.16 

The  rapid  beam  scanning  capability  of  an  array  can  be  exploited  In 
several  different  ways  In  a  quantitative  NDE  application.  If  the  array  Is 
sufficiently  long.  It  can  be  used  In  a  contour  scanning  mode  for  Inspecting 
circularly  symmetric  parts  such  as  engine  disks.  An  example  of  this  mode  of 
operation  Is  shown  In  Fig.  4.1.  Here  the  beam  Is  always  entering  the  part 
normal  to  the  surface.  The  beam  agility  feature  can  be  used  to  compensate  for 
focusing  effects  caused  ty  curvature  of  the  metal  surface.  The  tilting  and 
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pig.  4.1  Contour  scanning  mode. 
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translation  of  the  beam  Is  done  electronically  and  avoids  problems  with 
mechanical  Inaccuracies  and  backlash  that  are  frequently  encountered  In 
mechanical  drives.  The  technological  capability  for  making  a  linear  array 
transducer  that  would  be  long  enough  to  scan  over  the  entire  radius  of  a 
turbine  disk  has  not  been  demonstrated  at  this  time.  At  the  Inception  of  this 
program.  It  appeared  likely  that  an  array  transducer  of  length  6  Inches  with 
the  performance  needed  for  this  application  could  be  purchased.  This  did  not 
prove  to  be  the  case.  To  assure  the  performance  of  the  array.  It  was  neces¬ 
sary  to  reduce  the  length  to  about  1/2  Inch.  To  preserve  some  of  the  capa¬ 
bilities  of  a  longer  array  transducer,  two  of  the  shorter  array  transducers 
are  used  with  a  mechanically  variable  spacing  between  them  as  shown  In 
Fig.  4.2.  Although  this  severely  limits  the  use  of  the  array  system  for 
contour  following  applications.  It  can  be  used  to  demonstrate  the  feasibility 
of  the  concept. 

If  the  use  of  the  array  system  In  a  contour  following  mode  is 
considered  the  equivalent  of  a  search  mode  for  a  single  element  transducer, 
then  the  use  of  the  array  In  a  sector  scan  mode  can  be  considered  as  the 
equivalent  of  a  detailed  examination  mode.  The  array  can  be  set  up  to 
generate  a  fan  beam  to  obtain  a  B-scan  Image  of  the  flaw  or  It  can  be  used  to 
obtain  scattering  data  from  the  flaw  that  can  be  used  with  the  Born  Inversion 
or  long  wavelength  techniques.  This  scattering  mode,  shown  In  Fig.  4.3,  can 
acquire  data  using  either  a  pulse-echo  or  a  pitch-catch  technique. 


4.2  Description  of  System 

The  electronics  for  driving  this  array  have  a  somewhat  different 
objective  than  some  of  the  array  based  systems  that  are  currently  available, 
which  acquire  and  display  an  Image  In  as  little  time  as  possible.  This  system 
must  be  sufficiently  flexible  to  be  used  In  a  wide  variety  of  situations  and 
be  capable  of  operating  In  many  different  modes.  It  must  acquire  a  waveform 
while  minimizing  distortion  due  to  crosstalk  or  nonllnearltles.  It  Is 
desirable  to  display  a  single  Image  and  be  able  to  recognize  the  Important 
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details  of  the  object  under  study.  The  system  must  process  the  signals  with 
sufficiently  low  distortion  that  when  scattering  data  Is  acquired,  the 
resulting  signals  can  be  used  with  the  various  inversion  techniques. 

The  block  diagram  In  Fig.  4.4  delineates  the  major  subassemblies  of 
the  array  system  and  their  functions.  The  S/200  minicomputer  sends  out  a  set 
of  codes  that  selects  the  desired  set  of  16  transmit  elements,  the  desired  set 
of  16  receive  elements,  and  the  set  of  time  delays  that  define  the  transmit 
beam  direction.  These  codes  are  stored  In  the  control  memory.  The  timing  and 
control  block  Interprets  the  codes  and  activates  the  transmit  and  receive 
multiplexers  so  that  the  appropriate  array  elements  are  selected.  This  block 
also  sends  suitably  delayed  trigger  signals  to  each  of  the  16  selected  trans¬ 
mitter  elements  to  synthesize  an  ultrasonic  beam  with  the  desired  direction. 
The  trigger  signals,  which  are  positive  going  TTL  level  pulses,  activate  the 
pulser  and  cause  a  5  to  400  volt  pulse  to  be  applied  to  each  selected 
transducer  element.  The  returning  ultrasonic  signal  from  the  object  under 
Investigation  Is  received  by  each  of  the  16  selected  receive  elements.  The 
received  signals  are  each  amplified  26  dB  by  a  low  noise,  low  distortion, 
hybrid  preamplifier.  They  are  then  passed  through  the  receiver  multiplexer  to 
a  line  driver  and  then  to  a  signal  conditioning  circuit.  The  signal  con¬ 
ditioning  circuit  serves  as  an  anti-aliasing  filter  and  adjusts  the  peak-to- 
peak  amplitude  of  the  signal  as  well  as  Its  dc  level  to  make  It  compatible 
with  the  A/D  converter.  The  next  block  Is  the  8  bit,  18  MHz  A/D  converter 
which  digitizes  the  signal  and  loads  It  Into  an  8  x  l  k  bit  fast  random  access 
buffer  memory.  All  16  of  the  A/D  converters  are  started  synchronously  after  a 
delay  corresponding  to  an  ultrasonic  round  trip  time  that  Is  selected  by  the 
operator.  The  contents  of  each  of  the  sixteen  memories  are  sequentially 
clocked  through  the  Interface  to  the  array  processor  at  a  rate  compatible  with 
the  operation  of  the  arrqy  processor  and  minicomputer  (approximately 
0.5  Wz).  The  array  processor  has  the  task  of  shifting  the  waveforms  and 
summing  them  to  synthesize  a  beam  from  a  specified  direction.  The  array 
processor  also  Interpolates  between  data  points  In  the  waveform  If  the  desired 
shift  Is  not  an  Integral  number  of  clock  cycles.  This  reduces  the  sldelobe 
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Fig.  4.4  Ultrasonic  phased  array  system  block  diagram. 


levels  In  the  array  response.  The  array  processor  is  also  used  to  correct  for 
amplitude  errors  In  the  response  of  a  particular  array  element.  After  the 
beam-forming  operation  Is  completed,  the  processed  signals  are  sent  to  the 
S/200  minicomputer  for  storage  on  the  disk  memory  and  to  be  displayed  on  the 
color  display.  The  unprocessed  waveforms  can  also  be  sent  to  the  S/200  for 
storage  on  the  disk  memory.  The  display  unit  can  be  programmed  in  a  general 
way  to  display  the  data  in  many  different  formats.  For  image  data,  a  B-scan 
display  is  the  most  useful,  although  C-scans  can  be  displayed  if  desired.  The 
scattering  data  is  not  displayed  but  only  recorded  for  use  with  the  inversion 
algorithms.  Detailed  schematic  drawings  of  the  circuits  in  the  hybrid  pulser/ 
receiver,  the  multiplexer,  the  signal  conditioner  and  the  blocks  contained  in 
the  dotted  outline  can  be  found  in  Appendix  B. 

4.2.1  Transducers 

The  first  major  subassembly  is  the  transducer  array.  Originally  this 
was  to  be  a  6  inch  transducer  containing  240  elements  which  could  be  used  for 
contour  scanning.  An  array  of  this  sort  must  possess  a  sufficiently  wide 
bandwidth  and  short  transient  response  that  it  can  be  used  to  detect  flaws 
below  metal  surfaces.  This  condition  requires  the  "ring  down"  from  the  front 
surface  echo  to  decay  by  about  20  dB  before  the  flaw  signal  occurs.  Further¬ 
more  the  cross-talk  between  array  elements  must  be  about  25  dB  below  the 
direct  signal. 

A  vendor  was  contacted  who  claimed  that  a  custom  transducer  satis¬ 
fying  these  requirements  could  be  fabricated.  A  sub-contract  was  issued  to 
build  the  array  and  a  unit  was  manufactured.  Preliminary  testing  of  the  unit 
revealed  that  the  transient  response  was  too  long  by  an  order  of  magnitude. 
Also  the  peak  amplitude  of  the  individual  element  respones  varied  by  more  than 
±  10  dB  about  the  average  level.  The  vendor  made  several  attempts  to  improve 
this  performance  without  success.  Conversations  with  other  experts  Indicated 
that  the  current  state  of  array  technology  is  not  sufficiently  advanced  to 
manufacture  an  array  this  large  meeting  the  required  specifications. 
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To  obtain  the  required  transducer  characteristics  the  system  was 
reconfigured  to  use  two  32  element  arrays.  These  arrays  have  a  nominal  center 
frequency  of  2.25  MHz.  While  a  higher  frequency  is  desirable  for  NDE  applica¬ 
tions,  no  commercial  sources  of  such  arrays  could  be  located.  The  peak  ampli¬ 
tude  and  the  duration  of  the  transient  response  of  each  of  the  32  elements 
were  carefully  measured.  The  peak  amplitudes  were  less  uniform  than  desired 
with  the  amplitudes  of  95%  of  the  elements  being  within  ±37%  of  the  mean 
value.  This  variation  is  less  than  that  encountered  with  many  commercially 
available  transducers  and  was  accepted.  The  data  on  the  duration  of  the 
transient  response  is  summarized  in  Fig.  4.5.  Since  the  exact  shape  of  the 
transient  response  varies  among  transducer  elements,  we  have  measured  the 
decay  times  relative  to  the  time  of  the  peak  of  the  transient  response.  In 
this  way  the  curves  shown  in  Fig.  4.5  are  obtained.  The  thin  line  represents 
the  average  decay  time  of  the  32  elements.  The  thick  lines  form  an  envelope 
that  contains  the  decay  times  of  95%  of  the  elements.  A  third  curve  shows  the 
theoretical  decay  time  for  a  similar  transducer  element.17  This  decay  time  is 
considerably  shorter  than  that  of  the  array  elements  and  is  approximated  by 
the  performance  of  a  carefully  fabricated  single  element  transducer.18  The 
array  performance  is  adequate  for  our  needs  and  far  exceeds  the  performance  of 
the  240  element  array  that  was  tested  early  in  the  program.  One  of  the 
transducer  arrays  is  shown  in  Fig.  4.6. 

4.2.2  Pul ser/Recei vers 

The  next  subassembly  after  the  transducer  in  Fig.  4.4  is  the  pulser/ 
receiver  package.  Several  design  limitations  dictated  that  the  pulser/recei ver 
be  on  the  array  side  of  the  multiplexer.  First,  a  high  voltage  transmit  pulse 
is  required  for  the  transducer  excitation.  Second,  the  receiving  circuits 
require  a  voltage  limiting  transmit/receive  switch.  Third,  the  received  signals 
have  to  be  provided  with  a  local  low  noise  amplifier  to  overcome  any  low  level 
multiplexer  noise.  This  means  that  it  will  be  necessary  to  have  a  pulser/ 
receiver  for  each  transducer  element  or  a  total  of  64  for  this  system.  Placing 
the  pulser/recei ver  at  a  remote  location  would  be  prohibitive  because  of  the 
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bulk  of  the  cables  that  would  be  required  for  the  Interconnection  and  the 
possibility  of  severe  signal  degradation  due  to  cable  capacitance.  These 
considerations  led  to  the  decision  to  package  the  pulser/recelver  as  compactly 
as  possible  and  to  build  the  unit  as  a  hybrid  circuit.  This  imposed  several 
constraints  on  the  circuit  design  relating  to  availability  of  components  In  chip 
form,  to  the  power  dissipation  of  the  units,  and  to  the  size  of  the  completed 
circuit. 

The  pulser  circuit  is  shown  In  Fig.  4.7a.  The  switching  element  Is  an 
SCR.  A  number  of  alternatives  were  Investigated  Including  a  transistor  operated 
in  the  avalanche  mode  and  one  operated  as  a  hi t»h  voltage  switching  transistor. 
The  switching  transistor  was  found  to  be  too  slow  for  this  application.  The 
avalanche  switching  speed  was  more  than  adequate  but  it  was  difficult  to  obtain 
repeatable  results  with  a  selection  of  transistors.  Furthermore,  the  units 
would  only  operate  reliably  over  a  small  range  of  applied  dc  voltages.  The  SCR 
has  an  adequate  switching  speed  and  operates  reliably  over  a  wide  range  of  dc 
voltages.  This  voltage  range  is  sufficient  to  vary  the  output  amplitude  of  the 
transmitted  pulse  by  more  than  30  dB.  Transistor  Qj  Is  used  for  shaping  the 
current  pulse  produced  by  the  TTL  Input  pulse  so  that  It  will  reliably  drive  the 
SCR  (labeled  01).  Transistor  Q5  is  used  to  speed  up  the  charging  rate  of  C2  so 
that  the  pulser  can  be  operated  at  repetition  rates  of  1  KHz.  The  pulser  has 
been  used  to  excite  a  2.25  MHz  commercial  transducer.  The  echo  returned  from  a 
block  of  aluminum  has  been  compared  with  that  obtained  from  a  commercial  pulser 
(Panametrlcs  Model  5052)  using  the  same  transducer.  The  results  are  compared  In 
Fig.  4.8.  It  Is  readily  seen  that  the  transient  response  of  the  transducer  Is 
essentially  the  same  In  both  cases.  Thus  we  can  conclude  that  the  switching 
speed  of  the  SCR  Is  sufficiently  fast  that  the  spectral  content  of  the  pulse 
produced  will  fully  excite  the  bandpass  of  the  transducer.  The  receiver  circuit 
Is  shown  In  Fig.  4.7b.  The  signal  path  Is  composed  of  transistors  Qj  and  Q3, 
while  Q4  Is  used  to  turn  It  on  and  off.  The  gain  of  the  circuit  Is  principally 
determined  by  the  ratio  of  Rx  to  Ry  and  thus  can  be  made  relatively  Independent 
of  variations  In  the  active  elements  and  the  temperature.  The  gain  has  been  set 
at  26  dB.  The  circuit  Is  ac  coupled  to  avoid  problems  with  dc  drift  and  has  a 
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low  frequency  3  dB  point  of  0.7  MHz  and  a  high  frequency  3  dB  point  of  about  30 
MHz.  The  complete  frequency  response  of  a  typical  receiver  Is  shown  In 
Fig.  4.9.  The  low  frequency  3  dB  point  is  a  compromise  between  the  need  to 
provide  a  sufficiently  broad  bandwidth  to  prevent  phase  distortion  of  the 
signals  and  the  need  to  provide  a  rapid  recovery  after  being  driven  into 
saturation  by  a  high  amplitude  signal.  The  recovery  time  is  observed  to  be 
1  nsec  and  no  signal  distortion  can  be  detected  when  commercial  transducers  are 
used.  The  power  dissipation  of  the  receiver  is  only  70  milliwatts.  The  pro¬ 
vision  to  turn  the  receiver  on  with  a  TTL  pulse  when  it  is  among  the  16 
addressed  receivers  reduces  the  cross-talk  from  receivers  that  are  not 
addressed. 

Two  complete  hybrid  pul ser/recei vers  were  mounted  on  a  single 
substrate.  This  choice  minimized  the  board  complexity  and  kept  the  hybrid 
package  from  becoming  unwieldy.  The  package  shown  In  Fig.  4.10  is  a  3.4  cm  x  2 
cm  x  0.5  cm  hermetically  sealed  can  with  24  pins.  Two  voltages  are  needed  to 
operate  the  unit:  five  volts  for  both  the  trigger  input  stage  of  the  pulser  and 
the  receiver  stage  and  a  transducer  excitation  voltage  that  can  be  varied  from 
about  five  volts  to  four  hundred  volts. 


4.2.3  Multiplexer 

The  multiplexer  section  consists  of  three,  64  to  16  multiplexers.  The 
total  of  64  transducer  elements  are  grouped  into  segments  that  each  contain  four 
elements.  The  four  elements  in  any  one  segment  are  switched  as  a  unit.  Thus  it 
Is  possible  to  digitally  select  13  distinct  groups  of  16  contiguous  pulsers 
(four  segments)  from  among  the  64  available.  Similarly  and  independently,  16 
contiguous  receivers  can  be  selected.  There  are  two  digital  multiplexers:  one 
for  supplying  triggers  to  the  selected  pulsers  and  one  for  enabling  the  selected 
receivers.  There  Is  one  analog  multiplexer  that  Is  used  to  convey  the  output 
signals  of  the  selected  receivers  to  the  16  receive  channels. 

The  coding  scheme  for  the  multiplexers  can  best  be  understood  using  the 
diagram  In  Fig.  4.11.  Each  of  the  numbered  rectangles  corresponds  to  four 
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elements  (l.e.,  pulsers  or  receivers).  In  the  actual  hardware  each  of  the 
switches  has  only  four  positions  rather  than  the  arbitrarily  large  number  shown 
In  Fig.  4.11.  The  four  switch  positions  provide  13  distinct  groups  of  elements 
that  can  be  selected.  The  switching  is  accomplished  by  initially  sending  a 
master  code  on  the  line  labeled  D  that  selects  the  same  corresponding  position 
for  each  of  the  four  switches.  Note  that  this  automatically  selects  four 
contiguous  segments  of  elements  containing  a  total  of  16  elements.  To  shift  the 
selected  group  of  four  segments  by  an  Integral  multiple  of  four  segments  (16 
elements),  the  master  code  Is  simply  advanced  or  retarded  by  n,  where  n  is  the 
appropriate  Integer.  For  shifts  of  less  than  four  segments,  the  auxiliary  codes 
are  used.  A  shift  of  one  segment  is  provided  by  changing  the  “A"  code  by  one,  a 
shift  of  two  segments  results  from  changing  the  "A"  and  "B"  codes  by  one,  etc. 

In  this  way  It  Is  possible  to  select  any  group  of  16  contiguous  elements  in 
groups  of  four  segments. 

The  multiplexer  subassembly  Is  contained  on  three  printed  circuit 
boards  which  are  partitioned  and  Interconnected  as  shown  in  Fig.  4.12.  The 
circuitry  on  board  1  serves  to  buffer  and  decode  the  five  bit  transmit  and  five 
bit  receive  codes  that  are  sent  from  the  timing  and  control  circuitry.  Also 
located  on  this  board  is  the  circuitry  for  generating  the  eight  bit  code  that  is 
required  to  control  the  pulser  and  receiver  selection.  The  remaining  two  boards 
are  Identical.  Each  board  contains  half  of  the  multiplexers  and  half  of  the 
pulser/recelver  packages.  The  64  transducer  elements  are  arranged  in  a  row  and 
for  the  purposes  of  switching  are  divided  into  segments  of  four  elements  each. 
The  16  segments  are  numbered  from  left  to  right.  The  multiplexers  and  pulser/ 
receivers  corresponding  to  the  even  numbered  elements  are  located  on  board  2  and 
those  corresponding  to  the  odd  numbered  elements  are  on  board  3.  Fig.  4.13  is  a 
photograph  of  one  of  these  boards  and  Fig.  4.14  is  a  block  diagram  of  the 
multiplexer. 

As  shown  In  Fig.  4.14,  the  transmit  elements  are  selected  with  the  XMT 
code  which  Is  decoded  and  used  to  switch  the  XMT  MUX  to  the  desired  positions. 
The  triggers  for  the  hybrid  pulsers  are  then  transmitted  through  the  multiplexer 
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Fig.  4.12  Partitioning  and  Interconnection  of  multiplexer. 


Fig.  4.13  Photograph  of  multiplexer  board. 


to  activate  the  selected  transducer  elements.  The  receiver  elements  are 
selected  In  an  analogous  way  with  the  RCV  code.  In  this  case  the  code  switches 
both  the  digital  and  the  analog  multiplexers.  The  digital  unit  Is  used  to 
enable  16  of  the  receiver  circuits  by  switching  the  receiver's  +5  volt  power 
path  Internal  to  each  hybrid.  The  analog  unit  connects  the  outputs  of  the  16 
enabled  receivers  to  the  16  line  drivers.  These  line  drivers  connect  the 
multiplexer  assembly  to  the  signal  conditioning  circuitry  through  75  ohm  coaxial 
cables.  The  specific  channels  selected  by  each  of  the  group  codes  are  listed  in 
Table  4.1. 


Table  4.1 

XMT/RCV  Segment  Select  Codes 


Array  Segment  Selection 
(4  elements/segment) 

Binary 

Code 

Octal 

Code 

MSB  LSB  C  B  A 

1  thru  4 

0 

0 

0 

0 

0 

00 

2  thru 

5 

0 

0 

0 

0 

1 

01 

3  thru 

6 

0 

0 

0 

1 

1 

03 

4  thru 

7 

0 

0 

1 

1 

1 

07 

5  thru 

8 

0 

1 

0 

0 

0 

10 

6  thru 

9 

0 

1 

0 

0 

1 

11 

7  thru 

10 

0 

1 

0 

1 

1 

13 

8  thru 

11 

0 

1 

1 

1 

1 

17 

9  thru 

12 

1 

0 

0 

0 

0 

20 

10  thru 

13 

1 

0 

0 

0 

1 

21 

11  thru 

14 

1 

0 

0 

1 

1 

23 

12  thru 

IS 

1 

0 

1 

1 

1 

27 

13  thru 

16 

1 

1 

0 

0 

0 

30 
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It  is  essential  that  the  analog  multiplexer  used  for  the  received 
signals  have  crosstalk  and  feedthrough  levels  that  are  negligible  compared  with 
the  direct  signals.  The  design  of  the  multiplexer  function  is  such  that 
there  are  16  individual  multiplexer  Integrated  circuit  packages  in  the  analog 
section.  Each  IC  connects  to  only  one  active  transducer  element  at  a  time.  The 
IC's  are  sufficiently  well  separated  that  the  crosstalk  and  feedthrough  between 
them  is  of  no  concern.  It  is  the  Interactions  between  the  four  inputs  of  the 
individual  IC  that  might  cause  a  problem.  Using  a  bread-board  model,  a  known 
signal  with  a  particular  source  resistance  was  applied  to  all  the  input  pins 
except  for  the  one  pin  that  was  connected  to  the  output  pin.  This  input  pin  was 
connected  to  ground  through  a  39  ohm  resistor.  The  signal  level  at  the  output 
was  then  compared  to  the  input  signal  level  applied  to  the  input  pins.  The 
output  was  found  to  be  35  dB  below  the  input  signal  level.  Since  an  8  bit  A/D 
converter  has  a  dynamic  range  of  48  dB,  this  crosstalk  level  is  not  negligible. 
However,  In  the  actual  circuit,  the  receivers  for  the  other  input  pins  are  not 
enabled,  providing  at  least  another  26  dB  of  Isolation  (corresponding  to  the 
receiver  amplifier  gain).  Therefore  the  overall  isolation  of  the  analog  output 
signal  from  signals  at  inactive  multiplexer  inputs  is  in  excess  of  60  dB.  This 
degree  of  Isolation  Is  acceptable  and  provides  a  comfortable  margin  over  the 
minimum  signal  that  would  affect  the  output  of  the  A/D  converter.  The  actual 
performance  of  the  system  will  not  be  this  good  since  the  weak  link  in  the 
effort  to  suppress  crosstalk  is  the  transducer  array.  The  crosstalk  of  the 
particular  array  used  in  this  system  was  not  measured,  but  the  level  is 
reportedly  less  than  30  dB  below  the  direct  signals. 

The  printed  circuit  boards  comprising  the  multiplexer  assembly  are 
mounted  In  a  waterproof  box  with  dimensions  of  14  Inches  x  15  Inches  x  16 
Inches.  This  box  is  attached  to  the  two  transducer  arrays  through  64  coaxial 
cables.  The  lengths  of  the  cables  from  the  transducers  to  the  connectors  on  the 
p.c.  board  are  between  17  and  25  Inches  depending  on  the  physical  location  of 
the  p.c.  board  connector.  The  multiplexer  box  and  the  transducer  arrays  are 
held  in  place  by  a  special  frame  that  will  be  described  in  Section  4.2.9. 
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1.2.4  Signal  Conditioning  Circuitry 

The  purpose  of  the  signal  conditioning  circuitry  Is  to  prepare  the 
signals  for  the  A/0  converter.  The  circuitry  Is  composed  of  the  main  amplifier, 
a  digitally  controlled  attenuator,  a  clipping  circuit,  a  filter,  a  dc  level 
shifter  and  a  low  output  Impedance  buffer  amplifier.  A  schematic  of  the  circuit 
is  shown  in  Fig.  4.15  a  and  b.  (Drawing  B-2  in  Appendix  B  provides  component 
values  for  this  circuitry.)  The  main  amplifier  has  a  gain  of  26  dB  and  a  band¬ 
width  of  10  mz  with  a  maximum  peak  to  peak  output  voltage  of  23  volts.  It  is 
based  on  a  commercially  available  operational  amplifier  with  external 
compensation  to  realize  the  wide  bandwidth  and  high  output  voltage. 

Following  the  amplifier  Is  a  digitally  controlled  attenuator.  The 
design  of  the  attenuator  is  based  on  a  tapped  constant  impedance  ladder  net¬ 
work.  There  are  15  taps  with  a  3  dB  change  in  attenuation  at  each  tap  to 
provide  a  total  attenuation  range  of  45  dB.  A  16:1  analog  multiplexer  Is  used 
to  switch  between  the  taps.  In  order  to  obtain  the  requisite  bandwidth  of  10 
MHz  and  reduce  feedthrough  problems.  It  was  necessary  to  build  the  multiplexer 
as  a  two  stage  unit  with  a  buffer  amplifier  between  the  two  stages.  It  also 
proved  to  be  convenient  to  put  the  signal  clipper  composed  of  diodes  01  through 
D6  at  the  Input  to  the  buffer  amplifier.  The  buffer  amplifier  consists  of 
transistors,  Qj,  Q2,  and  Q3.  It  Is  a  constant  current  biased  emitter  follower 
where  the  current  source  Q3  Is  switched  from  one  emitter  follower  to  another 
depending  on  which  attenuator  tap  Is  being  addressed. 

The  signals  are  then  passed  through  the  filter.  The  filter  Is  designed 
as  a  6  pole,  low  pass,  Butterworth  filter  with  a  3  dB  frequency  of  6.6  MHz. 

This  type  of  filter  was  chosen  because  of  its  sharp  cutoff  characteristic.  It 
serves  as  an  anti-aliasing  filter  and  removes  digital  noise  originating  from  the 
18  MHz  clocking  signals.  The  output  of  the  filter  Is  ac  coupled  to  the  dc  level 
shifting  circuit  (852*  #53*  R54)  and  then  coupled  Into  the  buffer  amplifier. 

The  buffer  amplifier  has  a  low  output  impedance  that  effectively  suppresses 
transients  caused  by  the  large  current  pulses  associated  with  the  strobing  of 
the  A/0  converter. 
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Fig.  4.15  (b)  Schematic  of  signal  conditioning  circuitry:  filter,  dc  level 

shifter,  and  buffer  amplifier. 


4.2.5  Array  Digitizing  Electronics 

Referring  to  the  block  diagram  In  Fig.  4.4,  the  set  of  blocks  contained 
within  the  dotted  lines  comprise  the  Array  Digitizer  Electronics  or  ADE.  These 
functional  blocks  are  interconnected  through  a  variety  of  lines  carrying  arm, 
trigger,  clocking  and  data  signals  from  the  Interface  and  the  timing  and  control 
circuitry.  The  purposes  of  the  ADE  can  be  summarized  as: 

1.  Accept  and  decode  commands  and  data  from  the  host  computer. 

2.  Switch  the  multiplexer  to  the  appropriate  sets  of  pulsers  and 
receivers. 

3.  Select  the  gain  of  each  receiver  channel. 

4.  Send  trigger  signals  to  the  pulsers  with  the  appropriate  time 
delays. 

5.  Digitize  the  signals  received  after  a  specified  delay. 

6.  Store  the  digitized  signals  in  a  buffer  memory. 

7.  Transmit  the  contents  of  the  buffer  memory  to  the  array  processor. 

In  order  to  accomplish  these  tasks  there  are  a  total  of  eight  command 
codes  that  are  sent  by  the  host  computer  and  a  starting  memory  address  code  that 
is  sent  by  the  array  processor.  The  host  computer  codes  are  listed  In  Table 
4.2.  When  the  codes  are  received  from  the  host  computer  they  are  Immediately 
identified  and  the  data  associated  with  the  command  Is  stored  In  the  appropriate 
register  of  the  control  memory,  schematically  Illustrated  In  Fig.  4.16.  Two  of 
the  commands,  the  host  trigger  and  the  host  clear,  are  unique  in  that  no  data  Is 
associated  with  them  and  they  are  executed  Immediately  after  being  Identified. 

The  function  of  much  of  this  circuitry  can  be  explained  by  tracing  the 
operations  that  take  place  when  an  array  pulse  and  receive  sequence  is  exe¬ 
cuted.  The  first  command  to  be  sent  by  the  host  computer  Is  10  CIR.  This 
Initializes  the  circuits  to  a  known  state.  The  exact  sequence  of  the  commands 
following  an  10  CLR  Is  not  Important  since  their  function  Is  to  store  data  In 
the  six  registers  of  the  control  memory  and  this  can  be  done  In  any  order.  A 
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Table  4.2 

Host  Command  Codes 


Functions 

Acronym 

Command 
Code  (Octal) 

Data  Field 

Bits 

Trigger  Delay  Coarse 

TDC  01 

00 

8 

TDC  92 

01 

TDC  93 

02 

TOC  04 

03 

TDC  05 

04 

TDC  06 

05 

TDC  07 

06 

TDC  08 

07 

TDC  09 

10 

TDC  10 

11 

TDC  11 

12 

TDC  12 

13 

TDC  13 

14 
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Fig.  4.16  Schematic  of  control  memory 


typical  sequence  Is  to  first  send  the  command  for  selecting  the  set  of  16 
transducer  array  elements  to  be  used  for  transmitting.  The  data  field  speci¬ 
fying  the  five  bit  code  for  the  multiplexer  Is  stored  In  the  XMTG  register  of 
the  control  memory.  It  Is  then  sent  to  the  decoder  board  of  the  multiplexer 
where  the  set  of  transducer  array  elements  are  selected  according  to  the  scheme 
Illustrated  In  Fig.  4.11.  A  similar  procedure  Is  followed  for  the  selection  of 
the  16  receiving  transducer  elements  using  the  RCVG  command.  Next  the  gain  of 
each  of  the  16  receiver  channels  Is  selected.  The  data  fields  specifying  the 
four  bit  codes  for  the  digital  attenuator  multiplexers  are  stored  in  the  16  RGV 
registers  of  the  control  memory.  These  codes  are  then  sent  to  the  attenuator 
multiplexers. 

The  remaining  three  data  commands  all  specify  delays.  There  are  both 
coarse  and  fine  delays  to  be  specified  for  the  transmit  triggers.  Also  there  is 
a  delay  corresponding  to  the  ultrasonic  water  path  transit  time  that  Is 
specified  before  starting  the  A/D  conversion  process.  The  procedure  for  pro¬ 
cessing  the  transmit  delay  commands  is  illustrated  in  Fig.  4.17.  There  are  16 
transmit  elements  with  delays  for  each  one  specified  separately.  The  fine  delay 
command,  TDF,  results  In  the  four  bit  data  field  being  stored  in  the  TDF 
register  and  then  being  used  to  select  the  desired  fine  delay.  The  fine  delay 
Is  obtained  with  a  commercially  available  digitally  controlled  delay  line.  This 
unit  Is  capable  of  selecting  16  different  delays  In  steps  of  4  nsec  using  a  four 
bit  Input  code.  The  fine  delays  are  useful  for  reducing  the  quantization  error 
associated  with  the  specified  time  delays  for  a  given  steering  direction  and  for 
trimming  each  of  the  channels  to  correct  for  variations  In  their  Inherent 
delays.  The  coarse  delay  command,  TDC,  results  In  the  eight  bit  data  field 
being  stored  In  the  TDC  register  which  controls  the  delay  generator.  This  delay 
Is  specified  In  integer  multiples  of  a  clock  Interval  or  «  55  ns  steps.  After 
the  specified  number  of  clock  pulses  have  been  received,  a  trigger  Is  sent 
through  the  fine  delay  circuit  to  the  transmitter  multiplexer  assembly  as  shown 
In  Fig.  4.17.  The  other  delay  command,  known  as  the  CONV  DLY,  Is  stored  In  the 
CONV  DLY  register.  This  data  Is  held  In  the  register  until  the  main  trigger  Is 
transmitted.  This  delay  Is  the  same  for  all  16  receiving  channels.  After  this 
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delay  a  gate  Is  generated  that  starts  the  A/0  converters  and  starts  storing  the 
waveform  samples  in  the  buffer  memory  as  shown  in  Fig.  4.18. 

After  all  six  control  memory  registers  have  been  filled  the  host 
command  10  TRG  is  sent  to  the  decoder  circuit  and  is  executed  immediately.  This 
results  in  the  entire  pulse  and  receive  sequence  being  executed  according  to  the 
data  In  the  control  memory.  Once  the  digitizing  process  starts.  It  continues 
automatically  until  1024  samples  have  been  stored.  The  address  generator  shown 
In  Fig.  4.18  keeps  a  count  of  the  number  of  samples  stored.  When  the  memory  is 
full  a  write  finished  (WR  FIN)  signal  Is  generated  that  stops  the  operation  of 
the  A/D  converter. 

The  Digitizer  section  Is  shown  In  Fig.  4.19.  It  consists  of  the  8  bit 
A/D  converter  (the  TRW  Model  TDC  1007J  operating  at  a  sampling  frequency  of  18 
MHz),  a  synchronizing  latch,  and  a  1  kilobit  *  8  bit  RAM  buffer  memory.  The 
clocking  speeds  being  used  and  the  close  tolerances  on  the  clocking  signals 
dictate  that  all  components  of  the  digitizer  be  operated  as  a  closely  knit  unit 
sharing  the  same  p.c.  board. 

The  next  step  In  the  operation  is  to  transfer  the  contents  of  the  16 
buffer  memories  to  the  array  processor  for  the  beam  forming  operation.  This 
transfer  is  Initiated  by  the  array  processer  which  places  an  address  on  the 
Interface  line  as  shown  In  Fig.  4.18.  This  address  plus  "handshake"  between  the 
ADE  and  the  array  processor  will  start  a  "memory  read"  operation  that  will 
transfer  the  contents  of  the  memory  to  the  array  processor.  The  sections  of  the 
memory  that  are  to  be  transferred  are  specified  by  the  array  processor. 

4.2.6  S/200  Software  Utilized  for  Control  of  Array  and  AP400 

The  interactive  S/200  resident  software  package,  ARRAY,  allows  a  user 
to  control  and  coordinate  the  activities  of  both  the  array  digitizer  electronics 
(ADE)  and  the  Analogic  AP400  array  processor  (AP).  The  package  supplies  the 
user  with  the  ability  to  easily  command  the  acquisition  of  a  full  set  of 
waveforms  with  a  minimum  of  user  Input,  as  well  as  to  single-step,  l.e.,  command 
by  single  Instruction,  the  Initiation  of  the  AOE  and  acquisition  of  data  from 
the  ADE  via  the  AP400.  The  existing  version  of  the  software  does  not  take 
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advantage  of  all  of  the  capabilities  of  the  ADE  hardware.  It  provides  access  to 
a  subset  of  the  hardware  that  will  permit  the  generation  of  sector  scans. 

ARRAY  consists  of  both  FORTRAN  5  and  assembly  language  modules.  It  was 
written  for  a  Data  General  ECLIPSE  S/200,  and  runs  under  RDOS.  The  modules  of 
ARRAY  are  listed  with  their  function  in  Table  4.3. 


( 


In  addition,  ARRAY  calls  many  routines  supplied  in  the  AP400  library 
for  communication  with  the  AP400;  further  Information  on  the  AP400  is  Included 
In  Section  4.2.7. 

Upon  execution,  ARRAY  Initializes  variables  as  well  as  both  the  AP400 
and  the  AOE.  The  ADE  Is  Initialized  at  this  point  from  the  perspective  that  It 
has  been  defined  to  the  software.  Upon  execution  of  ARRAY,  the  user  Is 
presented  with  several  options,  ^ach  specified  by  a  2-letter  mnemonic  as  shown 
In  Table  4.4.  Selection  of  options  "AD",  "AP",  "PL",  and  "ST"  pass  control  to 
modules  ADEC,  AP400,  DATAPLOT,  and  WISP  respectively;  necessary  user  Input  to 
each  of  these  options  Is  detailed  below. 


Table  4.3 

Primary  ARRAY  Modules 


ARRAY.FR 

System  Initialization,  operator  Interface,  procedure 
controller,  and  Initiation  of  other  modules. 

adec.fr 

Operator  Interface  for  direct  communication  to 
initialize  the  ADE. 

ADE. SR 

Device  handler,  ADE;  10  is  programmed. 

AP400.FR 

Operator  Interface  for  acquisition  of  a  single  set  of 
data  from  the  ADE. 

DATAPLOT.FR 

Plots  most  recently  acquired  waveform  on  TEKTRONIX 
screen. 

wisp.fr 

Writes  ISP  format  data  files. 

! 

1 
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Table  4.4 
Commands  to  Array 


Mnemonic 

Action 

AD 

Communicate  with  ADE  in  single  command  mode. 

AP 

Communicate  with  AP  in  single  command  mode. 

HELP 

List  available  options. 

PL 

Plot  on  terminal  the  latest  waveform  acquired. 

PR 

Initiate  the  procedure  of  acquiring  a  set  of 
data. 

ST 

Store  latest  waveform  acquired  on  disk. 

XX, BYE 

Exit  ARRAY. 

For  the  standard  acquisition  of  data  from  a  range  of  angles  the  command 
"PR"  is  specified.  Upon  entry  to  this  option  of  ARRAY,  the  user  is  queried  for 
sector  scan  definition  which  includes  the  total  number  of  waveforms,  the  number 
of  degrees  in  the  sector,  and  the  starting  angle.  In  addition,  the  output  disk 
file  name  for  the  storage  of  the  waveforms,  the  initialization  parameters  for 
the  ADE,  and  the  number  of  channels  to  sum  (1-16)  must  be  specified.  The 
transfer  of  these  Input  parameters  to  the  ADE,  the  subsequent  acquisition  of 
data,  and  the  activities  of  the  array  processor  are  transparent  to  the  user. 
ARRAY  notifies  the  user  of  the  completion  of  signal  processing  for  each  of  the 
waveforms  corresponding  to  alternate  lines  in  the  sector  scan  and  also  of  the 
completion  of  the  entire  sector.  A  sample  screen  display  for  the  operation  of 
ARRAY  is  shown  in  Fig.  4.20. 

The  angles  corresponding  to  the  Individual  lines  of  the  sector  scan  are 
used  to  calculate  the  trigger  delays  for  each  channel  which  are  then  stored  in 
the  control  memory  of  the  ADE.  There  are  four  other  sets  of  parameters  that 
need  to  be  stored  In  the  control  memory  before  the  data  acquisition  can  begin. 
These  consist  of  the  specification  of  the  group  of  transmitting  transducer 
elements,  the  group  of  receiving  transducer  elements,  the  attenuator  setting  for 
each  of  the  receiver  channels,  and  the  delay  before  the  A/D  converters  start 
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ADC  INITIATION 


ACTUAL  SAMPLE 


RANGE  OF 
VALID  RESPONSES 


PR 

■  WFMSIN  FAN. 


r  DEG  IN  FAN 


BEGIN  ANGLE  I  -45  ]  :  1-51 

FILE  NAME  -  fPOlEXAM.CCl 
'INIT  ADE  W/FILE  (D  OR  KB(2l  1  1  \ 

CHANNEL  CHOICE 
CODE  CHANNELS 
0  ■  1-16 

1  ■  17-32 

2  •  33-48 

3  •  49-64  _ 

ENTER  CODE  FOR  XMIT  AND  RCV  ]  2.  2| 

FOR  RCVR  VLTG  GAIN. 

ENTER  LOW  CH=  OF  PAIR  (0.2  .14).  AND  2  STEP  VALS  10-15) 
EX  0.13,0  YIELDS  13  STEPS  FOR  CHO  AND  0  FOR  CHI 
ENTER  ICH.V1.V2  [0.1,1 
ENTER  ICH.V1.V2  2.1.1 
ENTER  ICH.V1.V2  3  4,1,1 
ENTER  ICH.V1.V2  6,1.1 
ENTER  ICH.V1.V2  8,1,1 
ENTER  ICH.V1.V2  10,1,1 
ENTER  ICH.V1.V2  ij.yi 
ENTER  ICH.V1.V2  14,1,1 
ENTER  ICH.V1.V2  -1.0.6 


•  WFMS  >  1 

*  DEGREES:  >  1 

-45  <  BEGIN  ANGLE  <  *  45 


1- YES.  ANY  OTHER  RESPONSES  ->  NO 


r:  0.  1.2.3  1 
0. 1.2.3  J 


INDEPENDENT 


A/D  DELAY  (USECS)  :|60| 
CHANS  TO  SUM  (16)  .  QF] 


DONE 

W/ANG 
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WANG 
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DONE 

W/ANG 

X 

8 
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3  68421 
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6.00000 
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0  <  DELAY  <  4095 
(55  NS  STEP) 

1  <t  r  CHANS  C  16 
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Fig.  4.20  Sample  screen  display  for  ARRAY. 
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operating.  The  specification  of  these  parameters  Is  termed  the  ADE  Initia¬ 
tion.  The  values  of  these  variables  may  be  left  unchanged  during  subsequent 
passes  through  “PR". 

The  existing  software  only  permits  the  selection  of  four  groups  of 
transmitter  and  receiver  transducer  elements  rather  than  the  sixteen  that  the 
hardware  is  capable  of  selecting.  These  groups  represent  the  transducer 
elements  1-16,  17-32,  33-48,  and  49-64  which  are  specified  by  the  code  numbers 
0,  1,  2,  and  3  respectively.  For  a  sector  scan  the  transmitter  and  receiver 
groups  are  identical. 

The  receiver  voltage  gain  must  initially  be  specified  for  each  of  the 
sixteen  receiver  channels;  In  subsequent  passes  through  PR,  the  user  may  update 
as  few  as  one  pair.  The  gain  Is  specified  for  two  channels  (one  pair)  on  each 
entry  in  terms  of  attenuator  steps  with  each  step  corresponding  to  an  attenua¬ 
tion  of  3  dB.  For  example  the  entry,  "0,  1»  3“  specifies  the  attenuation  for 
channels  0  and  1.  The  first  digit  specifies  the  even  number  of  the  pair,  0,  In 
this  case.  The  second  digit  specifies  the  number  of  attenuator  steps  for  the 
first  or  even  channel  of  the  pair.  In  this  case  1  Is  specified  and  results  in 
an  attenuation  of  3  dB  for  channel  0.  The  third  digit  specifies  the  number  of 
attenuator  steps  for  the  second  or  odd  channel  of  the  pair.  In  this  case  the  3 
results  In  an  attenuation  of  9  dB  for  channel  1.  Similarly  the  specification, 
"14,  0,  0“  sets  the  attenuator  In  both  channels  14  and  15  to  0  dB.  The  user 
should  only  enter  even  numbers  as  his  first  channel  specification;  an  entry  of 
"13,  0,  0"  would  yield  unpredictable  results. 

The  final  parameter  needed  to  Initialize  the  ADE  Is  the  A/D  delay  which 
Is  specified  In  microseconds,  and  may  vary  from  0  to  225  nsec.  This  Is  the 
conversion  delay  described  In  Section  4.2.5  and  permits  an  Interval  of  time 
corresponding  to  the  water  path  delay  to  elapse  before  digitization  of  the 
waveform  begins. 

After  the  user  specifies  the  number  of  receiver  channels  to  process,  an 
10  TRG  Is  automatically  sent  to  the  ADE  and  the  data  acquisition  begins.  After 
this  step  Is  completed,  the  data  Is  automatically  transferred  to  the  array 
processor  for  the  "beam  forming"  process.  For  sector  scans  the  angles  specified 
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for  the  receiver  are  the  same  as  those  for  the  transmitter;  no  additional  user 
Input  Is  required.  The  details  of  this  beam  forming  operation  are  provided  In 
Section  4.2.7.  After  the  beam  forming  operation  is  completed,  the  resulting 
waveform  Is  passed  back  to  the  host,  and  written  to  disk.  The  user  Is  notified 
as  this  occurs,  on  alternate  waveforms. 

Both  of  the  commands  "AD"  and  "AP"  provide  a  diagnostic  capability  to 
the  user.  The  command  "AD"  allows  the  user  to  send  octal  Integer  format 
commands  directly  to  the  ADE,  while  verifying  their  operation  by  use  of  a  scope 
or  by  acquiring  data  by  use  of  "AP".  The  commands  are  5  digits  In  length  and 
are  of  the  format  IIDDD  where  "I"  represents  a  destination  code  and  "D" 
represents  data. 

There  are  approximately  35  commands  available  to  control  the  ADE.  The 
correct  ordering  requires  one  10  CLR  to  precede  any  other  commands  after  ADE 
power-up  or  reset.  The  normal  sequence  Is  then  a  set  of  codes  followed  by  a 
trigger  that  actually  activates  the  requests.  Table  4.2  in  Section  4.2.5 
details  the  valid  commands  which  may  be  used.  The  three  octal  digits 
representing  the  data  field  bits  are  appended  to  each  of  the  command  codes  in 
Table  4.2  according  to  guidelines  in  Table  4.5 

If  desired  the  user  may  acquire  data  from  a  user-specified  location  In 
the  ADE  memory.  The  value  of  this  option  Is  Its  ability  to  acquire  data  from  a 
particular  ADE  receiver  channel  to  verify  Its  operation.  It  Is  used  In 
conjunction  with  the  "AD"  option.  The  Input  arguments  specify  the  memory, 
offset  address,  and  number  of  points  to  acquire.  The  memory  Is  specified  by  a 
number  from  0  to  15,  the  offset  address  by  a  number  from  0  to  1023,  and  the 
number  of  points  by  a  number  from  1  to  1024.  If  the  user  specifies  a  non-zero 
offset,  the  number  of  points  that  can  be  acquired  from  that  channel  Is  1024 
minus  this  offset. 

Finally  the  user  may  plot  the  most  recently  acquired  waveform.  This 
option  Is  straightforward  with  one  exception:  the  user  must  specify  the 
vertical  axis  limits  (YMIN,  YMAX).  If  the  user  enters  these  two  values  as 
equal,  the  routine  will  search  the  data  for  actual  minimum  and  maximum,  and 
auto-scale  the  plot  to  represent  the  entire  range  of  data. 
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Table  4.5 

Range  and  Interpretation  of  Octal  Data  Codes  for  Host  Commands 


Functions 

Range  of  Octal 
Data  Codes 

Interpretation 

Trigger  Delay  Coarse 

377g-0 

0  ns  to  255  *  55  ns. 

Delay  In  55  ns  increments. 

Trigger  Delay  Fine 
(2  channels  merged) 

0-3778 

The  four  LSB's  correspond  to 
the  low  channel  and  the  four 

MSB's  correspond  to  the  high 
channel.  For  either  channel 
the  range  is  0  ns  to  15  *  4  ns. 
Delay  In  4  ns  Increments. 

Receiver  Voltage  Gain 
(2  channels  merged) 

0-37 7g 

The  four  LBS's  correspond  to 
the  low  channel  and  the  four 

MSB's  correspond  to  the  high 
channel.  For  either  channel 
the  range  is  -45  dB  to  0  dB. 

Gain  In  3  dB  Increments. 

Transmit  Group 

00 

o 

o 

o 

Channels  1-16 

and 

oo 

o 

H 

o 

Channels  17-32 

Receive  Group 

020g 

Channels  33-48 

o 

CO 

o 

00 

Channels  49-64 

Conversion  Delay 

7777g-0 

0  ns  to  4095  x  55  ns. 

Delay  in  55  ns  Increments. 

4.2.7  Array  Processor  Control  and  Operation 

The  array  processor  manufactured  by  Analogic  and  designated  the  AP400 
Is  used  for  the  high  speed  acquisition  of  data  from  the  ADE  and  serves  as  the 
"beam-former"  for  the  phased  array  system.  The  AP400  communicates  with  the 
S/200  minicomputer  or  host  via  the  calling  routine  ARRAY  described  In  Section 
4.2.6.  This  calling  routlre  controls  the  sequence  of  events  occurring  In  the 
AP400  through  a  set  of  subroutines  known  as  "K-calls".  Most  of  these  assembly 
language  subroutines  were  provided  by  Analogic  In  the  form  of  a  library  although 
It  was  necessary  to  modify  some  routines  to  perform  specific  activities.  These 
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“K-calls"  control  the  operations  performed  by  the  AP400  through  a  set  of  AP400 
resident  subroutines  known  as  "Q- routines".  The  "Q-routlnes"  were  also  provided 
by  Analogic. 

Although  the  array  processor  can  be  used  In  an  Interactive  program 
which  allows  the  user  to  specify.  In  order,  each  function  that  Is  to  be 
executed,  the  unit  Is  used  with  ARRAY  In  a  mode  that  executes  a  set  of  commands 

specified  by  the  "PR"  section  of  ARRAY.  This  module  acquires  the  data  from  each 

of  the  16  channels  of  the  ADE  and  combines  them  In  a  specific  way  to  synthesize 

a  beam  returning  from  a  specified  angle.  This  synthesis  is  accomplished  by 
shifting  each  of  the  waveforms  by  a  time  delay  that  will  permit  the  waves  to  sum 
coherently  in  the  direction  specified.  For  a  sector  scan,  this  direction 
corresponds  to  that  of  the  transmitted  beam.  As  each  waveform  Is  shifted,  it  is 
added  to  a  running  sum  until  all  16  waveform  have  been  combined.  There  are 
several  ways  available  to  shift  the  waveforms.  Perhaps  the  most  straightforward 
technique  Is  to  shift  the  waveform  In  the  time  domain  by  an  Integral  number  of 
sample  intervals.  The  sample  Interval  for  the  ADE  is  55  ns.  When  a  beam  is 
formed  with  this  coarse  quantization  of  the  delay,  the  side  lobe  level  becomes 
considerably  higher  than  desired.  Figure  4.21a  shows  a  beam  that  has  been 
formed  from  an  apodlzed  2.5  MHz  array  of  16  elements  with  no  quantization  of  the 
time  delays  associated  with  the  elements.  When  the  same  array  Is  used  to  form  a 
beam  with  a  quantization  of  50  ns,  the  beam  shown  In  Fig.  4.21b  results.  The 
side  lobe  amplitudes  have  Increased  by  approximately  40  dB. 

As  the  quantization  of  the  time  delay  Is  decreased,  the  side  lobe 
amplitudes  are  reduced.  The  best  way  to  decrease  the  quantization,  which 
Implies  some  sort  of  Interpolation  between  samples,  Is  to  perform  the  time 
shifting  in  the  frequency  domain.  The  signals  are  Fourier  transformed  and  then 
multiplied  by  the  complex  phase  factor  exp  [lux]  where  x  Is  the  desired  time 
shift.  The  signal  Is  then  Inverse  Fourier  transformed  and  summed  to  synthesize 
the  beam  In  the  desired  direction. 

The  array  processor  can  do  the  Fourier  transforming  operations  In  about 
the  same  time  as  the  simple  time  shifting  operations.  Thus  the  unit  permits  the 
rapid  synthesis  of  a  beam  with  low  sldelobes.  The  current  algorithm  completes 
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Fig.  4.21  Ultrasonic  beam  formed  by  a  2.5  MHz  array  with  Hanning  amplitude 
weighting,  (a)  No  time  delay  quantization,  (b)  50  ns  time  delay 
quantization. 
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the  acquisition  and  manipulation  of  16  waveforms,  each  with  512  sample  Inter¬ 
vals,  in  400  msec.  Thus  for  a  sector  scan  containing  100  angles,  40  seconds  are 
required  to  process  the  data  exclusive  of  the  time  necessary  to  generate  a 
display.  The  manner  In  which  this  time  Is  partitioned  among  the  various 
computational  activities  Is  shown  In  Fig.  4.22.  It  Is  believed  that  this  time 
can  be  significantly  reduced  by  shifting  all  of  the  computational  activities  to 
the  AP400  (which  requires  additional  data  memory  in  the  AP400)  and  overlaying 
the  disk  storage  activities.  The  processing  time  for  one  angle  can  then  be 
reduced  to  80  msec,  l.e.,  by  a  factor  of  five. 

The  operation  of  the  AP400  under  the  control  of  the  S/200's  RD0S 
operating  system  had  not  been  tried  before  ARRAY  was  written.  This  has 
permitted  a  much  better  understanding  of  the  limitations  of  the  AP400  with 
regard  to: 

1.  Required  size  of  program  memory 

2.  Required  size  of  data  memory 

3.  Ordering  of  calls  to  perform  specific  functions. 

There  are  Improvements  that  can  be  made  in  the  operation  of  this  unit  but  they 
have  only  been  realized  as  a  result  of  having  designed  the  first  system. 


4.2.8  01  splay  of  Data 

The  general  features  of  the  Genlsco  display  processor  were  described  In 
Section  2.3.  Recall  that  the  unit  contains  a  480  x  512  pixel  memory  which  Is 
mapped  onto  the  3:4  aspect  ratio  CRT  of  a  television  monitor.  The  memory  Is 
filled  with  data  derived  from  data  on  the  disk  memory  of  the  S/200  mini¬ 
computer.  The  method  used  to  store  the  data  in  the  memory  of  the  display 
processor,  l.e.,  the  angles  and  positions  of  the  vectors  corresponding  to  the 
ultrasonic  waveforms  Is  controlled  by  a  software  routine  that  Is  specific  to  the 
application.  Each  of  the  pixels  can  be  characterized  by  an  8  bit  word  to 
provide  a  level  of  gray  or  a  color  according  to  a  code  that  Is  chosen  by  the 
operator. 


The  scan  pattern  of  the  phased  array  transducer  Is  in  general  some  sort 
of  a  compound  B-scan.  For  the  contour  following  mode  the  pattern  of  scan  lines 
can  become  complex  with  some  lines  crossing  over  others.  The  task  of  displaying 
a  B-scan  of  this  sort  has  not  been  pursued  because  the  32  element  arrays  are  not 
long  enough  to  readily  lend  themselves  to  this  mode  of  operation.  For  a 
detailed  examination  of  a  flaw  using  imaging,  it  is  appropriate  to  use  a  sector 
scan.  An  example  of  this  sort  of  scanning  pattern  is  illustrated  in  Fig. 

4.23.  The  data  obtained  with  the  program  ARRAY  can  be  displayed  in  this  format 
by  a  program  called  WEDGE.  The  software  for  operating  the  array  allows  the 
operator  to  specify  the  maximum  and  minimum  angles  of  the  sector,  the  number  of 
lines  within  the  sector,  the  number  of  degrees  in  the  sector,  the  starting 
angle,  and  the  water  path  distance  before  the  waveform  digitization  begins. 

This  in  turn  implies  that  the  size  and  shape  of  the  displayed  segment  of  the 
sector  will  change  as  the  input  variables  are  changed.  To  insure  that  all  of 
the  data  is  available  for  display,  the  sector  segment  is  centered  on  the  TV 
screen  and  scaled  so  that  there  is  a  guard  band  corresponding  to  4%  of  the 
screen  height  at  the  top  and  bottom  of  the  display  as  indicated  in  Fig.  4.23. 

For  certain  sets  of  input  variables  this  scaling  technique  does  not  fill  the  CRT 
very  well,  so  there  is  a  provision  for  the  operator* to  change  the  scale  factor 
and  thus  magnify  the  image. 

With  metal  specimens  there  is  a  significant  refraction  of  the  acoustic 
beam  as  it  enters  the  metal.  There  is  the  possibility  that  the  beam  will  not 
enter  at  all  because  it  is  incident  at  an  angle  larger  than  the  critical 
angle.  Once  inside  the  metal  the  acoustic  waves  propagate  faster  than  in  the 
water.  So,  equal  time  Increments  of  the  waveform  will  correspond  to  different 
spatial  distances  than  they  did  In  the  water.  The  complications  that  result 
from  the  refraction  caused  by  a  4  to  1  velocity  difference  are  Illustrated  In 
Fig.  4.24  for  a  sample  that  has  a  flat  front  face  and  a  fan  beam  scanning  from 
-10  degrees  to  +10  degrees.  Note  that  the  Internal  angles  are  -44  degrees  to 
+44  degrees.  A  considerable  distortion  of  the  Image  will  result  if  these 
Internal  angles  are  not  properly  taken  Into  account  when  the  data  Is  dis¬ 
played.  In  principle  this  sort  of  correction  presents  no  difficulty,  although 
In  practice  one  must  accurately  know  the  position  of  the  transducer  relative 
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to  the  surface.  If  the  surface  has  some  curvature  this  also  has  to  be  known. 

For  an  inspection  task  on  many  identical  parts  appropriate  jigs  can  be  fabri¬ 
cated  to  orient  the  part  and  the  curvature  of  the  surface  can  be  programmed  into 
the  display.  One  of  the  advantages  of  a  phased  array  is  the  ability  to  adjust 
the  phases  of  the  elements  to  compensate  for  the  focusing  or  defocusing  at  a 
curved  metal  interface.  To  properly  display  the  refracted  beam,  the  operator 
enters  the  ratio  of  the  velocity  of  the  metal  to  that  of  the  water  and,  for  a 
flat  surface,  the  perpendicular  distance  between  the  transducer  and  the  part. 

The  rays  are  then  displayed  as  being  bent  at  the  water  metal  interface  as  shown 
in  Fig.  4.24.  To  account  for  the  change  in  the  spatial  dimension  along  the 
beam,  the  length  of  each  displayed  ray  is  stretched  out  by  a  factor  equal  to  the 
ratio  of  the  metal  and  water  velocities. 

4.2.9  Packaging  of  Array  Subsystems 

The  components  for  the  phased  array  system  are  packaged  in  several 
separate  groups.  The  minicomputer,  array  processor,  and  display  processor  stand 
by  themselves.  There  are  two  custom  assemblies.  A  relay  rack  which  houses  a 
card  cage  that  contains  the  16  boards  for  the  receiver  channels,  the  boards  for 

the  timing  and  control  circuitry,  and  the  boards  for  the  control  memory.  Also 

housed  in  the  relay  rack  are  the  power  supplies  for  most  of  the  system.  The 
second  custom  assembly  contains  the  multiplexers,  the  hybrid  pulser/  receivers, 
and  the  transducers.  A  photograph  of  this  unit  Is  shown  in  Fig.  4.25.  The 
circuit  boards  are  mounted  In  a  waterproof  box  made  from  acrylic  sheet  stock. 

The  box  is  attached  to  a  frame  mounted  on  the  scanning  bridge  of  the  Test  Bed. 

The  frame  rides  on  roller  bearings  and  is  mechanically  linked  to  the  X-axis 

drive  of  the  scanning  system.  This  allows  It  to  be  moved  laterally  in  X  and  Y 
under  computer  control.  The  box  can  be  manually  raised  and  lowered  using  a  rack 
and  pinion  drive.  The  spacing  between  the  two  transducer  array  assemblies  can 
be  changed  manually  through  a  threaded  drive  assembly.  The  precision  Is  such 
that  one  complete  turn  o*  the  adjusting  knob  corresponds  to  a  change  in  spacing 
of  0.100  Inches.  The  spacing  can  be  varied  from  zero  Inches  (i.e.,  the  array 
transducer  housings  are  physically  touching)  to  a  maximum  of  6  inches.  The 
various  assemblies  are  Interconnected  via  a  combination  of  ribbon  cables  and 
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Fig.  4.25  Multiplexer,  pul ser/recei ver  and  transducer  assembly  mounted 
on  Test  Bed  scanning  bridge. 
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coaxial  cables.  The  ribbon  cables  are  used  for  power  supply  voltages  and 
digital  signals.  The  coaxial  cables  have  a  characteristic  impedance  of  75  ohms 
and  are  used  for  transmitting  the  16  analog  received  signals  from  the  analog 
line  drivers  in  the  multiplexer  box  to  the  signal  conditioning  circuitry. 

4.3  Operation  of  Phased  Array  System 

4.3.1  Array  Radiation  Pattern  Characteristics 

The  design  of  the  transducer  array  determines  the  angular  radiation 
pattern  produced  by  the  system.  Of  particular  interest  are  the  characteristics 
of  the  grating  lobes  and  side  lobes  generated  by  the  array.  Recall  that  the 
radiation  pattern  or  beam  pattern  of  an  array  is  obtained  by  summing  the 
contributions  of  each  of  the  elements.  In  the  far  field  the  sum  will  be  the 
product  of  two  factors.  One  factor  will  be  the  radiation  pattern  of  the 
individual  element  and  the  other  will  be  the  radiation  pattern  of  an  array  of 
point  transducers.  The  diagram  in  Fig.  4.26a  illustrates  the  distinction 
between  the  two  factors.  The  array  factor  can  become  large  in  several  direc¬ 
tions  other  than  the  direction  of  the  main  beam.  These  are  termed  grating  lobes 
and  the  angles  at  which  they  occur  relative  to  the  main  beam  are  a  function  of 
the  spacing  between  elements.  If  the  phases  of  the  signals  applied  to  the 
elements  are  adjusted  so  that  the  main  beam  Is  deflected  through  an  angle  as 
Illustrated  in  Fig.  4.26b,  then  the  grating  lobes  are  deflected  through  the  same 
angle.  If  the  array  factor  Is  considered  alone,  the  amplitudes  of  the  grating 
lobes  are  equal  to  that  of  the  main  beam.  The  single  element  factor  will  modify 
this,  but  It  is  Important  to  realize  that  the  amplitude  of  the  grating  lobes  can 
become  equal  to  or  larger  than  that  of  the  main  beam  for  sufficiently  large 
deflection  angles.  For  a  periodic  array,  the  only  ways  to  suppress  the 
amplitude  of  the  grating  lobes  are  to  make  the  spacing  between  the  elements  less 
than  a  half  wavelength,  measured  In  the  propagating  medium,  at  the  highest 
frequency  of  Interest,  or  to  use  short  pulses  to  excite  the  array. 

The  characteristics  of  the  side  lobes  associated  with  the  main  lobe 
(and  replicated  on  each  of  the  grating  lobes)  are  determined  by  the  overall 
width  of  the  array  as  well  as  the  amplitude  and  phase  distribution  applied  to 
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the  elements.  For  the  Test  Bed  array,  the  width  is  fixed  by  the  choice  of  using 
16  contiguous  transducer  elements  at  a  time.  The  phase  and  amplitude  distribu¬ 
tions  can  be  varied  and  Fig.  4.21  is  an  example  of  changes  that  can  be  effected 
through  altering  the  phase  distribution.  An  example  of  the  radiation  pattern 
from  the  Test  Bed  array  is  shown  in  Fig.  4.27  for  a  uniform  amplitude  distribu¬ 
tion.  This  should  be  compared  to  Fig.  4.21a  where  Hanning  amplitude  weighting 
is  used. 

The  angular  scanning  range  of  an  array  Is  determined  by  the  radiation 
pattern  of  the  individual  element  as  indicated  in  Fig.  4.26.  The  lateral  extent 
of  this  radiation  pattern  is  chiefly  determined  by  the  width  of  the  element. 

For  NOE  applications  involving  flaws  in  metals,  the  critical  angle  at  the  water/ 
metal  interface  is  about  15  degrees.  The  array  scanning  range  need  not  be  much 
larger  than  this  if  only  parts  with  flat  surfaces  are  to  be  examined.  If 
contour  following  is  to  be  considered,  the  range  of  angles  must  be  increased  to 
accommodate  the  tilting  of  the  front  surface  of  the  part.  The  Test  Bed  array 
transducer  can  scan  the  beam  over  ±  45  degrees  with  a  decrease  in  beam  amplitude 
of  about  3  dB.  Thus  the  tilting  of  surfaces  at  angles  as  great  as  30  degrees 
can  be  accommodated  and  still  permit  a  sector  scan  to  be  generated  that  would 
enter  the  metal  part. 

Although  the  timing  of  the  signals  transmitted  or  received  by  the  array 
can  be  adjusted  to  obtain  a  focused  beam,  this  is  of  no  value  because  of  the 
small  width  of  the  aperture  formed  by  16  contiguous  elements  (about  7.2  mm). 

The  far  field  distance  of  an  aperture  of  this  size  is  about  22  mm  in  water.  For 
focusing  to  be  effective,  the  working  distance  would  have  to  be  less  than  this 
distance.  It  should  be  noted  that  this  distance  scales  with  the  reciprocal  of 
wavelength  and  would  therefore  be  about  four  times  shorter  in  a  metal.  Signifi¬ 
cant  focusing  could  be  obtained  by  combining  signals  from  two  or  more  16  element 
groups  spaced  at  distances  of  up  to  150  mm  at  the  expense  of  decreased  angular 
separation  between  the  main  lobe  and  the  grating  lobes.  This  could  well  be  a 
useful  approach  In  situations  where  It  has  been  ascertained  that  the  grating 
lobes  will  not  hit  any  flaws.  It  would  effectively  be  a  detailed  examination  In 
which  the  resolution  was  being  Increased  via  a  synthetic  aperture  technique. 
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Fig.  4.27  Array  radiation  pattern  for  uniform  excitation  of  16  elements 
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The  angular  width  of  the  beam  produced  by  the  16  element  arrays  Is 
essentially  the  same  as  that  produced  by  a  single  element  transducer  having  the 
same  aperture.  For  a  7.2  mm  aperture,  this  can  become  rather  large  as  indicated 
in  Fig.  4.28.  To  show  that  this  lack  of  resolution  Is  the  same  for  both  an 
array  transducer  and  a  single  element  transducer.  Images  of  the  same  1200  un 
flaw  were  made  and  are  compared  in  the  next  section  after  the  details  of  the 
imaging  technique  have  been  explained. 

4.3.2  Imaging  of  Flaws  in  Metals 

The  concerns  associated  with  grating  lobes,  side  lobes,  etc.  relate  to 
the  ability  of  the  array  to  distinguish  between  targets  in  the  lateral  direc¬ 
tion.  For  the  imaging  of  flaws  in  metals  the  loss  of  longitudinal  resolution 
due  to  the  Increase  in  velocity  Imposes  severe  constraints  on  the  use  of  rela¬ 
tively  low  frequency  arrays.  Particularly  troublesome  is  the  low  level  of  the 
flaw  signal  relative  to  the  front  face  echo.  In  the  case  of  a  1200  micrometer 
diameter  spherical  void,  the  front  face  signal  must  decay  to  a  level  that  is  40 
dB  below  Its  peak  value  before  the  flaw  signal  Is  detectable  without  some  sort 
of  post  processing.  Although  the  transient  response  measurements  on  the  array 

in  Fig.  4.5  do  not  extend  to  -40  dB,  It  can  be  seen  that  even  at  -30  dB  the 

elapsed  time  is  4  us  corresponding  to  a  flaw  depth  of  one-half  Inch  In  a  metal 
like  titanium.  This  severely  restricts  the  use  of  real  time  arrays  for 
detecting  flaws  In  metals  except  in  cases  where  the  flaw  signals  are  very  strong 
or  the  flaws  are  located  one  Inch  or  more  below  the  surface. 

One  possibility  for  minimizing  the  front  surface  echo  Is  to  propagate 

through  It  at  an  oblique  angle.  The  extent  to  which  this  is  useful  depends  on 

the  magnitude  of  the  angle,  whether  the  part  can  be  tilted  relative  to  the 
array,  and  the  strength  of  the  target.  It  has  been  noted  that  the  reflection 
from  a  flat  part  remains  reasonably  strong  over  the  entire  range  of  ±10  degrees 
that  has  typically  been  used  for  sector  scans.  This  results  in  an  Imaging 
artifact  that  appears  as  an  arc  of  a  circle  that  Is  tangent  to  the  surface  of 
the  part  at  the  point  where  the  angle  of  Incidence  of  the  bean  is  zero  degrees. 
Since  the  amplitude  of  the  backscatter  from  the  oblique  Intersection  of  the  beam 
with  the  flat  surface  Is  very  small  relative  to  a  specular  reflection.  It  Is 
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Fig.  4.28  Polar  plot  of  radiation  pattern  within  a  metal  part  when  there  Is 
50  m  water  path. 


usually  not  visible  in  the  Image.  Thus  a  typical  image  obtained  from  a  part 
appears  like  the  one  shown  in  Fig.  4.29. 


As  it  stands  this  image  is  not  very  useful.  The  “ring  down"  from  the 
front  surface  echo  is  obscuring  the  flaw  that  is  present  and  the  artifact 
created  by  the  reflection  from  the  front  surface  is  creating  the  illusion  that 
the  surface  is  curved.  Most  of  these  problems  can  be  overcome  by  using  a 
differential  imaging  technique.  This  means  that  a  reference  image  is  obtained 
of  a  region  of  the  part  that  contains  no  flaw  but  has  essentially  the  same 
surface  reflection  characteristics  as  the  region  containing  the  flaw.  The 
reference  image  is  then  subtracted  from  the  flaw  image  to  eliminate  the  problems 
that  are  created  by  the  front  surface  echo.  When  this  technique  is  used  with 
the  image  in  Fig.  4.29,  the  differential  image  shown  in  Fig.  4.30  is  obtained. 
The  subtraction  has  suppressed  the  front  surface  echoes  by  about  25  dB  and  the 
image  of  the  1200  pm  spherical  void  stands  out  clearly. 

The  lateral  extent  of  the  flaw  reflects  the  width  of  the  ultrasonic 
beam  that  results  from  a  7.2  mm  wide  transducer  when  it  propagates  through  a  50 
mm  water  path  and  a  12.5  mm  titanium  path.  (Some  idea  of  the  scale  in  Fig.  4.30 
can  be  obtained  from  the  fact  that  the  distance  between  the  front  and  back 
surface  of  the  titanium  Is  25  ran).  This  does  not  result  from  any  character¬ 
istics  of  the  array  other  than  its  aperture.  Figure  4.31  shows  a  B-scan  image 
of  the  same  flaw  that  was  obtained  using  a  differential  imaging  technique  with  a 
single  element  transducer  whose  aperture  was  7.2  mm.  The  streaks  which 
demonstrate  the  difficulty  In  obtaining  a  differential  Image  when  the  transducer 
has  to  be  mechanically  moved,  are  caused  by  slight  Irregularities  in  the  motion 
of  the  transducer.  Note  that  the  width  of  the  flaw  Is  about  the  same  as  In  Fig. 
4.30. 


4.3.3  Born  Inversion  of  Array  Data 

A  feature  of  the  array  electronics  was  to  process  the  received  signals 
with  sufficiently  low  distortion  that  the  array  could  be  used  for  acquiring 
scattering  data.  To  show  that  this  Is  feasible,  the  zero  degree  waveform  that 
makes  up  one  line  of  the  Image  In  Fig.  4.30  was  used  as  an  Input  signal 
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Fig.  4.29  Unprocessed  image  obtained  from  a  titanium  part  with  a  flat 
front  surface. 
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Fig.  4.31  B-scan  of  1200  nm  spherical  void  obtained  with  a  single 
element  transducer  having  the  same  aperture  as  the  array. 
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for  the  Born  Inversion.  An  unsubtracted  backsurface  echo  was  used  as  a 
reference  waveform  for  performing  the  deconvolution.  The  bandwidth  was  centered 
at  about  ka  =  1.2  and  the  ratio  of  [ka]max  to  [ka]min  was  about  3:1.  The  center 
of  the  flaw  was  determined  by  extrapolating  the  phase  at  low  frequencies.  The 
Born  inversion  estimate  of  the  radius  was  520  urn  which  is  low  by  13%. 

Considering  that  this  is  the  first  Born  inversion  data  that  has  been  obtained 
with  an  array,  the  agreement  is  excellent. 

In  principal  the  array  provides  a  convenient  technique  for  acquiring 
scattering  data  as  a  function  of  angle,  although  time  has  not  permitted  this  to 
be  done  with  the  1200  urn  void.  Despite  this  potential  advantage  it  is  unlikely 
that  an  array  transducer  of  the  type  described  here  will  be  useful  for  the 
detailed  characterization  of  flaws  via  the  Born  inversion  algorithm.  The  Born 
inversion  work  described  in  Section  3.2.2  established  the  bandwidth  requirements 
of  this  algorithm  as  being  from  [ka]min  s  0.75  to  [ka]max  2  1.8.  From  this  it  is 
apparent  that  the  array  transducers  used  here  are  only  adequate  for  a  narrow 
range  of  flaw  sizes.  To  measure  an  unknown  flaw  a  selection  of  array 
transducers  whose  collective  bandwldths  spanned  a  wide  range  of  frequencies 
would  be  necessary.  The  array  approach  Is  much  more  likely  to  be  useful  in  an 
imaging  role  where  various  synthetic  aperture  techniques  can  be  used  to  increase 
the  spatial  resolution. 


5.0  SUMMARY 


The  Ultrasonic  Test  Bed  Program  has  been  an  ambitious  program  that  has 
had  to  cover  a  wide  range  of  technological  arenas.  On  the  one  hand  there  is  the 
hardware  intensive  multiaxis  system  that  has  provided  the  scanning  capability. 
This  is  coupled  with  the  data  acquisition  capability  that  permits  the 
digitization  of  entire  waveforms.  These  units  along  with  a  color  display 
processor  work  interactively  under  the  control  of  a  minicomputer.  Most,  if  not 
all,  of  these  units  are  commercially  available;  although  the  total  package 
equals  or  exceeds  the  capabilities  of  most  conventional  NDE  systems. 

Using  this  system  it  is  possible  to  study  the  utility  and  practica¬ 
bility  of  the  inversion  techniques  that  have  arisen  from  the  Interdisciplinary 
Quantitative  NDE  program.1  Here  the  object  of  the  work  shifts  to  translating 
theoretical  and  experimental  tools  Into  techniques  that  are  useful  in  practical 
NDE  programs.  For  the  first  time  an  integrated  effort  has  been  made  to  combine 
the  bits  and  pieces  of  information  from  a  flaw  that  are  available  over  a  very 
wide  frequency  range.  The  object  Is  to  obtain  a  significantly  better 
characterization  of  a  flaw  then  has  heretofore  been  possible  when  only  straight¬ 
forward  measurements  of  the  signal  amplitude  produced  by  the  flaw  were  used. 

Beyond  this  the  Test  Bed  program  has  produced  a  flexible  state-of-the- 
art  phased  array  system  that  can  be  used  to  determine  the  feasibility  of  array 
systems  In  NDE  applications.  The  system  represents  a  hardware  intensive  project 
that  Involves  new  ultrasonic  array  technology  at  one  end  and  digital  array 
processor  technology  at  the  other.  In  between  are  16  custom  A/D  converters  and 
64  hybrid  pul ser/recel vers.  The  synthesis  of  this  system  requires  many  hardware 
and  software  skills  to  be  focused  on  techniques  and  systems  that  have  just 
become  available. 

The  Science  Center  has  been  able  to  bring  together  In  a  relatively 
small  team  of  people  the  skills  required  for  this  technologically  diverse 
project.  The  major  goals  of  the  program  have  been  met.  However,  as  with  any 
large  program  there  are  areas  that  will  require  more  work.  Further,  there  are 
some  things  that  would  be  done  differently  in  a  second  generation  system.  The 


discussion  below  summarizes  the  work  that  has  been  accomplished,  points  out 
technological  needs  for  improved  capabilities,  indicates  desirable  design 
changes,  and  discusses  areas  requiring  additional  work. 

Recall  that  the  major  hardware  subsystems  associated  with  the  Test  Bed 
are  the  minicomputer,  multiaxis  controller,  the  data  acquisition  system  and  the 
display  processor.  Two  of  these  units,  the  multi  axis  controller  and  the  A/D 
converter  portion  of  the  data  acquisition  system,  represent  areas  of  rapidly 
changing  technology.  The  discussion  in  Section  2.0  of  the  report  covered  many 
of  the  problems  encountered  with  these  units.  It  should  be  noted  that  a  much 
better  mechanical  system  could  be  realized  by  placing  each  of  the  axes  in  a 
closed  loop  control  system  using  encoders  linked  as  closely  as  possible  to  the 
transducer  motion.  This  permits  a  rapid  assessment  of  the  location  of  the 
transducer,  permits  more  accurate  positioning  without  the  need  of  excessively 
rigid  mechanical  structures,  and  provides  signals  that  are  directly  linked  to 
the  motion  of  the  transducer  for  triggering  the  pulser/recelver.  The  multiaxis 
controller  suffers  from  not  having  a  more  powerful  diagnostic  capability  that 
would  make  it  friendlier  to  the  operator.  Robustness  in  the  operation  of  this 
unit  would  be  a  desirable  goal.  In  a  system  depending  on  several  asynchronous 
digital  interfaces,  it  Is  desirable  to  link  the  timing  of  the  data  acquisition 
system  tightly  to  the  motion  of  the  transducer  to  ensure  that  the  area  scanned 
Is  sampled  uniformly.  This  objective  requires  some  design  changes  in  the 
present  system  but  could  be  readily  achieved  In  a  second  generation  system. 

The  technology  of  high  speed  A/D  converters  Is  rapidly  changing  with 
the  emphasis  on  both  higher  sampling  speeds  and  higher  accuracy  of  the  sampled 
data.  For  most  current  NDE  needs  the  sampling  speed  is  adequate  but  the 
accuracy  needs  to  be  Improved.  It  would  be  desirable  to  have  an  accuracy 
equivalent  to  eight  bits  at  a  data  frequency  of  25  MHz  without  the  need  to 
temporally  average  signals.  Although  newer  units  that  have  become  available 
commercially  are  more  accurate  they  still  fall  short  of  this  goal.  The  digital 
data  acquisition  device  used  with  an  NDE  system  should  also  have  a  high  speed 
bus  capability  to  convey  data  to  a  host  computer  at  a  high  bit  rate.  A  rate  of 
500,  1024  point  waveforms  per  second  would  satisfy  many  needs. 
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Section  3.0  on  inspection  procedures  covered  many  topics.  The  first  of 
these  were  concerned  with  the  search  mode.  The  Test  Bed  Program  has  not  devoted 
too  many  resources  to  the  detectability  problem  nor  to  the  reliability  of 
detection.  This  is  recognized  as  an  important  area  of  study  which  should  be 
pursued  with  the  same  diligence  as  the  techniques  for  the  characterization  of 
the  flaw.  There  are  a  number  of  generic  detection  techniques  that  can  be 
included  in  future  work  as  well  as  specific  problems  that  relate  to  particular 
geometries.  The  Test  Bed  system  uses  a  straightforward  amplitude  threshold 
technique  to  detect  signals  of  interest.  The  reliability  can  be  increased  by 
making  the  threshold  a  function  of  depth  within  the  metal.  Future  work  should 
be  addressed  towards  subsequent  analysis  of  the  detected  segments  of  waveforms 
to  compare  them  with  known  properties  of  flaw  signals.  In  this  way  the 
reliability  of  detection  can  be  Improved. 

The  sensitivity  of  the  detection  techniques  can  be  improved  by  simply 
using  a  focused  transducer  rather  than  an  unfocused  one.  To  the  extent  that  the 
flaw  size  is  less  than  the  focal  spot  size,  each  reduction  in  the  focal  spot 
size  by  a  factor  of  two  will  Increase  the  flaw  signal-to-noise  ratio  by  a  factor 
of  16.  This  implies  an  increase  of  the  scanning  time  by  a  factor  of  four  to 
cover  the  same  area.  Detailed  studies  of  these  Important  areas  of  detection 
sensitivity,  scanning  time,  and  reliability  of  detection  have  yet  to  be  done, 
but  the  Inherent  scanning  and  signal  processing  capabilities  of  the  Test  Bed 
system  will  permit  these  studies  to  be  carried  out  efficiently. 

The  first  aspect  of  the  detailed  examination  of  a  flaw  that  was 
considered  was  the  Imaging  techniques.  The  techniques  used  were  C-s can  and  B- 
scan  techniques  which  utilized  the  display  processor  for  presenting  a  gray  scale 
or  color  coded  Image.  Future  work  could  Improve  these  capabilities  In  two 
areas.  First  It  would  be  desirable  to  display  a  larger  dynamic  range  In  the 
Image  then  Is  currently  possible.  The  limit  Is  currently  set  by  the  available 
dynamic  range  of  the  A/0  converter  of  about  48  dB.  This  can  be  extended  by 
using  nonlinear  processing  of  the  signals  prior  to  digitization  so  that  the 
dynamic  range  of  the  signals  Is  compressed  into  the  48  dB  range  of  the  A/D 
converter.  After  digitization  the  signals  must  either  be  compressed  further  or 
color  coded  to  fit  Into  the  approximate  30  dB  dynamic  range  of  the  display 


monitor.  The  use  of  these  techniques  would  permit  a  viewer  to  simultaneously 
observe  signals  produced  by  both  the  flaws  and  the  grain  scattering.  It  would 
also  permit  the  color  coding  capability  of  the  system  to  be  used  to  full 
advantage. 

The  second  area  that  could  be  pursued  is  Improved  resolution  of  the 
focused  transducers.  At  the  present  time  commercially  available  transducers  are 
used  that  have  relatively  low  numerical  apertures  (N.A.)  and  consequently 
produce  focal  spot  sizes  that  are  considerably  larger  than  the  theoretical 
minimum  of  one  wavelength.  When  the  beam  is  focused  inside  a  metal  part  they 
suffer  additional  degradation  due  to  spherical  aberration  at  the  metal /water 
interface.  Performance  in  this  area  can  be  improved  with  specially  designed 
focusing  lenses.  It  should  be  noted  that  the  inspection  protocol  assumes  that 
imaging  can  be  used  for  k a  values  as  low  as  6,  but  the  capabilities  of 
commercially  available  transducers  imposed  a  minimum  ka  value  of  between  20  and 
30.  Therefore  the  resolution  capabilities  of  a  QNDE  system  should  be  improved 
to  permit  the  minimum  imaging  performance  requirements  to  be  met. 

The  next  segment  of  the  detailed  examination  procedure  was  concerned 
with  the  use  of  the  Born  inversion  and  the  long  wavelength  techniques.  The  Born 
inversion  was  demonstrated  to  be  useful  for  the  characterization  of  flaws  in  the 
range  of  ka  values  that  were  beyond  the  capabilities  of  imaging  techniques.  It 
worked  both  for  known  flaws  of  ellipsoidal  shape  and  also  for  naturally  occurr¬ 
ing  flaws.  The  bandwidth  requirements  were  shown  to  be  within  the  capabilities 
of  available  transducers  and  the  accuracy  of  the  estimated  flaw  sizes  was  found 
to  be  robust  in  the  face  of  noise.  Techniques  were  demonstrated  for  determining 
the  location  of  the  center  of  the  flaw  using  long  wavelength  procedures.  The 
conclusion  is  that  the  Born  inversion  should  be  incorporated  into  the  procedure 
used  for  conventional  NOE  inspection.  Before  this  can  be  done  the  technique 
should  be  tested  on  a  large  group  of  representative  flaw  types  as  is  being  done 
in  the  Air  Force  program  to  develop  a  reliable  QNOE  module. ^  In  addition,  and 
this  is  perhaps  the  most  important  step,  the  procedures  for  applying  the  Born 
inversion  to  a  flaw  signal  must  be  automated.  This  will  permU  the  technique  to 
be  carried  out  with  minimal  input  from  the  operator.  Until  this  is  done  the 
technique  is  likely  to  remain  a  laboratory  tool. 
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With  the  advent  of  this  technique  and  the  long  wavelength  techniques, 
the  transducer  requirements  have  changed.  It  is  no  longer  essential  that 
different  but  nominally  identical  transducers  have  closely  matched  sensitivities 
or  bandshapes  since  these  characteristics  will  be  eliminated  from  the  flaw 
signals  during  the  deconvolution  procedure.  It  is  essential  that  the 
transducers  have  very  wide  bandwidths,  not  so  much  to  produce  short  transient 
responses  and  increase  longitudinal  resolution,  as  to  permit  the  analysis  of 
flaws  with  a  wide  range  of  diameters  and  to  obtain  long  wavelength  data.  Some 
work  has  been  done  on  the  problem  of  obtaining  signals  with  two  transducers 
having  different  but  overlapping  bandwidths;  and  then  after  deconvolution, 
splicing  the  two  specta  together.  This  has  worked  in  some  cases  but  requires 
additional  work  to  establish  a  robust  procedure  for  matching  the  two  spectra. 
Other  approaches  to  realizing  ultra  wide  bandwidth  transducers  require  tradeoffs 
of  sensitivity  but  may  be  viable  for  some  applications.  The  question  of 
transducer  requirements  that  is  raised  by  the  need  to  have  data  over  10:1  or 
20:1  bandwidths  needs  to  be  considered  in  a  comprehensive  way. 

After  the  tools  for  characterizing  flaws  have  been  applied  the  results 
are  combined  in  an  algorithm  that  uses  all  the  information  to  provide  a 
probabilistic  estimate  of  the  flaw's  character.  This  algorithm  was  used 
successfully  to  characterize  flaws  In  ceramic  materials.20  Although  it  is 
available  as  part  of  the  Test  Bed  software.  It  has  not  been  used  for  flaw 
characterization.  As  the  use  of  the  Born  inversion  and  long  wavelength 
techniques  becomes  better  established,  the  need  for  this  combining  algorithm 
will  increase.  At  the  present  time  the  data  is  combined  in  an  informal  way  by 
the  operator  to  arrive  at  a  complete  characterization  of  the  flaw,  an  example 
of  this  is  provided  by  the  characterization  of  the  naturally  occurring  flaw  that 
was  provided  by  Rolls-Royce  Ltd. 

The  last  part  of  the  protocol  is  to  subject  the  total  findings  about 
the  flaw  to  an  accept/reject  criteria.  This  criteria  Is  beyond  the  scope  of  the 
Test  Bed  Program  and  has  not  been  available  to  the  program.  The  establishment 
of  such  a  criteria  is  essential  to  a  QNDE  program  and  must  arise  from  the 
results  of  future  work. 


The  final  aspect  of  the  Test  Bed  Program  to  be  discussed  was  the  ultra¬ 
sonic  phased  array  system.  The  dominant  feature  of  this  system  Is  Its  flexi¬ 
bility.  The  system  utilizes  two  32  element  phased  array  transducers  with  a 
variable  spacing  between  them.  There  are  16  separate  receiver  channels;  each 
with  its  own  A/D  converter  and  buffer  memory.  The  beam  forming  for  the  array  Is 
done  by  an  array  processor  under  control  of  the  minicomputer.  The  scan  patterns 
can  be  changed  in  software  and  a  wide  variety  of  synthetic  array  techniques  can 
be  investigated. 

The  system  has  been  demonstrated  with  ellipsoidal  flaws  In  the  metal 
specimens  and  produces  images  with  the  same  characteristics  that  would  be 
obtained  with  a  scanned  single  element  transducer  having  the  same  aperture  and 
center  frequency.  The  quality  of  the  signals  from  the  array  Is  exceptionally 
good.  It  has  been  possible  to  process  these  signals  with  the  Born  inversion 
algoritm  and  obtain  an  accurate  estimate  of  the  flaw  size.  Since  the  system 
just  came  on  line  at  the  conclusion  of  the  Test  Bed  Program,  there  has  been  no 
opportunity  to  fully  exploit  it  and  exercise  all  of  its  modes  of  operation. 

The  immediate  plans  for  the  array  system  are  to  develop  software  for 
investigating  its  capabilities  in  a  synthetic  aperture  mode  of  operation  that 
will  permit  increased  lateral  resolution.  In  general  there  needs  to  be  a 
systematic  investigation  of  NDE  applications  that  might  be  addressed  more 
effectively  with  phased  arrays.  Since  this  system  is  coupled  to  a  minicomputer, 
extensive  post-processing  of  the  data  is  possible.  This  permits  practical 
difficulties  related  either  to  the  transient  response  of  the  elements  or  to  non 
uniformities  of  the  single  element  response  to  be  eliminated.  Thus  the  true 
capabilities  of  the  array  can  be  evaluated  and  system  needs  can  be  identified. 

In  conclusion  the  Test  Bed  Program  has  addressed  a  wide  range  of  topics 
that  are  Important  to  a  quantitative  NDE  program.  As  was  stated  In  the  work 
statement  for  this  program  it  will  remain  set  up  and  operating  for  at  least  one 
year  beyond  the  conclusion  of  the  program.  During  this  period  personnel  from 
outside  the  Science  Center  are  invited  to  use  the  system  to  address  specific 
problems. 
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APPENDIX  A 

DIFFUSION  BONDED  SAMPLE  DOCUMENTATION 


Documentation  for  Simulated  Turbine  Bore  Samples 

Included  In  this  Appendix  is  the  complete  documentation  for  the 
simulated  turbine  disk  bore  samples.  A  complete  listing  of  the  27  defects  and 
the  sample  that  contains  them  is  given  in  Table  Al.  In  Fig.  A1  the  three 
samples  are  shown  along  with  the  sites  of  the  defects.  In  Fig.  A2  through 
Fig.  A28,  the  detailed  documentation  for  each  of  the  defects  is  provided. 

Each  of  these  figures  is  laid  out  In  a  similar  manner.  The  sketch  in  the 
upper  left  corner  shows  the  configuration  of  the  sample  and  also  shows  which 
diffusion  bonding  plane  contains  the  defect.  The  drawing  on  the  upper  right 
shows  the  exact  location  of  the  defect  in  the  plane  with  respect  to  the  center 
line  of  the  sample  and  the  cyllndrically  curved  surface.  The  drawings  in  the 
middle  and  lower  left  show  cross  sectional  and  plan  views  of  the  defect.  The 
dimensions  given  are:  1)  the  nominal  defect  dimension,  2)  the  measured 
dimension,  and  3)  in  the  case  of  an  inclusion,  the  actual  dimensions  of  the 
inclusion.  The  photomicrograph  in  the  lower  right  shows  the  actual  plan  view 
of  the  defect  prior  to  diffusion  bonding. 
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List  of  Flaws  Contained  In  Simulated  Turbine  Bore  Samples 
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Fig.  A1  Simulated  turbine  bore  samples. 
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Sample  88-6:  400  um  WC  spherical  inclusion 


g.  A8  Sample  88-7:  1200  wm  simulated  fatigue  crack. 


Fig.  A10  Sample  88-9:  400  ym  AI0O3  spherical  Inclusion 


Fig.  All  Sample  89-1:  800  ym  AI0O3  spherical  Inclusion. 
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Fig.  A12  Sample  89-2:  100  urn  x  1200  ym  penny  shaped  crack 
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Fig.  A13  Sample  89-3:  800  gm  spherical  void 
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Fig.  A16  Sample  89-6:  1200  ym  spherical  void 
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Fig.  A17  Sample  89-7:  120C  pm  WC  spherical  inclusion. 
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Fig.  A19  Sample  89-9:  800  um  WC  spherical  inclusion 
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Fig.  A20  Sample  90-1:  200  um  x  800  um  oblate  spheroidal  void 
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Fig.  A23  Sample  90-4:  800  urn  x  100  um  prolate  spheroidal  void 
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Fig.  A24  Sample  90-5:  100  urn  x  1200  ym  penny  shaped  crack 
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Fig.  A25  Sample  90-6:  800  um  spherical  void 


1200  nm 
ACTUAL  1178  *im 


ACTUAL  1220  nm 
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APPENDIX  B 

DOCUMENTATION  FOR  PHASED  ARRAY  HARDWARE 


The  documentation  for  the  Phase  Array  acquisition  System  (PAAS)  is 
composed  of  system  drawings  and  electrical  schematics.  The  PAAS  Interconnect 
drawing  (Dwg  B-l)  shows  the  electrical  paths  between  the  main  units  of  the 
system. 

The  Array  Digitizer  Electronics  unit  (ADE)  is  housed  in  a  standard 
19-in.  relay  rack.  Two  card  cages  contain  the  system  control  electronics  and 
analog-to-digltal  conversion  circuits. 

Sixteen  digitizer  boards  handle  the  parallel  conversion  and  memory 
functions.  The  digitizer  electrical  schematic  (Dwg  B-2)  shows  the  details  of 
this  board.  The  control  functions  are  shared  on  two  boards.  A  ten-page 
timing  and  control  electrical  schematic  (Dwg  B-3)  contains  the  details  of 
these  boards. 

The  multiplexer  unit  (MUX)  Is  housed  in  a  custom  Plexiglas  enclosure 
mounted  on  the  bridge  of  the  scanning  tank.  The  MUX  interconnect  drawing  (Dwg 
B-4)  shows  the  electrical  connection  within  this  assembTy.  The  details  of  the 
buffer  control  board  are  shown  In  the  MUX  control  board  electrical  schematic 
(Dwg  B-5).  A  five-page  MUX  board  electrical  schematic  (Dwg  B-6)  displays  the 
details  of  the  multiplexer/demultiplexer  functions  contained  on  two  Identical 
boards.  The  hybrid  pulser-recelver,  Pr-SCR400V-1,  contains  two  complete 
circuits.  The  electrical  schematic  (Dwg  B-7)  and  electrical  parts  list  (Dwg 
B-8)  contain  the  details  for  one  circuit. 
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DRAWING  B-7 


TOLERANCE  INCH! 
DECIMAL  XX 
-XXX 


SURFACE 

ROUGHNESS 


UNLESS  OTHEI 


NO. 

NEXT  ASSY 

B40/41 


NEXT  ASST  I  USED  ON 


APPLICATION 


TOLERANCE:  INCHES,  METRIC 
DECIMAL  J(X  =t  Z'' 

joa  ±  / 


y  SURFACE 
ROUGHNESS 


UNLESS  OTHERWISE  NOTED: 


TEsr  BED 


PULSER  -  RECEIVER 

tiYBR.' d:  PR-  SCE4-00V-I 


WEIGHT  —  I  SHEET 


O  3  9  3  8  o  o  a  E  S  ooocSoooo 


PARTS  LIST  -  PULSER  RECEIVER 


TOLERANCES 

DCCIMAL:  «=-010 
»»»=  2  005” 
»«**  =  2:  .001” 
.  f II ACTIONAL:  2  14." 
'ACt:  y  I  ±15' 


PR-SCR400  V-1 


DWG  B-8  SHEET  1 


COMMENTS 


2N2907 

2N918 

2N42S8 

2N2907 

UPTB540 


400  V,  UNITRODE 


CHIP  #  UCT-0050.  $0.55 
OR  SIMILAR  REPLACEMENT 


C106D 

1N5818 

1N5818 

1N914 

1N914 


MOTOROLA  SCR,  400  V 
SBD,  MOT. 


10  A  PULSED,  500  NSEC 
10  A  PULSED,  500  NSEC 
10  A  PULSED,  500  NSEC 
1.3  A  PULSED,  500  NSEC 


0.01  pf.  +/-  10%,  50  V 

2000  pf.  +/-  10%,  500  V 

1E0  pf,  +/-  5%,  50  V 

0.033  pf,  +/-  5%.  50  V 

0.33  pf ,  +/-  20%,  50  V 

3300  pf  ,  +/-  20%,  50  V 

100  pf,  +/-  20%,  500V 

5.6  pf,  +/-  5%.  50  V 


;  10  A  PULSED  (DISCHARGE) 


NOTES: 

1.  C2  STABLE  OVER  TEMPERATURE  AND  VOLTAGE.  MUST  BE  LOW 
INDUCTANCE  AND  RESISTANCE. 

2.  D1  SCR  HAS  SPECIAL  SPECIFICATION:  IGT  MAX  -  200  pA,  IGT  MIN  -  200  pA 
FOR  TURN  ON  CONDITION. 

3.  SUBSTITUTION  CRITERIA:  05, 400  V;  D2,  D3,  PULSED  CURRENT  CAPABILITY. 
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Science  Cantor 
Rockwall  International 


■VT  R.B.  HOUSTON 


PARTS  LIST  -  PULSE R  RECEIVER 


10/2/79 


REF 


mm  tct  .  *> 


TOLERANCES 

OCCIMAL:  ««=~010 
««»=  ^  00b  ' 

*«**=  T  001' 
FRACTIONAL:  *  W*" 

:  ±J5' 


PR-SCR400V-1 


jjjjlT 


DWG  B-8  SHEET  2 


COMMENTS 


R1 

300  fl 

5%  ,  1/8  W  ,  50  V 

R2 

20  K  ft 

5%  ,  1/8  W  ,  50  V 

R3 

750  n 

5%  ,  1/8  W  ,  50  V 

R4 

150  SI 

5%  ,  1/8  W  ,  50  V 

R5 

12 

5% ,  1/8  W ,  50  V 

R6 

iooo  Kn 

5% .  1/4  W .  500  V 

R7 

1.8  K£2 

5%  .  1/2  W  ,  500  V 

(300  V  ,  300  NSEC  .  T  -  1  MSEC) 

R8 

1.8  KS2 

5%.  1/2  W  ,500  V 

(300  V ,  300 NSEC  ,  T=  1  MSEC) 

R9 

300  n 

5%.  1/2  W,  500  V 

(300  V  ,  300  NSEC  ,  T  =  1  MSEC) 

RIO 

120  SI 

5%  .  1/2  W  ,  500  V 

(300  V  ,  300  NSEC  ,  T  =  1  MSEC) 

R11 

5.1  Kfi 

2% ,  1/8  W ,  50  V 

R 12 

13  KS2 

2%,  1/2 W.  50V 

R13 

470  SI 

2% ,  1/8  W  ,  50V 

R14 

15S2 

1%  ,  1/8  W  ,  50  V 

R15 

1.2  K« 

2%  ,  1/8  W ,  50  V 

R16 

435  Si 

1% ,  1/8  W  ,  50  V 

R17 

220  n 

2% ,  1/8 W ,  50V 

R18 

4.7  KS2 

5%  ,  1/8  W ,  50  V 

R19 

4.7  KS2 

5% ,  1/8  W ,  50  V 

R20 

2.2  KS1 

5%  ,  1/8  W  ,  50  V 

NOTES: 

1.  R7-R10  OPERATE  IN  A  PULSED  MODE.  RESISTOR  SPECIFICATIONS 
MUST  BE  COMPATIBLE  WITH  OVERSTRESSED  CONDITION. 

2.  R14  AND  R16  NEED  TO  TRACK;  USE  SAME  MATERIAL,  PLACE 
ADJACENT. 
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